Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques - p 11 
Solution de l'équation de Laplace en coordonnées cartésiennes 2D à l'intérieur d'un domaine Q 


Solution de l'équation de Laplace à deux dimensions avec conditions aux limites sur un rectangle 
par la méthode de séparation des variables 


Sans réduire la généralité du problème nous pouvons positionner le rectangle à l'origine des 
coordonnées du plan cartésien (0,0), de dimension (k,1), où k est la longueur et 1 sa largeur. 


Le problème aux limites se présente donc sous la forme : 
Gite, Ty) 


0 
ox? ду? 

a, Dia y) f ED =й 
п х=0 

a, Т(х,у)+8, TE) = falx, у) 
On |, 

a, Tx. y) В, С = f,(x,y) 
n |, 

a, TG, y) f, TE =й бй 


Appliquons la séparation des variables x et y pour les solutions recherchées : 


TO ЕЕ ec GG] 

Ox Oy 
=> X" Q)YG) ХОУ" (y) 20 
x"G) Y'O) 


Х(х) = YQ) 
ЖО) = Y'O) 15 
Х(х) YO) 


La constante peut être soit positive, soit négative, soit nulle. 


2 
Dans le cas +^ , les solutions sont : 


X"(x) _ 42 — X (x) = ae* tbe * = C Cosh(Ax) +d Sinh(Ax) 
X(x) 

Y") =-2? = Ү(у) = e Cos(Ay) f Sin(Ay) 

Y(y) 


Dans le cas — Ё, les solutions sont 
X"(x) _ 
XQ) ` 
Y") 
Y(y) 


-2 2 X(x) = e Cos(Ax) + f Sin(Ax) 


=+% => Y(y) = ae? +be ” = с Cosh(Ay) +d Sinh(Ay) 
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Dans le cas 0, les solutions sont simples et pourront toujours être associées avec l'un ou l'autre des 
cas de valeur positive ou négative. 


A =0— X(x)=ax+b 
X(x) 

Y"O) _ _ 

Yo) =0—Y(y)=c y+d 


L'application des conditions aux limites permet de trouver les valeurs du paramètres À admissibles, 
soit une série en général infinie, ce qui donne la solution à rechercher sous la forme. 


T(x,y) = X(x,2,)Y(v,4.) + У Х(,А,)У(у,А,) 


An #0,n=1,00 


X(x,2) et Y(y, ^) solutions pour À = 0 


II faut préciser que le choix du signe de la constante se fait lorsque l'on peut identifier dans l'une 
des équations différentielles ordinaires un probléme de Sturm-Liouville avec des conditions aux 


limites homogènes. Dans ce cas les valeurs de À sont toutes les valeurs propres solution du 
roblèmes de Sturm-Liouville. L'identification de la direction dans laquelle les conditions aux 


limites sont homogénes détermine donc en général de maniére univoque la solution à choisir. 


Par le principe de superposition et du fait que les dimensions du rectangle ne sont pas identiques, 
nous pouvons réduire l'étude des solutions à deux problémes: 


2 2 2 2 
ав LD rene ТОУ) 0 
Ox Oy Ox Oy 
ÔT (x, ОТ (x, 
a Dis mie, TE) = доу) отоу) С] =0 
On |. On | 
ÔT (x, ÔT (x, 
a, T(x;y)+ В, E o a, T(x.y)+ p, EZ =0 
On x=k On x-k 
OT(x, OT(x, 
a T(x.y)+ f, EN =o a; Тоу) B, EZ == f(x) 
n | дп |, 
ÔT (x, ÔT (x, 
a, T(x.y)+ p, ZED Lg a, Тоу), ZEN o 
дп E On = 
y et y 
Notation 
Dans la suite de l'exposé, on note les dérivées sous la forme réduite : 
EAI rey) а ED ру) 
Ox Ô 
GT, y) ^ d T(x, y) ` д?Т (x, y) 
— © =T (x, ———— = Г, (x, — T, (x, 
ad x (х,у) ду? y (x, y) дхду „ (x, y) 
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Exemple : rectangle soumis à des conditions de Dirichlet de valeurs fixées 


d T (x, y) n Ô T(x, y) _ 


ôx? ду? | 
T(x, y), , = 0 
T(x, y)| um 
Tx, y), , =0 
Tx, y) , =1 


Les conditions aux limites sont homogénes en x, la solution à rechercher est donc pour les variables 


A ' 
—^ ets'annule enO et k: 


mE =-А > X(x) = e Соѕ(А,х)+ f Sin(A,x) 


X(0)=0 X(k)=0= À, = X(x) = f, Sin(À,x) 


Ү(у) = с, Cosh(A, y) - d, Sinh(A, y) 
T (x, y) = У (с, Cosh(A, y) d, Sinh(A, y)) Sin(A,x) 


À,#0 
Le problème de Sturm-Liouville est le suivant : 


p(x)=1s(x)=0;w(x)=1; D", (х) +4, (x) =0 
хє [0,4] 
C.L. Ф (х), =0 Ф (х) 


х= = 0 
| dx D (x) (x)=0 sinzm 


k k ‚пл k 
Í dx D, (x)? =Í dx Sin(— x)? = — 

0 0 k 2 
Les coefficients de la solution sont à rechercher en respectant les conditions qux limites en y. 


Comme la solution est nulle en y=0, nous avons cn=0, soit les seuls fonctions en sinus hyperbolique. 
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On propose une démonstration classique de l'orthogonalité des fonctions propres : 
(1) 0,@x)+4, Ф,(х)=0 хФ„(х) 

(2 D,(x)+4, D,(x)=0 хФ (x) 

C.L. 


' ' а ' ' а 
а,Ф „(х,) + 5,Ф „(х,) =0= Ф (x) zs <= Q) a,D,(x,)+0b,®,(x,) 2-0 j, (x) = = Os) 


aD,(x)+b,®,(x) =0 => D, (x) = 0,1) а,Ф „(х„) + b, (x;) =0— D, (x) = o, (а) 
€» o, (x)0, (x) - A; Ф, (x), (x) - D (х)Ф, (r) - 4, 6, (x), (х) = 0 
e - X, Jp, 099,0) +0, 0)9, (ғ) - 7 (х)Ф, (ғ) - 0 
e (д - X, o, Geib (х) + (e, 036, 09 -;, 00,09) 2 0 
intégrant sur l'interval [х„х„]= Ee -A Ја Ф (x), (x) + kb. CE, (х) –Ф (х)Ф, col se D 
be 690,6) - 6,699, Al 
(2-2) 
CL = [Ф (х)Ф, (х) Ф, (х)Ф,„(х)] =0 


soit [ax o, (x), (x) = 


(дд Ја D (х)Ф (x) =0 


Ce qui démontre également l'orthogonalité. 


D, (x) = aCos(A, х) + bSin(A, x) 


Avec les fonctions sinusoidales alors les normes se calculent à 


l'aide de la formule : 


ER 


n 


o, cf =f" gier: dauer +0" er | - 40,00, Z 


1 
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On peut retrouver cette formule de la manière suivante : 
posons À, = А, + AA dans l'équation => А2 — А2 = АДА, + A. 


[dx Aë Ai: poe, (s £ A 


substituons Ф (x) = Ф (x) (28. et en passant à la limite NÀ — 0 


en considérant les fonctions dépendant de r et А => = 4 -> - 
r r 
2 
o. (x e. (x) sa (x) Ede ) La GI. OG, (x) 0b CIA жк=ф б) ош. P A1 
Ox dÀ Ox Ox OxOA 
UA aes OG, (x) _ O (x 29,0) Ф (x) 
дА дх OxOA 
2 » 
E GE _ („уд ne 
а. CHE дА Ox OxOA 
x N 
Lt. (2, 4.) 
lorsque A. > À, 
2 » 
E GA) A reif P 
* à CO дА дх OxOA 
xY = 
(0, 24. 
2 
Е 1 дФ,(х) 1 0D (x) P Ô Z =L 64 
x дА À €x Ox 
2 2 
P 0°, (x) _ Ô (25). Ô | x o. 1 ôD, (x) 2 p.) ET (LD x 6,6) 
OxOA Ох\ дА Ох\ À дх À Ox Ox À Ox 
| x дФ„(х) дФ„(х) D, (x) (2.09 А8, о) 
“y o W = L дх Ox À, Ox E 
Í * GR 22 


n 


at = o, o et 


^s s.c) po, o «4/0, 6 )- o, 009, c); 
212 22, 


+ 1:9, er | à hew] 


n n 


= |, (0,09 = = 
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La condition en y=1, donne : 
1= Уа, Sinh(2,) Sin(À,x) 


À, #0 


En appliquant l'orthogonalité des fonctions propres 


1 k 2 k 

= dx Ф Ф =" 

T suas p to Puer; 
* dx Si 2l [^ a smt = LE ae lq -lay 
Í, dx MET [ dx Sin(x) pal Cos(x) | 2 (1— Cos(A, К)) " (1-(=1)") 


2 n 
d, = À,Sinh(À, )k a-C» 


La solution T(x,y) est donc : 


_ (2n« 1) л 
š k 
4 Sinh(A y) Sin(A x 
T(x, y) = > бу): (А, ) 
л Ze, (2п +1) Sinh(À,) 


À 


Pour le problème de Dirichlet plus général : 
Т(у), Ty) _ 


ôx? ду? É 
T(x, уу y = 
TQ, y), =0 
Tx, y), , =0 
DEE ETC 


Les fonctions propres sont : 


Ф (x) = Sin — x) 


X 
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et la solution s'écrit sous la forme : 


À, = T Ф (x)-Sin(Ax) ®,(0)=0 Ф, (/)=0 Ф' (х) = 1,С05(4, х) 


x 


pa = (aro, dauer: Aner) Aaen! 


n n 


0 


1 
l| 1 ` 1 l l 
=—|x— D", (x) | =1 Ф" (Ly == Cos? (A1) = 
> Ar) d OM 2 os (А) > 


В,= [` ах /,(х)Ф„(х) 


2 В 
T(x, y) = " Sinh(A, у) Sin(A, x 
(x, y) s зк (2,3) $їп(А„х) 


В, = IN dx f. (x) Sin(À,x) 


Si le probléme est le suivant : 


T Т 

д (y) ,2 (y) 2p 

Ox ду 
T(x, y), , =0 
T(x, y) ia 0 
ТО, у а= Лб) 
T(x, y) = 0 

. Alors la solution est la suivante, par le changement de variable y-> ly-y 

_AT fh А 2 _ L 

Anc В,= ELOD) [ do => 


2 B, ^ _ А 
T(x, у) = 1 2, 5їпї (АТ) Sinh(, (1, y) Sin(AÀ, x) 


B,- IN dx f. (x) Sin(A,x) 
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Exemple : rectangle soumis à des conditions de Dirichlet de valeurs fixées 


d T(x, y) Л Ty). 


Ox ду? v 
T(x, y), ., =0 
T(x, y), , 71 
Tv], =0 
T(x, y), =0 


Les conditions aux limites sont homogénes en y, la solution à rechercher est donc pour les variables 


2 
+Z et s'annule en 0 et 1 : 


ro =-X > Y(y) =e Соз(А„у)+ f Sin(À, у) 
Ү(у) 

Ү(0)=0 Ү()=0= 24, =пл Х(х)= f, Sin(A,x) 
X(x) = c, Cosh(A, x) * d, Sinh(A, x) 

T(x, y) = 3 (c, Cosh(A,x) * d, Sinh(À,x)) Sin(À, y) 

DEZ 

Le problème de Sturm-Liouville est le suivant : 
p(y)=1;s(y)=0;w(y)=1; Ф", (y) +4®,(y)=0 
x € [0,1] 

C.L. d,0) o = 0 GE = 0 


La P,(7)®,(7)=0 sinzm 


1 1 1 
[0 o, = |, dy Sin(n x yY => 
Les coefficients de la solution sont à rechercher en respectant les conditions aux limites en x. 
Comme la solution est nulle en x=0, nous avons cn=0, soit les seuls fonctions en sinus hyperbolique. 


La condition en x=k, donne : 
1= Y d, Sinh(A,k) Sin(À, y) 


À,#0 


En appliquant l'orthogonalité des fonctions propres 


1 1 
d = dy D 
" Sinh(A o, Gf [o 9.0) 


disi) c PS: 
[| dy Sin( ee ly Sin(y) 


1 Q= EL: _ ab. (тр 
ОЕ CDS (I) 


2 
d = ы ү, 
== L Sinh(a, pl C» 


n 
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La solution T(x,y) est donc : 


А, =(2n+1) Z 

4 1 
pO H = Sinh(À,x) Sin(À 
(x. Y) 128 Qn Siap E) DL) SCD) 


Pour les problémes de Dirichlet plus généraux : 
OT») Ty) TC»), Tx, y) _ 


ôx? ду? Or" ду? n 
Т(х, уу y = T(x,y) „= 30) 
TG, yy. , = f,Q) T(x, y). , =0 
T(x, y), , 70 T(x, y), , 70 
Т(х, у), =0 TG, y), = 


T(x, y) T, (x, y) T(x, y) 2 T; (x, y) 


Les solutions sont les suivantes 


ly H 
¿= B,=[ dy f, 0) SA») 
Ж 
2 В 
T(x,y)- "  Sinh(A,x) Sin(A. 
G7 lo кк 
et 
Dia y) 2 Y) В Sinh(A,(1, — x) Sin, y) 


1, Ze Sinh(A,L,) 
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Exemple : rectangle soumis à des conditions de Dirichlet de valeurs fixées 


OT y), T(x, Y) _ 
Ôx? ду? 

Ty), = To 

T(x,y), = T4 

T(x, у) _‚ 7 To 

Т(х,у) =T 


0 


On applique d'abord le principe de superposition, puis celui des symétries de la solution en x et en y 
et enfin le principe de linéarité. 


T(x, y) - T4, *TiGo y) € 1 y) + Tu, *T,(x, y) + To * T, (k — x, y) 


Pour ce qui est du principe de symétrie de la solution. Il est aisée de voir que pour le probléme 
OT y) PTY) _ 


ôx? ду? ; 
T(x, y), „=0 
T(x, y), = 0 
Tæ, y), „=! 
Tx, y), , =0 


la solution est T o») = T,(x,1 y) 


et que pour le probléme : 
Т(у), Ty) _ 


дх? ду? ; 
T(x, y), , 71 
T(x, y), ., 70 
T, y), , =0 
Tx, y), , =0 


la solution est ТОУ) = T,(k x, y) 
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Exemple : rectangle soumis à des conditions de Dirichlet et de Neumann de valeurs fixées 

O°T(x, y) x O^T(x, у) ü 
ôx? ду? 

TG»), =0 

Ty). =1 

Т(х,у) o =0 


0 


Les conditions aux limites sont homogènes en y, la solution à rechercher est donc pour les variables 


2 
+Z et s'annule en 0 et 1: 


YO) =-X = Y(y) =e Cos(A,y)* f Sin(À, у) 
Y(y) 


Y(0)=0 Y()20o4A,2nz Х(х) = f, Sin(A,x) 
X(x) = c, Cosh(A, x) * d, Sinh(A, x) 
T(x, y) = Y (c, Cosh(A,x) d, Sinh(A,x)) Sin(A, y) 


Le probléme de Sturm-Liouville est le suivant : 


p()-kbs(y)-20;w(y)-L Ф", (y)-40,(y)-0 
x € [0,1] 
CL. ®,(»),,=0 0,0) ,=0 


[ dy P,(»)D,(y)=0 sinzm 


1 2 _ 1 К 2 1 
Jd o, 0 [| dy Sinn» yY => 


Les coefficients de la solution sont à rechercher en respectant les conditions aux limites en x. 
Comme la dérivée de solution est nulle en x=0, nous avons dn=0, soit les seuls fonctions en Cosinus 
hyperbolique. 
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La condition en x-k, donne : 
T. y), 702 YA d, Sin(4,y) -0— d, =0 
dE 1,40 
T.(x, ») SE 1 D Ac, Sinh(A,k) Sin(2,y)=1 
т À, 20 


En appliquant l'orthogonalité des fonctions propres 


P — di o, (y) 
А, Sinh(A,k)|, (у) 
D) = Sinay) [EO => 


vss [aps 
[| у Sin( Schei у Sin(y) 


1 EN Е -(=1)" 
= [F Coso) =7-(1-Cos(2,)) = -a- C1) 


n n n 


2 
= 1-(-1)" 
2 6 тс I), c 


La solution T(x,y) est donc : 
А„=(2п+1)л 


4 1 
T(xy)-— > — 
T° Ze (2n+1) Sinh(A,k) 


Cosh(A, x) Sin(A, y) 


Pour les problémes de Dirichlet plus généraux : 
^T ^T 2 2; 
Ty) Т(у) o Ty), 0 TG») = 


ôx? ду? ôx? Oy 
TG, y) , =0 TQ y)... = f, O) 
TG»), 7,0) Lol, =0 
T(x, DE -0 TG») =0 
T(x, DM = 0 T(x. y)... -0 
Les solutions sont les suivantes respectivement : 


d A 
Ac Bc) 7,0) Sin(A,y) 
ВА 


2 В 
Ме S Cosh(A, x) Sin(À 
_ L 2 À, Sinh(A,l.) ( Z ) | ( „У) 


2 В 
Plesi n CoshlA =) Sintà 
(x, y) x 27 wan (A, , =x)) Sin(2„y) 
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Exemple : rectangle soumis à des conditions de Dirichlet et de Neumann de valeurs fixées 


ТО, у) , ё Ту) _ 
ôx? ду? 

T. y). „=! 

Tœ y), , =0 

Т(х, y). o =0 


0 


Par raison de symétrie la solution est déduite d'un exemple précédent, par changement de variable 
x en k-x: 


А, = (2п+1) л 


4 1 
T(x, y) = Cosh(A, (k — x)) Sin(À 
(x, y) 22 2 Qn +? Sinh(À, 0) osh(A, (k — x)) Sin(A, y) 


Exemple : rectangle soumis à des conditions de Dirichlet et de Neumann de valeurs fixées 

T(x, у) Р T(x, y) _ 
x” ду? 

T;(x.y)) _, =0 

Tœ y) =0 

Т(х, у), = 0 

T(x. y). =] 


0 


Dans ce cas la solution se trouve immédiatement dans la solution à constante nulle de l'équation 


séparé T(x,y)-a y+b 


Т0, y) =Y 


Si inversement la condition de Dirichlet en y=0 est nulle, alors TG. y) 21—-». si les deux conditions 
de Dirichlet sont fixé à 1, il est encore plus évident que la solution est Т(х, у) =1 
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Donc pour trouver la solution analytique du problème 


d T(x, у) К ӘТ (х, y) 


дк? ду? а; 
T. y). . = 9, 
T. (x, y). =Q 
TG»), = Ту 
Т(х,у) = T, 
La solution s'écrit 
Т,(х, y) = 4 A : Cosh(A, x) Sin(2, у) 


eni. (2n +1? Sinh(2,k) 
A Т,(х, y) +Q, T,(k x, y) -T,, 1 y) * T4 y 


Soit 
А, = (2п+1) л 
Т(х,у) = 40, > | Соѕћ(А, х) Sin(A, y) + 


л? Zei GE Sinh(À, k) 


E Le | 
2 ! (2n +1)? = k) Cosh(A, (k — x)) Sin(À, y) + T,, (1 y)*T, y 


Exemple : rectangle soumis à des conditions de Dirichlet et de Robin de valeurs fixées 


d T(x, y) " ОСОБ) 


ôx’ ду? i 
Т(у) - hTQs у) „=- 
T(x, V| =0 
T(x, Du =0 
T@(x,y)| =0 


Les conditions aux limites sont homogènes en y, la solution à rechercher est donc pour les variables 


2: 
+Z ets'annule en 0 et 1 : 


Y(y) 
Y(0)=0 Ү()=0= 24, =пл X(x)= f, Sin(A,x) 
X(x) = c, Cosh(A, x) * d, Sinh(A, x) 

T(x, y) = Y (c, Cosh(A,x) * d, Sinh(A,x)) Sin(À, y) 


DEZ 


A => Y(y) =e Cos(A, y) * f Sin(A, у) 
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Le problème de Sturm-Liouville est le suivant : 


p(y)=1;s(y)=0;w(y)=1; Ф", (у) +40,(y)=0 
xE [0,1] 
C.L. DO) o =0 D, O)| a =0 


La P,(7)®,(7)=0 sinzm 
1 


[dv 6,0 =f dy Sin(n x у) => 
, dv 9, Q)' =| dy yes 


Les coefficients de la solution sont à rechercher en respectant les conditions aux limites en x. 


Les condition en x=0 et x-k donnent : 
T(x, у), 202 Y (c, Cosh(A,k) * d, Sinh(A,k)) Sin(À, y) = 0 
An #0 


= c, Cosh(2,k)+d, Sinh(A,k) = 0 


e» = СОЮ = T(x, у) = za Cosh(2,x) - 


И Cosh(A, k) 
Е " Sinh(A,k) 1,70 


Sinh(À,k) 


О) $їп(А„у) 


=-1 


х=0 


OT(x, у) 
————— — hT(x, 
Ss (x, y) 


= Y` Che, +1, d,) Sin(A, y) 2-1 


DE 


= Yc|-h-1, M ену а 

mer] Sinh(A, k) 

En appliquant l'orthogonalité des fonctions propres 
1 


= [ dy O, (y) 
Cosh(2, k) SCH 
H Sinh(A,k) Jo. 1 
Ф„(у)=5ї(А„у) Ф, = ; 


1 | Ir 
ne p ius 


1 dg. E _ = (ту 
-4 b Cos(y)]" =—0- Cos(4,) =2-01-(-1)) 


n n n 


2 

>c, = (1-(-1)") 

А " Gen z 
Sinh(À,Kk) 


4 


€, = 0 > Cu = 


(2n+ LL i. cer | 


Sinh(À,k) 
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La solution T(x,y) est donc : 


А„=(2п+1)л 
Т(х, у) = ы py : à зна} сона, – AGE) Sn, | 

GE E Ze 

Sinh(A,k) 
Pour les problèmes de Dirichlet plus généraux : 
2 2 2 2 
e UE PH Tœ) _ qg 
Or ду Ôx ду 
Ty) -ATY „=-/,О) ТОУ) = 0 
T(x, y], =0 T. y)* hTQ у), = 7,0) 
T(x, y), , 70 T(x, у), =0 
T(x. y)... -0 T(x, DM -0 
T(x, y) T; (x, y) T(x,y) T; (x. y) 
Les solutions sont respectivement les suivantes 
пл L | 
Ac В, = |а 7,0) Sinay) 
у 

T(x, у) = z ` B, Sin(A, y)| Cosh(A, x) — COSA CARD) ACA) 

[ч Cosh(A,1.) Sinh(,l,) 

? h4AÀ,—————- 

Sinh(2,1.) 

2 B А Cosh(2,1.) x. 

T(x,y)= 2 Sin(A, y) Cosh(A, M, =x))_-= Atnbt 4 M, — 
Gy) XZ in( » ОНА, (0) r Sr, >) 


yo uade Cosh(A, |l.) 
Sinh(A,l,) 


nx 
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Exemple : rectangle soumis à des conditions de Dirichlet et de Robin de valeurs fixées 
2 2 

д AE < E 

T.(x,y)-hT(x,y)| =0 

T(x, DE =] 

T(x, Y| =0 


T(x, У), = 0 
Les conditions aux limites sont homogènes еп y, la solution à rechercher est donc pour les variables 


2 
+Z et s'annule en 0 et 1: 


Y''O) _ —А = Y(y) =e Соз(А„у)+ f Sin(À,y) 
YQ) 


Y(0)=0 Ү()=0= 24, =nr X(x)= f, Sin(A,x) 
X(x)=c, Cosh(A, x) * d, Sinh(2, x) 
T(x, y) Y (c, Cosh(A,x) + d, Sinh(A,x)) Sin(A, y) 


DEZ 


Le problème de Sturm-Liouville est le suivant : 
pQ) =1;sQ) =0;w(y)=L; Ф", (у) +40, (y)=0 
x € [0,1] 

C.L. o, (y). o -0 Ф, (y)... -0 


[dy Ф,(у)Ф„(у)=0 sinzm 


1 2 _ 1 | 2 1 
[о o, Y = |, dy Sinn x yY => 
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Les coefficients de la solution sont à rechercher en respectant les conditions aux limites en x. Les 
condition en x=0 et x=k donnent : 
TQ, y) - hT (x, y), =0=> Y Che, +1, d,) Sin, y) -0 


An +0 


=> h c„ =2„ d, >d, = p = T(x, y) = Def costa) + I Z) Sin(À, y) 


n А, +0 п 


T(x, DM =] = Seen + E ОО) Sin(,y)=1 
À,*0 n 


— En appliquant l'orthogonalité des fonctions propres 


1 1 
с, = [dy 0,0) 


[cosa sina JOO 


| 1 

D.O) = Sin») [Ф.Г =5 
1 . l ph, , 

[ dy Sin(A, y) = 24 dy Sin(y) 


1 Е РВ _ 1 BEEN 
gu oM 


n n 


D c, = : 1- (-1)) 
| 


А + N Sinh(A, К) 


4 


Con = 0; Cu = 


(2n + Dal cosi + I ST 
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La solution T(x,y) est donc : 


Д, = (2п+1) a 
4 1 
бн pr 


WEE EK [coi o $ z Sinh(A, k) 


Pour les problémes de Dirichlet plus généraux : 
2 2 2 2 
0 T(x y), Ô TG») _ Тоз») O0 TG, y) _ o 


ôx? ду ôx? ду? 
T. y) - АТО, y] „=0 Туй = 50) 
Ta»), = f, 0) TG, y) +AT, y), =0 
Т(х,у) o =0 Tx, y). -0 
T(x. y)... = 0 T(x. y), =0 
T(x, y)= Ts, (x, y) T(x, y) = Tg (x, y) 


Les solutions sont respectivement les suivantes : 


ly . 
Ai В,= | dy f,O) Sinay) 


T(x, y) = 7 Kä = Sin,» сано, + R STE 
di [Cosi E STM | 


et 


Ï o0. 
y "=l, [Cosi + 1 Sinh(A,l,) 


n 


n 


| Sin] Cos + u SIM 


T(x, y) = - SS B, | Sin(À, os D - x) I Sinh(A, (1, — 2 


Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques - p 30 
Exemple : rectangle soumis à des conditions de Dirichlet et de Robin de valeurs fixées 


D'après le principe de superposition la solution du problème 

Т(у) Ty) _ 
ôx? ду? 

T. y) -hT(x, y) = AT, 

T(x, y) ., =T, 

T(x, DE =0 

T(x, »} à = 0 


0 


est la suivante : 


T(x,y)=h Т, T,(x, y) +T, Т,(х, у) 


А, =(2n+1) Z 
T. (x, y) = 4 2. : Sin] снна, — NU sinh, 
Mr КЕ у... ee? 
Sinh(2, k) 
4 1 А Hou 
R(x, y) = Kä — cost + 7 Z 
"E (2n +1) С + 2 GT " 


7 geen EST _ И sou + ` Z 
in 
=> T, (x, y) = — D Sin(A,y) ВНЕСЕ P; = 


m Qn +D) | ^ „g Соз) 
" Sinh(A,k) 


n 


= 
С + ; Z 


T(x, y) = : 2 e ) (c, Cosh(A, x) d. Sinh(A,x)) 


pa SPACE) 
h T, Т А, 
a= BCE" h i L 
pap CAGE) ed pts би) | | Cosh(4,0 ид) 
Sinh(A, К) À, Z, 
‚ Cosh(A,k) r^ (T. -T,Cosh(A,k)) 
m Sinh(À,k) À, ___М 
AA) | | Cosh(a,k)+ 2 Sinh(a | | Cosh(A,k)+- Sinh(A,) 
Sinh(A, К) À, А, 


Sans utiliser le principe de superposition et de linéarité, on peut également tenter de résoudre 
directement le problème aux limites. 
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Les condition en x=0 et x= 


k donnent : 


T,(x. y) -hT(x, y) _ =-hT, => Y. (=h c„ +1, d,) Sin(A,y) =—hT. 


DEZ 


= T, => Y (c,Cosh(A,k)+ d, Sinh(2,k)) Sin(A, y) =T, 
À,#0 


En appliquant l И des fonctions propres 


—hc +À d 


n n n 


lag 


avec Ф (у) = Sin(A, у) 


[ dy Sin, y)  — n dy Sin(y) = ZE Cos(y)]; = =a -Cos(À,)) = q- СЄ”) 


-2hT, 


—-hc,*A,d,- 


n 


c, Cosh(A,k) * d, Sinh(2, k) = 


=> c,Cosh(A,k) * d, Sinh(2, k) = 


mu a d 


h Cons 


2п C2n4 T 


С, Cosh( À 
=> = 


E Cosh(2,„. k) + 


| 


2п+1 


d 


2п+1 À 


h 


2п+1 


А 


2n+1 


h 


Cosh(2,,..k)+ 


2n+1 


k) +d, Sinh(à 


1, сона 


ART – Ty Cosh(A,, Kk) 


n 


(1-(-1)") 


pop ito 


A a-(-1)o e, 05d, = 0 


Ee 7 


d 
h 


AT, 


2n+1 EE 


k) = 


2n+1 
2n+1 


k) +d, Sinh(A, 


2n+1 2n+1 п+1 


4T, 
k) = => 
)=7 


2n+1 


4 


Ta Cosh(A,, k) 


2n+1 


SI = 


2п+1 


SI = (T, -T,Cosh(A,„.k)) 


2n+1 


4 hT 


n 


- TCosh(2,,.k)) 


` Aami (Aon Cosh(2, 
4 h 


2n+1 


п+1 


k)+h Sinha, K) z(2n+1) (A, Cosh(A 
(T, -T,Cosh(A,„.k)) 


EIERE | 


AT - Ty Cosh(4,, ,k)) 


Cosh(A, , EI 


GEN 


2n+1 


AT, _ 4 


k)-- h Sinh(2,, EN 


2n+1 


(Т, - T,Cosh(2,„.k)) 


Connu = 


(4,, Cosh (4 


2п+1 


4 


К) + h Sinh(A,, D À 


| z(Qn«Dn| ° 


2n+1 


c 


E d o) 
dë 


"imb nasi 
| сона, 


k)+ 


п+1 Sinh(2,,. 


| 


2п+1 


h, k) + 


2n+1 


Z) 


2n+1 
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La solution s'écrit donc en renommant l'indice n des valeurs propres qui est identique à la 
précédente : 


А„=(2п+1)л 
4 Sin(A, y) А 
Т(х, у) => — ^> " (c, Cosh(À x) - d. Sinh(A. 
(x, y) = > (Ол +1) (c, osh(A, x) + d, Sinh( „x)) 
É Nu Gre 


c = 


+ A soun! 


n 


h 


d, = ^r 
+ К ST 


n 


(T, — T,Cosh(A,K)) 


Exemple : rectangle soumis à des conditions homogénes en y mixte de Robin/Dirichlet ou 
Robin/Neumann 


AT (x, y) 

T(x, у)| =0 
Ta, yy. , = f,Q) 
Ty) o =0 


aT, Gs у)+ В,Т(х,у)] _, =0 
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Les fonctions et valeurs propres sont les suivantes : 
Ф (у) =е, Cos(A,y) - f, Sin(A,y) avec ye lo, 1 | 


CL. 0,()),,=0 аФ,(у)+ Bal, =0 
C.L1 e, 20— Ф (у) = Sin(A,y) Ф (у) = A,Cos(A,y) 
C.L.2 — a,2,Cos(2,l,)+ B,Sin(A,L,) = 0 = А, solution de l'équation a,„A„Cos(2„1 ,) + B,Sin(A,L,) — 0 


n'y 


L, 


[/4v9,006,0)-0 sinzm fave oy =» p g (у) . D' (Y) H D > 9,000, Z 


1 2 1 ' 1 ' ' 
65 tx 07} x — (e, (0)o", (0) — , (1,)Ф", (1 2 


avec O,(y)-Sin(A,y) > ®,(0)=0 ®,(0)=1, 
a, @'„ (C)+ B,6,(1,) -0— a,@'„ (Ll) = -B,,(,) ®,(,)=Sin(,l,) 


Bal E ДИЕ 
= Le 25): 7. b, (L) | 


Ф„@°[1,+1, В, + B. € Sin(A,l.) joe pur ee 
E "a, A. ra, 2 d NE e ” a, 
"e Sin(2,1,) = _ 22 соб) 
[peanti fs] 

B; QA, 0 
Cos(A,1.) | L, + es 
d os( ves M ` 8? 4) 
о + afi H. 0А, J 

| ` В, ` В; 


expression alternative de 
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Les coefficients de la solution sont à rechercher en respectant les conditions aux limites en x. Les 
condition en x=0 et x=Ix permettent de trouver la solution : 


T(x, DEM =0= X(x)e Sinh(A,x) Ф (у) = Sin(A, y) 
T(x, y) = У e, Sinh(À,x) $їп(А„у) 


À,#0 


T(x, y) 


a =/,0) => D e,Sinh(A,1.) Sina, y) = fO) 


À,#0 


=> En appliquant l'orthogonalité des fonctions propres 


l 7 | 
= d Sin(À 
Sinh(A,t If EE 


À, solution de l'équation 0,А,С05(А,1,)+ B,Sin(2,l,) = 0 


c, = dy /,О)5їп(А,„у) 


m n T 
Ip Cv) FS [ + a À [ SS || 


п 


T(x,y)=2 >` c, Sinh(À,x)Sin(A, y) 
"© Sinh(A,l,)Sin(A,L," ` р Al me | 
X aA, а, 


Comportement de la solution à des valeurs limites de a2 et 82 


Lorsque 62-0, a2 =1, le probléme devient 


AT (x, y) 

тосу) 5 =0 
Т(х,у) = f,Q) 
T(x, y). =0 

Т, (х, vi. = 0 


qui admet la solution par passage à la limite 82->0, a2 ->1 et qui est également celle que l'on 
obtient lorsque l'on résout directement le probléme énoncé : 


À,Cos(A,,) = 0= À, = (n+ Dr — Sin, =1 


c, = [dy f.G)Sin(A,y) 


eee У " Sinh(A, x)Sin(A, y) 


Na Sinh(A.L.) 
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Lorsque a2=0,82=1, le probléme devient 


AT (x, y) 

T(x, y), = 

TG, yy. , = f,Q) 
T(x, y), , 70 
TQ, y), , =0 


qui admet la solution par passage à la limite a2->0,62->1 et qui est également celle que l'on 
obtient lorsque l'on résout directement le probléme énoncé (pour on utilise l'expression alternative 
de la norme des fonctions propres): 


Sin(A,1,)=0= À, = = < Cos(A„1,)? =1 


2 
о," = 2 Cos(A,1,) |1 + |1+[ CAA < expression alternative 
2 d у В 2d p 


2 2 


c, Sinh(A, x)Sin(A, y) 


T(x,y)=2 > A 
n=1,00 HE + AU Q^, | 
B; В, 
1, i 
c, = [ dy f,(y)Sin(A, y) 


2 
а, =0 3 =1= T H = n 
E EDT Zi Sei? 
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Pour le problème : 


AT(x, y) 

Т(х, у)| =0 
TQ, y), , =f,0) 
T (x, y) "s 0 


aT, (x, y) Biz, =0 


les fonctions et valeurs propres sont les suivantes : 
D (y) - e, Cos(4, y) * f, Sin(A, y) c, а =0> f, =0= @,(y)=Cos(A,y) Ф,(у) = -A,Sin(A,y) 


C.L.a,®,(y}+ В.Ф, (у), = 0 –0,4,51п(1,1,)+ B,Cos(A,L,) 2 0 — À, tq aA, Sin(A,L,) = B,Cos(A,L,) 


I, 


1 


[а Ф,0)Ф,0)=0 sinem Гаф, (у) Ada += O, url? 


D 2 p X o 


0 


1 2 1 ' 2 E ' = t 
65 tq 00) Е (Ф®„(0)Ф', (0) эм.) 


avec Ф„(у)= Соз(А„у)=> Ф„(0)=1 @,(0)=0 a,®',(,)+8,0,(4,)=0= a, (1,) =-B,0 „(1,) 


2 _1 2 В В, 2 m 2 B B, 
lo. G) = ES (d e G, (I..) | j| Ont) g [14 " || 


expression alternative de | Ф (y) 
À 


a E 
Cos(A,l,) = 23 Sin(A,l,) 


2 


212 
le, oY = Е 5 sna stu 2) 
2 B; i AL а, 
242 242 
zL Sin(A,L, Y À 23^, + E eA. [ 2.) 
2 B; a,^, B; Kë 


d. Gef = | smi eer я || 


2 


>| 
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Les coefficients de la solution sont à rechercher en respectant les conditions aux limites en x. Les 
condition en x=0 et x=Ix permettent de trouver la solution : 


T(x, DEM =0= X(x)e Sinh(A,x) Ф (у) = Со5(А, y) 
T(x, y) = dc, Sinh(A,x) Cos(A, y) 


HES 
Tan), = 7,0) > > c,Sinh(A,1.) Соз(А,у) = f, (v) 
i A, 0 
= En appliquant l'orthogonalité des fonctions propres 
1 L 
с, = dy ]„Су)Соз(А„у) 
Sinh(A,! |o, (y) La: 


À, solution de l'équation aA, Sin(A,1,) = B,Cos(A,L,) 


n'y 
1, 
c, = LU dy f, o) Cos, y) 


2 E NET 
le, Cv) = j| Ont) [ + ET [144 з J 


2°°п 


Т(х,у)=2 Y C, — ANE 7 
n-l,oo 
" Sinh(A 1 )Cos(A. LY | L 4 —2 | 14-1, ES 
inh(A.l.) 05(4,/,) ` L Al 
ou bien 
Т(х,у)=2 Y c, Sinh(A, x)Cos(A, y) 


2 2 


2 
SS HE Й | E B | 


Comportement de la solution à des valeurs limites de a2 et 62 


Lorsque 62-0, a2-1 le probléme devient 


AT (x, y) 

T(x, у)| =0 
Ta, y). , = f,Q) 
T, (x, y) i =0 


T,(x, »|., -0 
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qui admet la solution (on doit ajouter la fonction propre de valeur nulle): 


À, solution de l'équation Sin(,1,)=0= А, = E 


y 


1, l, ; 
Ф,(у) = Соз(д,у) Jo, => 0,0)=1 a [/4уУФ„(у) < Sin(A,1,)=0 
= D, et Ф, orthogonales 
1 A 1, 
е, = dy f, Qo, (y) — e, = | dy f,(y)Cos(A,y) 
Sinh(A,L.) | g | Í 


PO 


e L L 
AL['dy рало) S A- [/ d f,O) 
Ax 2 x ` Sinh(A, x) Cos(A, y) 


Pug 
VER 2 Sinh(4,) 


La solution du probléme de départ tel que calculée donne effectivement zéro lorsque 62-0, il faut 
donc rajouter la solution à valeur propre nulle. On peut toutefois remarquer que lorsque 62-50, la 
première valeur propre non nulle, tend également vers 0: 


a, Sin(A,l,) = B,Cos(A,l,) 


TUE 212 
= ы, 5 SC * O(4, ] = P, (1- "x +0(4,)) 


473 272 


2 AL, 6 АГ, 4 
= ail, —a; 6 +0(4,) = B, — В, 2 +O(Z,) 


à T D MP 
eA GR = В, +G, n -O(A,) 


° 5 , 
< À, BP = B, + О(Л„) 


E. cubo 
1? 0,1 
Las + B, J № < 


< А = 


et ce faisant, avec par exemple les conditions aux limites f,(y)=1, dans la solution: 
c, Sinh(A, x)Cos(A, у) 
T(,y)-2Y - 


"=> Sinh(A,1,)Cos(A,L, Y" [ + e = в.) 


nx 
Kë 9; 


п 
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Le premier terme de la série devient le terme dominant et se comporte de façon linéaire en x 
en tendant vers la solution lorsque a2->1 et 62->0 
cl, nzl 


2c, Sinh(À,x)Cos(À, y) 21, (¿x+ 00x) oy) H AX x 
т ` 02217100257), +1,+0(B,)) 2A, 1 
AIO 
Sinh(A,1.)Cos(A,L,. Y Ë + B. [ a) | УУ ( nor nY y у 2 хуп x 
| „А; d, 
soit 
x 
T(x, y) e 


l 


x 


Pour une fonction f,(y) quelconque, il faut isoler la première valeur propre qui tend vers 0, des 
autres valeurs propres et nous retrouvons la solution recherchée. Le premier terme est linéaire et 
les termes suivants reproduisent la décomposition de Fourier de la fonction f, (y). 

À, solution de a,2,Sin(A,1,)= B,Cos(A,l,) 


В, DB. 


t ul 
[; < aA S 
al, + p, E 
пл 


Pour n>1= À, solution de Sin(2,1,)=0 par passage а la limite B, — 0 et В, > 0> À, > E 


9 


Pour n =1= А] = 


1, 
ауес À, => c, = Al dy f,() 
Pourn>1=c,=2,[ dy f,()Cos(,») d, = |2, (Вла, +a,8,)Cosh(a,1,)+(B,8, + aa, )Sinh(2,1,) 
Premier t 7 ay tte 
>F — 
remier terme [1 ^ у f,O 


2c, Sinh(A,x)Cos(A, y) 
1 Sinh(2,1,) 


Ce qui correspond à la solution trouvée précédemment : 


Termes suivants => 


c, = | ау f,()Cos(,») 


I, 

A- [' dy 5,0) 

TG, y)- Ах 2 Y E SUINAGA) Cos(A, y) 
Db. d us Sinh(A,l,) 
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Lorsque 62-1, a2-0, le probléme devient 


AT (x, y) 

T(x, y), , = 0 
TG, yy. , = f,Q) 
T, (x, y) "S 
T(x. y),., =0 


qui admet la solution par passage à la limite a2->0,82->1 et qui est également celle que l'on 
obtient lorsque l'on résout directement le probléme énoncé (pour cela on utilise l'expression 
alternative de la norme des fonctions propres): 


À, solution de l'équation Cos(A,L,) = 0 = Sin(A,l,) =] 
А = (2n+1)r 
21, 
L, 

c, = [ dy /,Оу)Соз(А„у) 

2 c, Sinh(2,x)Cos(À 
EE 

pe $їп (А l.) 
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Exemple : rectangle soumis à des conditions homogènes en y mixte de Robin 
AT(x, y) 20;T(x, y), _, = 0; T (x, y) 


e 050) AT@,y)=0;T(x,y)| „= fO); TOs y), =0 
TG, y) - hTGs y). | =0 T,Œ y) AT (>). 0 
TG, y) Ate, =0 , Do) TG. vi. =0 
» e y 


Les conditions aux limites sont homogènes en y, la solution à rechercher est donc pour les variables 


EN et respectent les conditions aux limites de robin homogènes en 0 et ly : 

y" 

— =-X = Y(y)= e Соз(А„у)+ f Sin(A,y) X(x)= c, Cosh(A,x)+ d, Sinh(A,x) 
У 


T(x, y) = > (с, Cosh(A,x) +4, Sinh(A,x)) (e, Соѕ(А,у)+ f, 5їп(А„у)) 


À,*0 
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Le problème de Sturm-Liouville est le suivant: 

p()-ks(y)-0;w(y)-L Ф", ()«20,0)-0 xebh] => Ф, (у) = е, Cos(A,y) + f, $їп(А„у) 
CL. Ф,(у)-АФ, (у) ,-0 Ф,(у)+љФ, (у), = 

= Ф, (у) = Ate, Sin(À,y)+ f, Cos(À,y)) = Ф,(0) = 4, f, ,Ф,(0) = e, => À, f, = he, 


=> Ф (y) < Cos(A, y) * ^ Sin(À, y) 
deuxième condition = А„(—е„ Sin(A,L,) - f, Соѕ(А,1,)) + h, (е, Cos(A„1,)+ f, Sin(A,1,)) = 0 
Cos(2,1,))+h,(e, Cos(A,L,) + Ac. 


n n 


= A, (e, Sin(2„1,)+ im 


Sin(2,1,)) = 0 


= A (7 A,Sin(A,l,) +h, Cos(2,1,))+h,4,Cos(A„,)+ hh, Sin(4,I,) = 0 
=> (hh, — A )Sin,l,) 4, Cos (4,1, X(h, + h,) = 0 

À, solution de l'équation (2,° =h,h, )Sin(4,1,) = A, (h, + h,)Cos(A,l,) 

1 


y 


[/49,0)6,0)-20 sinam f gyo, O= dau + => 0, o») 50,09, Z 


0 


ELC ex. DELIS 999, o 


ES ES 
"dy O, (= 1| t| o, ue Jour le (0 O, (0) -, 0. weu )) 
[40/97/0000 0 ep, 020-0070) 
avec D, (9) = Cos(À y) * 2. Sin) PO) = Sin y) Z Cos(À,y) Ф (0) =1 
d, (0) -АФ,(0) =0 = Ф,'(0) = АФ, (0) Ф", (1,)+л,Ф (1) 209 Ф", (1) =-hD,„(,) 
1 h] 1 1 h,(,h,+1)] h 
Jo, o» = Juan dr 3 (10,07 ^e ay) - eur: pom D = 
2 

d "Op И hh, +D) h 

I, G) D c 2 — ) A 


Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques - p 43 


Les coefficients de la solution sont à rechercher en respectant les conditions aux limites en x. Les 
condition en x=0 et x=Ix donnent : 


n 


T(x, y) , =0=>c„=0 Ф„(у)= [Con + à Sin») 


Т(х,у) = $c, Sinh(À,x) oan у)+ x Sin(À, »| 


А„*=0 п 


T(x, y) 


a 0 ЛУО) => > c,Sinh(A,1.) @,(y) = f,Q) 


4,+0 


=> En appliquant l'orthogonalité des fonctions propres 


1 1, 
= "d D 
SE [а f,0)0,0) 


n 


o, (yy 


La solution T(x,y) du premier probléme est donc : 


А, solution de l'équation (2,2 — hh, )Sin(2,1,)= A, (hy + h; )Cos(4,1,) 


c, = [| 4 f, 0,0) 


2 
h,(1 h, +1 
le, o = Den T À Sin) Ü + 2 = d # x 
2 h 
TG, y) = n Sale += Z 
2 Sinh(2, 1.) db. Gel Ee 


Lorsque f(y)=1 alors nous avons 
À, solution de l'équation (2,2 — hh, )Sin(2,1,)= A, (h, + h,)Cos(A,1,) 


ou [ - d 501(4,1,) = — Cos(A,l,) 


n n 


c, 7 ['4уФ„(у) =| dy | Cos ya Sia 9) 
„=|; 49,0) =| dy +7 Sat 


1 
hE: h h 
2 Jana Т Cos(2,1,) + А 


y 
0 n 


2 
h, (lh, +1 
db Oo Henn" à sn] c ©» JA 


n 


c, = H зна, z Can) 


n 


n 


Sin(A,L,) - А Со5(4,1,)+ E 


T(x, y) = У: | : = | Saal Cost PEG Sin») 
n=1,0 А, Sinh(A,1,) |, (y)| À 


n 
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Pour mieux étudier les solutions dans des cas particuliers nous allons écrire le problème sous la 
forme : 


2 2 
Ô TUS. iB TU) 50 
Ox ду 
T(x, y). „=0 
T, у), =f (y) 


aT, (x y) - ВТО, у) _„=0 


oT, Gs y)* ВТО, у) =0 


les fonctions et valeurs propres sont les suivantes : 
Ф (у) =е, Cos(A,y)+ f, Sin(À,y) avec ye lo, L, | 


CL. оФ,„(у)- ВФ, (у), ,=0 аФ,(у)+ 8,Ф,(у) 


C.L.1 => Ф (у) = À, (7e, Sin(A, y) * f, Соз(А„у)) > Ф (0) = à, 5 ,Ф (0) =е >21, fa, = Ве, 
=> Ф (y) oC À,,&Cos(A, y) + B, Sin(À, y) 


у=1, = 0 


CL2—ajÀ, са Sin(A,1,)+ 7, Cos(2,1,))+ He Cos(4,1,)+ f, Sin(2,1,))=0 
1 1 


= а,4,(-А,0, Sin(A,l,) * B, Cos(A,1,)) + B, (A„a, Соѕ(4,1,)+ В, 51п(4,1,)) = 0 
=> À, solution de l'équation DHA — В,В,) Sin(A,1,) = À, (e; B, + Ba )Cos(A,,) 


, | , WE bl ig | 
[Га Ф„(у)Ф„(у)=0 sinzm |'ауФ,„(у) Hang tg T0) Jeu. Z 
leur, Lan url (Ф, O, (0)-Ф, (1)! (1,)) 

2 y n y A n y 2 n n n y n y 


ауес ,(y) = A,a,Cos(A, y) B, 5їп(А„у) => Ф,(0) = 4,0, ®,(0)=2,8, > $, (0) = Daum 
a, 

a @'„ (0)- Bb, (0) =0= a, @'„ (0) = B,P „ (0) 

o0, (L,) - 8,0 (L,) 0— а,Ф', (1) = —8,Ф ) 


o, Cf = СЕ L. } В É LD, (0) Bau) 


n 


D.(V)| = o. et 2 + Ё. Jean 
a n 


o. (y = Jaen) + B, HE y A et A +a, 1 
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On peut écrire une forme alternative des fonctions propres du problème : 
=> À, solution de l'équation DEE — B.B) Sin(A,L,) = A, (a5 B, + В,а,)Со5(2,1,) 


ауес Ф (у) = A,a,Cos(A,y) * B, 5їп(А y) => ®,(0)=2,a, Ф (0) = À В, 
a @'„ (0) — B,, (0) 20— a, @'„ (0) = BD, (0) 


o, = [еа + z9. vr + => (0,00, (0) Ф, (1, o, a») 


a,0',(1,)+ B, (1,) eh, Bebe dt Dé ou Ф,(1,)= SE 0) 
| Í d, 2 


o, = 5 1С E vr Lan tg, e) 


Ф (у) -3 0.04 stu, p, 1 


2/5 Kë 


Ф (у) = (баса) * f, suh x Aj +1, 3 +a, 1 


2 "n 2 

ou 

2 i| a? i E la Р 
Ф (у) ДЕСЕТО ТЕСТА 

| o l 2 
9.0 =>] “> DE 9) 
п 2 2 
P.Q) | -1 ( | Cosg -ASA (re en 
zl B B, 


2 


loy LA; 
D,» => x - E 


si p = В, =0 œ =a, =1=®,(y)= À,Cos(À,y) > | 


Ï 
si b, (y) = Соѕ(4,у) > |o, (| => 


Yy 


22 21; -2 


Lën, OH LA! d 


sia=a,=0 В = 8, =1= Ф, (у) = Sin(2,y) > |o, 0 = 
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Qui admet donc la solution suivante toujours le premier problème : 
А, solution de l'équation (a,a,2,° — B.B, }Sin(à,1,)= 4, (a, B, B, ) Cos (4,1,) 
o, (y) = À „a Cos(2„y) + P, Sin(À, y) 


0,9) = T + f, UE A A +1, 2) +a, J 


2'"n 2 


2 


ou |, of = Is Cos(4,1,) CU š afi T zl +a, J 


Së | dy f,(V)P„() = IN dy f, (y) (,a,Cos(À, y) + B, Sin(A, y) 
Tæ y) 22Y! c, Sinh(2,x) (À,eCos(A,y) + B, Sin(A, y) 


wis зин.) (uto + B, Sin(2,1,)) ү + А + I, 2) LI 
QA ` a 


2^"n 2 


ou T(xy)-2Y c, Sinh(A,x) (A,oCos(À, y) + B, Sin(A, »)) 


n-lo ,.. A aÇ [04 
КО | Cos(2,1,)- TD y aL EE =d J +a, A. 
2 


2 


si f,(y)-1 


c, = =. [A o Sin(2, у) – В, Cos(2, v) = = [2,a,Sin(A,1,)-B, Cos(2,1,)+ В] 


(2,aSin(2,1,) - В, Соѕ(4,1,)+ В, ) Sinh(2„x) (4,a,Cos(, y) +В, Sin(2, y)) 


T(x,y)=2 > 
a: So.) (nto + B, Sin(A,1,)) ` + A +1, 3) +a, Д 
Жунду: (A, Sin(2,1,)- B, Cos(A,1,) + B,) Sinh(A,x) (2„a,Cos(A,y) + B, Sin(A,y)) 


n=1,0 A A [04 aÇ 
А, КО | Cos(4,1,) - азы, |, + aL +22 т 3) +a, A. 


2 2 


Comportement de la solution à des valeurs limites de ol a2, 61 et 62 


Lorsque a1=a2=1 61=62=0, les conditions aux limites (Neumann) sont de dérivées nulles aux 
extrémités x=0 et x-lx. La solution précédente devient : 


А, solution de l'équation a,a,A, Sin(A,1,) -0— À, = T 


Ф (y) = Cos(4,y) 


" ; Iq 
c, = f; dv 0,0) - |, dy Costa, y) = aal =0 


n 


=> T(x, y) =0 
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Or la solution du problème suivant : 
АТ(х, y) 


‚ = 4,0) 


TG»). = 


avec fy(y)-1 admet la solution triviale 10У) = x Dans le développement de Іа série il faut donc 
envisager d'ajouter la solution de valeur propre nulle sous la forme du terme A x, d'autant qu'il est 
facile de vérifier l'orthogonalité de la fonction propre de valeur nulle avec les autres fonctions 
propres : 


А, solution de l'équation Sin(A,1.) - 0 = À, = 


n'y 
y 


$,0)-1 9,0)-Cos,») [P,O = 2. [dy 0,9) œ Sin(A,1,)=0 
=> D, et Ф, orthogonales 


== : sl dy f,(y)b,(y) = c, = = dy f,(y)Cos(A, y) 
Sinh(A,l,) 


0,0) 
l, 
Т(х,у)= Ах+ 2 C, UE „у) 


nx 


A[ dy = рало) = A= 


Supposons que la fonction soit de la forme AUD AO а Ж е о. nonne 


Alors la solution devient 


k 
e dy 
= nT a = Íi 20450. 3) 2 Sin(2,k) adi _k 
Г SE 1,Sinh(A,1,) A, p L 
k 2 Set? `> ots. 
fuss as). Sinh x) Cos (4; 
UB pes л n Sinha (+) COSA Y) 


> 
et si la fonction limite est constante (k=ly) alors on retrouve cn=0 et A=1 soit T=x. 


On voit que l'on peut suspecter la présence d'une fonction propre à valeur nulle dans le 
développement de la solution : 

- si la fonctions propre à valeur nulle est orthogonale aux autres fonctions propres à valeurs 
positives 

- lorsque la superposition des deux problèmes aux limites complémentaires conduits à une solution 
triviale 

- cela permet enfin dans certains cas de résoudre les solutions de problèmes aux limites 
discontinues, comme c'est le cas ici. 
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Lorsque a1-a2-0 61=62=1, les conditions aux limites sont de Dirichlet aux extrémités x=0 et x-lx. 
La solution précédente devient par simple passage à la limite 


À, solution de l'équation В,В,А, Sin(A,l,)=0— À, = E 


y 


GO, (y) = Sin(A, y) 


~x 


1 : a a 
|Ф, ОП zs 2 ( 1 Cos(A,l,) Sc? ZU T ol À, 9) + z = » 


2 2 
E . 
c, = |? dy f,O) Sin(A, y) 


2 c. Sinh(A x) Sin( À 
TG, y) = Z S © Sinh(À,O) Sin) 
l, ауД Sinh(A,l,) 


Elle coincide avec la formule de l'exemple comportant des conditions homogénes de Dirichlet sur y. 


Pour le second probléme présenté dans l'exemple : 


AT(x, y) 0 
T(x,y) = f,Q) 
T(x, DM =0 


Ty) -AT _, =0 
TG, y)* TQ, у) =0 
solution Т(х, y) 
La solution s'écrit : 
А, solution de (eo, — BB, Sin(,1,) = 2, (a B, + B. ) Cos (4,1, ) 
d 1, б 
c, = Аа LU dy /,(у)Соз(А„у)+ B [dy f, Q^) Sin(A, y) 
Tay) =2 Y с, Sinh(A,(1, — x) (4,0, Cos(A,y) + B, 5їп(А„у)) 


SS эин) (nto + B, 5п(2,1,)) б + SJ +I, 3) j ap 
[04 


2“ "n 2 


жуу c, Sinh(A, (l. — x) (А„«,Соз(А„у)+ B, Sin(À, y)) 


С зин.) Cos(2,1,)- HE D zl +4, zl F ал, 


2 


si f,0)=1= e, = 1- Dassin у) В Соу); = [a Sm) В, Cos(A,4,)+ B] 


n n 


(A, Sin(4,L,) — B, Cos(A,1,)+ B,) Sinh(A, (I, — x) (4, Cos(A, y) + B, Sin(À, y) 


CH Sn) baran +B, SH + SJ +1, 3) +a, J 
E а, n a 


2 


T(x, y) =2 Y. 


n=1,00 


(2,aSin(2,1,)- B, Cos(À,1,) + B) Sinh(A, (1, — x) (A, Cos(A, y) + B, Sin(2,»)) 


УЙ зна) | Соѕ(4,1,)- азы, |, i: 5 [ pa eU. ) +a, J 


ou T(x,y)=2 у, 


n=1,00 


2 2 
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Exemple : rectangle soumis à des conditions homogènes mixte de Robin en y et inhomogène 
mixte en x 
AT(x, y) 20 


T.(x,y)-h;T(x, y). =0 
Тху) + Тху) | = 0) 
T, Gs y) - ATQ, y). , =0 
TG y) + TGsy)| =0 


La solution cherchée est de la forme 

T(x, y) = Y` (c, Cosh(A,x) * d, Sinh(A„x)) (e, Соз(А„у)+ f, Sin(A„y)) 
са Avec un problème de 

Sturm-Liouville sur l'axe des y qui donne la même solution que précédemment : 

D, (v) = Cos(4, y) + ŽE Sin, y) 


n 


À, solution de l'équation D = hh, Jsin(À,1,) = À, (h, +h,)Cos(A,1,) 


2 
h (lL h, +1 
Io, o» Aes Si] Ü AU у 1 
. Les coefficients de la solution 


2 
À, 
sont à rechercher en respectant les conditions aux limites en x. Les condition en x=0 et x=lx donnent 


n n 


Т,(х, y)- TQ, y). = 0 = > (=h, Cn +, d,) [cosa Sin) = 0 


HEH n 


=h c, =À, d = d, = + > T(x, y) = Y« comaa) + £ Z (cosa, y) + i Sin(À, »| 


n An #0 


Ty) + Тху), =f,0) 


> Ус, BEE + H Cash $ Ол + 


À, #0 


- sans) Ф,(у) = f, (y) 


Zele +h,)Cosh(A,L,) +a, EA ОО) Ф„(у)= f,() 


EH n 


=> En appliquant l'orthogonalité des fonctions propres on obtient 


2 
1 h <œ h,(Lh,+))\) h 
lo, О = ET + Sin) 0 + É Se | 2 | 


1 
е 


C + h,)Cosh(A,1,) + Ë + Gi sni.) jo. ( D 


n 


[а f, 009,0) 


T(x, y) = У; С + T Z (cosas) + £ Lr) 


n-l,oo 
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Pour mieux étudier les solutions dans des cas particuliers nous allons écrire le problème sous la 
forme : 
AT(x, y) 0 


o T, Go y) - BT( у) _, =0 
aT y) + ВТО, у), =f,0) 
aT,(x.y)- ВТО, y). , =0 


aT, y) BST Guy) | =0 


les fonctions et valeurs propres sont les suivantes (d'aprés l'exemple précédent ): 
Ф (y) = А,а,Соѕ(А,у) È В, 5їп(А, y) 


À, solution de l'équation (a,a,4,? — p,p, )Sin(2,1,) = А, (ei B, + Ва )Соз(2,1,) 


ои 


o, = d (hant + f su e f z 


1 [es E Jean 


o, = jf | Cos(4,1,) - Aa, Sin(A,L,)f [ т zl +22 zl +a, 1 


2 2 


yn 


2 2 


| xy DO A 
si B, -fp,-0 a =a, =1> D „(y) - A,Cos(A, y) => lo, Cv) = = 


I 
si Ф„(у)= Соз(А„у) > |b, Gel = 5 


i 2 
1,0'„(,) 5А, 1 


у 


2% 227 2 


sia =0,=0 В =f, =1= Ф, (у) = Sin(2,y) > |b, = 


L'application des conditions aux limites en x=0 et x=lx donnent : 


оТ, Qs y) - B,TQ, y). _,=0= Ale, d,- B, с,)Ф,(у)=0 Vy d, = — 
HES 3 À, 
= T(x, y) = У с,(0,4,СоѕтА, х) + B, Sinh(2,x))®, (y) avec D, (y) = А,0,Соѕ(А,у)+ B, Sin(À, y) 
À,#0 


aT (xy) +AT = 0) 


> $ c, (A, (A.G Sinh(A,1.) + В, Cosh(A,1,))+ B,(A„a,Cosh(2„1.)+ B, Sinh(A,1,))) E, O = 7,0) 


An #0 


Y c (A (Ba, +a,8,)Cosh(à,1.)+ (8,8, oo, 22) іп.) Ф, (у) = f, C) 


An #0 


> En appliquant l'orthogonalité des fonctions propres on obtient 


DN Lä teg, 
" ` (2, (Ba, +a,B,)Cosh(,1,)+ (8,8, +оо )Sinh(2,L)) 


o, (yy 
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Qui admet donc la solution suivante : 
À, solution de l'équation ЖУ = B.B.) Sin(2,1,) = À, (e, B, +a, B, )Cos(A,1,) 


dh. (y) = A,a,Cos(A, y) + B, Sin(, y) 
В, В, 
aA, [ И ^ a; | t d 


1 
| 
o.) = i6 Cos(4,1,) - TIU š: zl +22 zl LI 


2 


ou | 


o, = i sacos + f, IE " 
^B 


c, = |? dy f, QA, a, Cos, y) + B, Sin(A, y) 
а, - [A (Ba, +a;B,)Cosh(4,1,)+(B,B, + 05044) Sinh(2,1,)] 
T(,y)-2y © (4,a,Cosh(2,x) + B, Sinh(A,x)) (A,aCos(4, y) + B, Sin(A,y)) 


Ke WEEN + B, Sin(À,L, v. ` SJ +1, 3) +a, À 
| a; À, Kë 


ou 


Ta, y) =2 Y) 5 
кое 20 Cos(A,l,) — 1,G Sin(A,1, Ji ` + " [ +1, a) ш d 


(4,a,Cosh(2,x)+ B, Sinh(A,x))(2,a,Cos(à, y) - B, Sin(A, y)) 


si f,(y)-1 
с, = = [4,a,Sin(A„y)- B, Cos(A,y)k = = [2,a,Sin(A,1,)- B, Cos(4,1,) + В] 
е, =[A„a,Sin(A,1,)- B, Cos (4,1, pl 


d, = "n (Ba, +a;B, JCosh(A,1,) Ë (8,8, Ж оза, ) Sinh(A,l. 1 
e, (4,a,Cosh(2,x)+ B, Sinh(2,x))(4,a,Cos(, y)+ B, 5їп(А„у)) 


T(x,y)=2 >` 
HERS A d, [e P, san) + P. +, IER 
| < a, À, ` Ü, 
ou 
Кз) (4,a,Cosh(2,x) + B, Sinh(A,x))(2,a,Cos(2,y)+ B, Sin(À, y)) 


ee? 20 Cos(4,l,) - TIU + =). 67 a) " d 
2 2 
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Comportement de la solution à des valeurs limites de ol a2 , 61, 62 


Lorsque 81=82=0, les conditions aux limites (Neumann) sont de dérivées nulles aux extrémités x=0 
et x=lx. Le probléme devient : 


ӘТ (х, y) М O^T(x, у) M 

ox? ду? 
Т,(х,у)- ТО, y)| =0 
T;(x. y) +h Ty) = 3,0) 
Tj y), , =0 


T (x, vi. = 0 


0 


D А : 2: ЛУ Т. (х,у) = Ах+В . 
Dans ces conditions apparaît une solution triviale linéaire еп х du type s (Xy) si par 


exemple la valeur fy(y)=1. Dans le développement de la série il faut donc ajouter la solution de 
valeur propre nulle sous la forme du terme A x +B : 


À, solution de l'équation Sin(A,1,)=0= А, = T 


у 
So 
Ф, (у) = Cos(A„y) lo, (| = 2 


2 
с = 


Â C + h, )Cosh(2,1.) + Ë + a Sinh(A,L.) 


nx 


| [; dy f, Coss») 


4А-ҺВ=0=В== et A+h(AL + В)= Fin Vyelou,] 


3 


— À [ + 2 L + Jr dy = E dy f,Cy) 


zc [av fo 14 [а fo) 


ІА all L, h, +h, + h;h,l, 


3 


3 n-l,oo 


T(x,y) A(x+ Ay + y С + = Z) Cos(2„y) 


Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques - p 53 


Supposons que la fonction limite est constante (fy(y)=1) alors on voit que cn=0 et donc la solution 
devient 
p h, EO h, бї Eus (Ax +1) 

h, +h, + hhl, h, +h, + hh, h, h, +h, +, ` 
Supposons maintenant que la fonction soit de la forme : 
f,(y)21-0(y-Kk) ke[O,] Ө fonction heaviside 


. Alors la solution devient 


À, = e» с, = — [ dy Cos(A, y) 
+ 1, C + h, )Cosh(2,1,) + Ë + fa. Sinh(A,l, ) 
2Sin(A,k) 
C, = 
KI. 
nz | (h, + h, )Cosh(2,1.)+| А, nr `" [ Sinh) 
_ k h, 
l, h; +h, +, 
Py) E ЕН ЕГИ 
Lh+h+hhl ` => h 
+= > HR) (снб) i R SIN Cos(À, у) 
T n=1,00 n 


n [6 + h,)Cosh(A,1,) + Ë + Gi smi 
On peut prouver que la solution du probléme initial tend bien vers la solution du probléme 
précédent par passage à la limite lorsque a1->0, a2->0 , 61-50 et 62-50. En effet on remarque que 
la premiére valeur propre non nulle, tend également vers 0 lorsque 61-50 et 62-50 : 

(210,4, – Bf.) Sin(41,) = 2, (a B, +a B.)Cos(A,1,) 
3 73 272 


e (A, aal = 2 * о = À, (ab, +a, B, 1 = ^ +0(А,)) 


213 272 


< RS — В.В, fi Jm + 24) = (aD, + e B,)(1— 2 * O(A; )) 


^r 4 
€ aa, A, 1, - pbl, + BB — +04) +0) = (a B, +a, B.) (а,в, + at B, ) 5 -* OU) 


< A, [aen * (a, B, +a, P.) z + BB; d = (a, B, +a, P.) + В.В, + O(,) 
2 (a. B, +a,B,)+ 3,8,1, 
< Z = 


Ñ I; Г 
аал + (a5 B, + GPI + 8.8, J 


A Le u A : 
< lorsque D, et В, 0; x si B, =0 et = = si B, =0 
В, al, P, ол, 
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et ce faisant, avec par exemple les conditions aux limites f,(y)-1, dans la solution: 
c, = |” dy f,O) (2,0 Cos, y) B, Sin(À, у)) 

а, - [A (Ba, 1) (8,8, oso, 21) Sinh(A,1,)| 

c, (4,a;Cosh(2,x)+ B, Sinh(A,x)) (4,0 PN y)- B, Sin(A, y)) 


°) 
T(x,y)=2 > ) 
ud WEE + B, Sin(À,1,) UE + E Ë „8 2) SS J 
ou 
T(,y)-2Y € (4,a,Cosh(2,x) + B, Sinh(A,x))(A,a,Cos(A, y)+ B, Sin(A, y)) 


GH "m Co^) Sagan) M [144 3 b Ikon 


2 


Le premier terme de la série devient le terme dominant et se comporte de façon linéaire en x 
en tendant vers la solution lorsque 01,02->1 et 61,62-50 
ici n =| 


ly . l; j 
c, = LU dy (Aa Cos(A,y)+ B, 5їп(А„у))= [7 dy (2,0, +B, A,» 003) 


l 
c, = АА, É + P > ) ~ А,1,01 


а, = |, (Во, + B, )Cosh(4,1,) + (B, B, + озо, 2) Sinh (4,1,)] 


d, = D (Lë, +a,B,)+À Baba + ЖЕЛЕ À, (ë.a, tas, +L, BB: ] 
_ 26, (A,o,Cosh(A, x) + B, Sinh(A,x)) (4,0 Cos(A, y) + B, Sin(A, y)) РЕР 


2 ә L, 

[04 

WIESEN +В, Sin(2,1,)) Ë + el +1, al +A, J rt 
Q^, [04 


2 


_ 24,10 (Аа A B + ORNA, (œu + Б,у)+ 002] _ (+в) (+В) 
¿[Ba +о, В, +1,8,8, Јо, + Bt, J 002) 2,] [Rn tap, LESE. Ja? Ва, nns BHL. 


soit 


T(x, y) = (as + Bx) 
Ba, +azB,+1,B.P4 
en posant В, = h, В, -h,,a,-1,o,-1 
(Ax +1) 


on retrouve T(x, y) =—— — —— 
EET 
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Pour une fonction f,(y) quelconque, il faut isoler la première valeur propre qui tend vers 0, des 
autres valeurs propres et nous retrouvons la solution recherchée. Le premier terme est linéaire et 
les termes suivants reproduisent la décomposition de Fourier de la fonction f,(y). 


А, solution de (оца, = B.B, ) Sin(4,1,) = 4, (a +a, B,)Cos(A,l,) 
(a, B, +a,B,) + 8,8,1, 


1? P 
[aa * (a5 B, +a, Ba) P + BB; J 


Pour n=1— À = et В, xl, si B, = 0 et = zl, si B, =0 
а 


2 2 
aA 1 


Pour n>1= À, solution de Sin(A,1,) = 0 par passage а la limite B, — 0 et B, > 0> À, > T 


avec À, 2» c, = (" dy O) (а +В A y 002) e x Аа, |" dy 7,0) 

Pour n» 12 c, - A [^ dy f,Q)CosA X). d, = A, (Buon, +a,8,)Cosh(a,1,)+ (8,8, + aa, )Sinh(2,1,) 
(a, + бух) | dy f y) 

B,a,+a,B,+1LB,B4 


24, (,a,Cosh(2,x)+ B, Sinh(A,x)) Соѕ(А, у) | dy f,(y)Cos(4,) 
d,l, 


Premier terme =< 


Termes suivants => 
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Exemple : Solution de l'équation de Laplace à deux dimensions sur des corps de dimension infini 


soit en coordonnées cartésiennes (x,y), la solution de l'équation de Laplace sur un ruban infini avec 
des conditions aux limites de Dirichlet. 
АТ(х,у)= 0 T(x,oo) fini 


T(x,y) = 0 
TQ y)... =0 
Т(у), = f.) 


La solution se développe en série de fonctions propres sous la forme : 
T(x, y) Y (с,5іп(А,х) + d,Cos(A, ol fe" ie, g^) 


n-l,oo 


CL=>d.=0 f,=0 Ae 


x 


с, = | dx f, GO Sin(A,x) 


2 
T(x, y) = = > c„Sin(2,x)e ^" 


x n=l,œ 


1-(-1)" 4 in(À. ud 
(-1) = ўе У Sin(À,,.,x)e 


si f (x)=1— e = 
f) 5 И Ж =й 2n+1 


Avec des conditions homogènes mixtes sur le ruban : 
АТ(х,у)= 0 T(x,œ) fini 


Т(х,у) , =0 
T. (x, y) 
Ty), =f.Œ) 

La solution se développe en série de fonctions propres sous la forme : 


аы G = ја f.) Sin(A, x) 


= 0 
x=l, 


А, = 


T, y) 2 Yc Sina ne” 


x n=l,œ 


. 1 4 Sin(A. Ale Set 
SI J (X =] =c = =T X, — 2n+1 
S. (x) UU (x, y) = 2 xs 
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Avec des conditions homogènes mixtes sur le ruban : 
АТ(х,у)= 0 T(x,œ) fini 


Т(у) = 0 
Т(х, y). =0 
Т(х,у) y = f.) 


La solution se développe en série de fonctions propres sous la forme : 


(2n + lr 


à, = ‚= [dx f.GO Соз(А„х) 


x 


2 
T, (x, y) = L > c,Cos(À,x)e ^" 


X n-l,oo 


1 (-1)" Cos(À, x)e ^" 
= 1 => =- > 1 = 2 5 Ш 
ЛЯ G À, 0%) = 0e 2n+1 


On vérifie facilement que : 
TG, y) - T, — x, y) 
4 — Sin(Al-Ax)e" —4«Sn(Ax-AL)e ^" 
T 1 — х, = H x n = n nx 
iG. cx y) E 2n+1 2 2n+1 


nx nx „x)) 


n= vk 
2 1)" Cos(A, x)e ^" 
)2 CI = zd, -x,y)- = ус САКЕ pe 
erp 2n+1 


Avec des conditions homogènes de Neumann sur ruban : 
АТ(х,у)= 0 T(x,oo) fini 


Tœ y). =0 
Tœ y|, =0 
TG. y), , 5A 


La solution se développe en série de fonctions propres, incluant la valeur propre nulle, sous la 
forme : 


Соѕ(4,1) 20 Sin(A,l 


nx 


јело) 
= Jam Cos(A,x) TZ 


x x 


T(x, y) =F += ` c,Cos(A,x)e ^" 


x n=1,0 


si f (x)=1>c,=0>T(x,y)=F=1 
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Lorsque le ruban est centré : 
АТ(х,у)= 0 T(x,œ) fini 


TD 8 
Туй] a =0 
T(x, y), = f. Q) 


La solution se développe en série de fonctions propres sous la forme : 


ly 
Б, 
с, = Í dx f (x) р Т(х, у) -2 Kä cn (+ ze" À, = ға 
ly 
2 


x n=l,œ 
; Lie 
1 1 n 4 Sin и кк e ^n) 


Р eg 2n+1 


si f (x) =1— е, = 


Avec des conditions homogènes mixtes sur le ruban : 
AT(x,y)=0 T(x,oo) fini 


TG, y), 220 
Т,(х, у) Led 
T(x, у), „= f. CO 


La solution se développe en série de fonctions propres sous la forme : 


Le 
да 2 e, = Í de f.c Si e) 76-2 Y евна [i Je 


2 L 1. 


2 


x n=l,œ 


$їп 20553 pw 
А 1 4 2 
аа KS 


, Tito 2n+1 


Avec des conditions homogènes mixtes sur le ruban : 
АТ(х,у)= 0 T(x,oo) fini 


(MES vi. n 9 
Т(х, Y) zc 
T(x, y], = f. GO 
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La solution se développe en série de fonctions propres sous la forme : 
1 


Quels | | 3 
de nun Гало) Cos] À, Nm 


k 
2 


TOU 2 S с Соз a DE К” 


Pepe 
, NR 
(—1)" Cos Lj, | x+ | le ^? 
ug Qno coc e tpe ji y > 
x n А, 54 z A SSC 


Avec des conditions homogènes de Neumann sur le ruban : 
АТ(х,у)= 0 T(x,c) fini 


(MES vi. n ? 
T; Gs y) -—— 
Fog. en 


La solution se développe en série de fonctions propres, incluant la valeur propre nulle, sous la 
forme : 


à à 


+= 


nz d l -5 
À,=— = | d C +211 F=—— —. 
, C, | x f (х) ofa (; 3) j 
72 


x 


T(x,y) =F + ү 25 аса» je 


x n=l,œ 


si f(x) -1 с, = 05 T(x, y) =F =1 
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Exemple : Solution de l'équation de Laplace à deux dimensions sur des corps de dimension infini 


soit en coordonnées cartésiennes (x,y), la solution de l'équation de Laplace sur un ruban infini avec 
des conditions aux limites de Dirichlet suivante : 
AT(x,y)=0 T(x,œ) fini 


= 


Т(х, V Ya] = 


[+ 4,12 = T, 


MOT" 
Deg 


On va s'inspirer du résultat de l'exemple précédent en orani le problème en deux sous- 
problèmes : 

Q={{x,y}/xe[-1,/2,H,/2|yel0,+}=Q-Q Ua, 

Q =x, y)/ x e [=1,/2,+1,/2] y el0,1, } 

Q, = Hx, yY x E [- I /2,+1, /2)у e li EJ 


duh me» 


Ta ha =T. Ta 


velo}, 


x-[-4,/2,4,/2] = 1с 
cl, 


x=[-1,/2,+1,/2] = Т, 
у=0 


I Те ‚ (x, VY) x 


Le deuxième sous-problème présente la solution suivante, oà Bn sont des coefficients déterminés 
par la valeur à l'interface des deux sous-domaines : 


T, 6,5) - У Ве” sa (e) AT 


Ld x 


Ta, Gc, y) 0] To (x, y) =0 


a +] 


s Ë l, )) 

2n+1 ES n 

Partition de l'unité 1= Í ка 2 =É >: Ше 1) loé 2623 
T n=0,00 2n+1 I À 2 


x n=l,œ n 


Le premier sous -problème présente la solution suivante, en se ramenant à une condition aux 
limites homogènes de Dirichlet : 


2 A Ï nza 
ls (x, y) =T, = 2 Sinh(2, y) Sin| А„| x+ — À, = 
ol y) 0 L, 2. Sinh(À,l,) l ( JY) 1 | | :)) 1 


х 
De méme les coefficients An sont déterminés par Іа valeur à l'interface des deux sous-domaines. Il 
vient en introduisant la partition de l'unité : 


Ta (x, y) =T, : Kä =E GE + A J — N Kä А, Sinh(A, y) SE + a 


x n=l,œ п 2. x n=l,œ Sinh(A,l >) 


E (-Cy) , smi,» |}, | L) 
То, CPE 2 É n A, ЗРНА) Sin| À,| x+ x 


x n=1,0 n 
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Pour trouver la solution, il ne reste plus qu'à assurer les conditions de continuité С, et C, à 


l'interface des deux sous-domaine, comme suit : 


ж e а zg 


n 


= a 4 Em A) a A OPAL) р CI) | Costa.) 
Sinh(À,l,) Sinh(A,l,) À, Sinh(A,1,) 
4, =, ECU gina, Je et pag CU Cosh(A,1, e ^^ 


ce qui donne la solution suivante du probléme : 
AT(x,y)=0 T(x,œ) fini 


"d m =T, T(x, WS -1.12,41,12] = Т, 


j думе Г 44 
sous la forme : 
Дыш Ш 
L 
о y =. I" "> Sinh(A, y) 16 (x+ d 
x n=l,œ n 
To, (x, ) =É Y n Deos AS "sia | x+ 3) 
x n-l,o n 
Soit encore : 
4 1 - a, 
Т, (x, y)  — T, - > -U-e "> Sinh(2 in 
Й | (2n l)e SE T Zoa Sak |; P a) 
ER d. 4 1 
PB GNETR 2, cm Сап) е" Je GH x+ 2 
Et également 
À, = g =» sifa, x+ 3) = SÉ + E > ei = Cos(i,x+nx)=(-1)" Cos(A,x) 


Ta, (x, y) = E > C 1} ( p Sinh(À, у)Соз(л,д) 


° (2n +1) 


To, (x, y) = E У ла Cosh(A,L, Je?^"Cos(A,x) 


° (2n +1) 
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Exemple : Soit en coordonnées cartésiennes (x,y), la solution de l'équation de Laplace sur un ruban 
infini avec des conditions aux limites de Dirichlet. 
AT(x,y)-0 T(x,o) fini 


Ty), = 0 
Т(у), =0 
Tx, y) „= f,0) 
La solution est : 
i 
nT I 

A e, EC) 

y 0 


тоу) Ye,Sin(A,y)e ^^ 


y n=l,œ 


Si Í, (V) =1= е, = [dy Sin(A, y) = 0) = — [i — ed 


n n 


4 Sin(2,, le 2" 
T х, = 2n+1 
(x. y) л 2, 2n +1 


Avec des conditions homogènes mixtes sur le ruban : 
AT(x,y)-0 T(x,o) fini 


TG, y). o = 0 
Ty) =0 
Tx, y) „= 30) 
La solution se développe en série de fonctions propres sous la forme : 
1 
2п+1)л f | 

a © ое ало) Sin» 

y 0 


TG, y) =Ê Ye, Sin, ye" 


y n=1,0 


À 


n n 


l 

f 1 1 
si 7,0) =1=> c, = [dy Sin(4,y)===>[Cos(2.y)]; == 

0 


4 Sin(2, y)e ^" 
T(x, y) = Nn = 
ә) T px 2n+1 
Avec des conditions homogènes mixtes sur le ruban : 
AT(x,y)-0 T(x,o) fini 


Tj y), =0 
Т(х, y), =0 
T(x, y), y 7 3,0) 
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La solution se développe en série de fonctions propres sous la forme : 


(2п+1)т 
A, = c,= [dy fO) Соз(А„у) 
21, d 
T(x, y) = E Y c, Cos(A,y)e ^" 
› n=1,0 
y 1, | Ë D" 
si 0215 e, = | dy т. Уша 
0 n n 


4 (-1)" Соѕ(А, vie ^" 
Т(х, у) = z 
у) 2. 2n +1 
Avec des conditions homogènes de Neumann sur le ruban : 
AT(x,y)-0 T(x,o) fini 
T, Gil. = 0 
T (x, D =0 


T(x,y), =J;Q) 
La solution se développe en série de fonctions propres, incluant la valeur propre nulle, sous la 
forme: 


1, 
„> L [4 £ toi 
Ai = |4 3,0) Cos.) Pris 
y 0 


y 
T(x, y) = Е+ 2. V c,Cos(A,y)e ^^ 


y n-l,oo 


1 
; А 1 | 
si 5,0) =1= c, = [dy Соз(д„у)= 4 [Sin(2,y)l; =0 
0 


=> Т(х,у)= Е =] 
Lorsque le ruban est centré : 
AT(x,y)-0 T(x,o) fini 


Т(х, y), à =Q 


Т(х,у)] 5 =0 
T(x, у) y = f,Q) 
La solution est : 
+1,42 
ASE eo [dr f,0) äi, (+012) To) 7 2 Le,Sina,(v+h/2)e ^ 
y —ly/2 y n=1,0 
+1,/2 1, 1 1 
0) =1=е,= [dy Sina (+112) = Í dy Sin(4, y) - –—-[Сох(а, у) - 7-b- c] 
-ly/2 0 n n 


) KRÉIE 
T(x, y) = 4 > Sin(4,,, g 2))e 


n=0,œ 


Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques - p 64 


Avec des conditions homogènes mixtes sur le ruban : 
AT(x,y)-0 T(x,o) fini 

TŒ. y). , =0 

Т, (x, Paa SC 

T, у) „= 30) 


La solution se développe en série de fonctions propres sous la forme : 


1,12 
А, = Qati с,= [dy fo) Sina (v +1, 12) 
2 L, -1,/2 
T(x, y) = = V c,„Sinl2, (y +1, /2)e^" 
; n=1,0 
i 1,/2 1, 1 1 
si f,(V)=1=>c„= [dy Sin(a, (y +1, /2))= | dy Sin(à,y)=- 4 [Cos(4„y)ls = — 
-1,/2 0 n n 
4 Sin, (y +1, D) ^ 
T — J 
0») 2 2n+1 
Avec des conditions homogènes mixtes sur le ruban : 
AT(x,y)=0 T(x,œ) fini 
T,(x, Ds SN 
TY) ain =0 
T(x,y), =J,Q) 
La solution se développe en série de fonctions propres sous la forme : 
1,/2 
к E [dy f,O) Cosa er 1,12) 
2 L -1,12 
T(x, y) = - Ye, Cos(, (y +1, /2)e^" 
п=1,00 
ï 1,/2 [P 1 (= 1)” 
si f,(y) =1= c, = [av Cos(A, (v PO /2))= [av Cos(1, y) 4 Sn, y) mer 
=L; 2 0 n n 


A — CI Cosa (y +1,12) 
TO. 2n+1 
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Avec des conditions homogènes de Neumann sur le ruban : 
AT(x,y)-0 T(x,o) fini 


ТУ) x 
T (x, у), = 0 


T(x, y), = 3,0) 


La solution se développe еп série de fonctions propres, incluant la valeur propre nulle, sous la 
forme : 


1,12 
gé [dy 7,0) 
a, =" a= [dy f,O) Cos. 1,72). F =>" 
L, EN l, 
T(x, y) =F + 2 У c„Cosla, (y +1, /2)e^" 


y n=1,0 
1,/2 1, 
si f,(y)=1> с, = [av Cos(à, (y+ L, /2)= [av Cos(A, y) = 2 (56а, y =0 
-1,12 0 


=> Т(х,у)=Ё=1 


n 


Exemple trivial : un ruban infini entre x=0, x-Ix, x=0, condition aux limites homogènes de 
Dirichlet 


Je ne résiste pas à l'illustration d'un exemple trivial, qui se trouve être parfois illustré dans des 
ouvrages publiés, afin de l'utiliser dans un contexte de changement de système de coordonnées. Il 
s'agit du probléme suivant : 

АТ(х,у)= 0 T(x,o) fini 

T(x, y), , =0 


T(x, y) х= = T, 


Solution | T(x, y) =T, - 


x 


Utiliser cette formule dans un contexte de changement de système de coordonnées 
(transformations conforme par exemple) permet de trouver des solutions à des problèmes en 
apparence plus complexe du point de vue géométrique dans le plan (x,y) mais qui se réduisent à ce 
problème ultra simple comme celui-ci. 
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Solution du problème de Dirichlet sur un domaine infini en coordonnées cartésiennes : une 


bande ou un demi-espace, représentation intégrale des solutions à l'aide des transformées de 


Fourier 
Cas général d'un domaine en forme de bande infini 


Sur une bande définie par x € [0,Ix] et y Є(-о°,+оо), Il est facile de construire la solution du probléme 
de er par le principe de superposition à partir des deux problémes suivants : 
— œ ee] ; ye [— oo, +оо] Q., = [0, L k [-œ, +оо] до „= {0)х [— oo, +оо] Ñ. }x [- 95,40] 


0) [e dat. Eil Salat, m0 0) Poney, o0 Tes): ЛО) 


—00 ,+œ] BS œ +оо] —00 ,+œ] —00 ,+œ] 
Supposons que bn de ces fonctions limites puisse se représenter à l'aide de sa transformée de 
Fourier, dés lors que la fonction limite est paire, par exemple pour le deuxiéme probléme : 


AO = Í da FAC y) es 0) Í dy f 00 Cos(À y) 


Et que l'on puisse de maniére équivalente représenter la solution avec son spectre de Fourier en y, 
alors en tenant compte de toutes les valeurs propres possibles venant des coordonnées y, on peut 
écrire la solution sous la forme d'un développement en intégrale de Fourier, comme suit : 


T(x, y) = {ал Т, (A)Cos(A y)Sinh(A x). Le respect de la condition aux limites en x-lx donne : 
0 


a Fh ә T, (1)Sinh( A1.) = F(A) = Т (А) = E En faisant de même pour le 
i Sinh(Al,) 


premier probléme et en introduisant également des conditions aux limites impaires, on arrive 
finalement à la solution suivante du probléme : 


ü (a) Sinh(A(L, =x)) |. FC EE бө 


de Sinh(A1,) Sinh(A1,) 
T(x, y) = [4А 
0 zo Sinh(A(1, — x, Sinh(A x) |. 
У (z M sinna) `? С ША; EI Sa 


f'G)- До) AC у) c()- AO) AC y) 


f ()= Ao f») R6)- AG) hi») 


F (a)=2 [av f (v)Cos(a y) F (a)= 2 |а у (узи) 


F (a)=2 [dy уу (s)Cosla y) À (a)= fa f (y)Sin(A y) 


La solution peut également s'écrire avec une transformée de Fourier complexe: 


_ 1 7? =.) Sinh(A(1. — x)) Sinh(A x) 
ж, LN, Aa де КАТ] = ЛТ. Je 


Dale Í dy fye” RUE Í dy (уђе 
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Sur une demi-bande entre x € [0,Ix] et y €(0,+°), 
xe[-o,] ; »e[-9] О,,=[0,/,]х[0,+%] 60,, = {0}x[0,+0]0 H }х[0,+о] ә [0. z, ]x {0} 


T(x, D ho, Е = f(y) T(x, y) ch vi = f(y) TG у) | = 0 >. P = 
Tis, =A0)= TG Teg den 
Т(х, y sal = f(x) chu) = 0 Т(х, у) р) 20 = f(x) 


Il est facile de construire la solution du probléme de Dirichlet par le principe de superposition à 
partir des trois problèmes données. Les deux premiers viennent d'être étudiés et correspondent à 
des fonctions limites impaires, car l'on impose la condition de Dirichlet homogène sur [0,1х] : 


т(х,у)= [aa k o ra À Co SP ыа ») 


ЖЕ 2 +œ | Es 2 +00 j 
avec О) = 2 [а fosna) Ё(4)= 2 |а Losna) 
0 0 7 
La solution du troisième probléme se développe en série, et a déjà été étudié : 


г [а f(x) 
pA È [ax f. (X) Cos(A, x) 9 T, (x, y) = 2—— —— 4 — Y` с„Соз(А„х)е е 


: l, 0 L x n=l,œ 
Soit la solution du problème : 


[aa ПЕ (1. -x), Ё (A) Sinh(A У sinta j 


h(A1,) Sinh(A.1,) 


кал (ez 2:6; Cos(A,x)e ^" 


L 


T(x,y)= ° 


L. 


avec F == Jos (y)Sin(2 y) F, DEE fos (y)Sin(A y) À, „=Í dx f,(x) Cos (A „x) 


L 0 
Pour le Кайл de Dirichlet sur un dent hlon où la fonction limite est donnée sur l'axe des 
abscisses : 
x e [=00,+œ] : ує[- оо, +оо] eo. M 0j Tx. у) a To) 
Il vient : 


T(x, y) = je e^" AF, (4)Cos(2 y)+ F (2)Sin(a y)} avec. f (y)= fo fex) f(v)= = 


F (a)=2 [ay у.фусо(д›) F (a)=2 ја у Gsm») 


Dans la littérature, il y a une condition pour que l'intégrale de Fourier existe sur la fonction limite 
f(x) (voir E. C. Titchmarsh, Introduction to the Theory of Fourier Integrals, second edition, Oxford 
University Press, London (1950), Théoreme 3, p. 13): 


f (y) continue, bornée sur tout interval [y,, y,], [otro est finie 
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La solution peut également s'écrire sous la forme d'une transformée de Fourier inverse (attention à 
la normalisation de la transformée de Fourier) 


т(х,у)= jaa cay "e. c(a)= Ta roe AE 70)- 1 Tay pfe 
S 2x 7, BE =Í 


=> T(x, y) = A [da f(a)e "Pe 


+00 +00 0 
Comme g. (y)= = [are ^ Ple? 2 = fas eh) + 12 GE d | Tr u: гы | 


Ол |іу+\х| i»y-|x| 


1 1 1 1 D e "ll _ d 
ur) 2: = d ү > 
8.0) Jim [nisl J2z x + y? = | Y 8. (ve 


—oo 


Théorème de convolution de Fourier  F'{F{f}xF{g}}= T xg attention à la normalisation 
V27T 


[x 


х®+(у-т) 


7970) эб) fees [ar во) р far f el far re) 
T(x,y)=F-\F\f }х Fey 


=> Т(х, y) = - [ar so) AL formule intégrale de Poisson pour le plan 


°+(у-т) 


D'ailleurs il s'agit bien de la même solution : 
Т(х,у)= [dA e^" AF, (a)Cos(â y)+ Ë (A)Sin(2 у) 
0 


T(x, y) = HT dà e ^" (F (A)Cos(A y)+ F (A)Sin(A y)}- Í ад е^“ (F (-A)Cos(A y)- F (-A)Sin(A A 


—00 


T(x, y)= Hj [da e ^F ЧЁ. ()Cos(À y)+ Ë (2)Sin(2 у)|+ +] dA e "P AE. (A)Cos(A y) F (2)Sin(2 2 


—iÀ 


ZZ y d Uu _ iAy 
T(x, y)=4 21 | di - Kë + = 


ai dà vllt „ay Ё.) (А) | us FO) HE) leen Det FQ) =iF Q) 
2 2 gu 2. 


E50) Е zlo fe?" 


F(A)*iE (A) 1" TN 
à SZ | dy f(y)e 


F,a)-— 3L f.(y)Cos(A у)= + fo fG)cos( y) => 


Р(д)== | dy Е_(у)5їп(А y) - — í dy f(y)Sin(A y) => 


FE(A)-iF(A) _ 1 f TR C. Y Pa ur 
o = 5555 dye theta) 
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Application de la transformée de Fourier pour l'établissement des formules de Poisson 
La formule intégrale de Poisson sur le demi-plan 2D, une bande et le quart de plan 


Problème aux limites de Dirichlet sur le demi-plan 


Soit un demi plan [0,+c/]xf-0°,+0], dont on donne en х=0 la fonction limite f(y), en utilisant la 
transformée de Fourier la solution peut aisément s'écrire : 


T(x, y) = [aa Je Par fa)- [dy f(y)e 7 = T(x,y)= [da F{f}x F{g}e 


+00 +00 0 

-2m|Ax| i p -27 |A x| i ] —2x|À i o 

Comme g. (y)= [але z| x| ,i27^3 = fare z| x| „1272 y + fare D 01272) 
—00 0 — 


—2хлА|х| _-i2n n -2mA|Xx —2лАу CS 
= [aA eren 4 [aa ette Ld е2", ? Лу Je 22 il, : А) _1 L | p _ 1 24 j 
d d 2x -|x| =i y -|x| +i y o 2m (dir |x-iy| Ze +y 


0 
1 M| 


5 F = e 2 = [dy g.) e 272 


—o0 


= g,(>)= 


Théorème de convolution de Fourier f *g =F '{F{f}xF{g}}=T(x,y)=F '{F{f}xF{g}}= [ar f(e)e(y—z) 


Дә fo) 26) пове [4 Ов 0) fa r0 


pl 


— 777 formule intégrale de Poisson pour le plan 
x «(v-7) 


=>T(s.y)=! far ste) 


Voici l'illustration concrète de l'application de cette formule de Poisson sur le graphe de contour 
suivant, avec une fonction limite en fonction palier d'amplitude 1 et de largeur 1 : 


T 


20r 


0.0 - 


19 == я 
00 0.5 


1.0 1.5 2.0 
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Problème aux limites de Neumann sur le demi-plan 


Soit le probléme aux limites de Neumann sur le demi-plan [0,* ee]x[-eo, * ee] dont on donne en х=0 


gro) = 70) 
le flux limite f(y) de la dérivée première selon la direction x de la fonction T : Ôx 1 . la 


solution T sera toujours définie à une constante près. La fonction T étant harmonique, sa dérivée 
première selon x est également harmonique dans le demi-plan. De ce fait on peut utiliser la formule 
de Poisson du problème de Dirichlet pour exprimer la dérivée première de T : 


Vx> 0 
TI suffit donc d'intégrer cette formule en x pour obtenir la solution : 
Vx > 0 ту) ас Ров +(у-г)) 
T —oo 


Si l'on sépare la fonction limite de flux en ses parties paire et impaire, la formule peut s'écrire : 


fG)= f `(z)+ f (z) 


Vx > 0 dE rz far Foole eoe ren ez fa пе QT md 
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Problème aux limites de Dirichlet sur une bande 
Si maintenant on se place sur une bande [0,1х]х[-°°,+°°] dans le plan, dont on donne en х=0 la 


fonction limite fi(y), et en x=lx la fonction limite f(y). En utilisant la transformée de Fourier la 
solution peut aussi s'écrire : 


y Sinh 
<| E 2 Sinh[ 


L) Sinh(2z|A|,) 


` Sim) ut) Sinh(27|A\!, - x) DIEN 
R p = L "-0- = L е 
PT A Sinhala) vorm is s a 
ИЗ Sinh(2z|A|x) RAS Nn ) 27 2 
= dà 2лАу i27m Ay 
g(x,y) Í ESO OR Eh SS S Ge e 


(rD). EE Stat, La ider: О) Гә n0» 


-i2zÀy 


lL +х—1у l,—x-—iy 
R. x x 
petes 


Comme sl x ,y)= far ш Ex) e^» = 2f À Sinh(2 À x) Cos(2x А y - 
Sinh(2x|Al1.) , Sinh(2z AL) л L, 
Sud us l, +x-iy isa l.—x-*iy als sys 1 Relv(z)-w(1-z)) ^ ReíCotan(zz) _ 1 Re Tan z x—iy 
21, 21, 2л l, 21: 2 S 


= glx, y) = i tells Ss ) = 2 Re [rad Call e(l. 2 х,у)= sr Re Coun z T ) 


_ Sinh[2z| äh). Sinh[ (1, = x)) _ jo z -x,y) n 


-i2zxAy 
7 пил 2]1, ) = јо Saile Sinh(2z|A|!,) 


Remarque l'intégrale ci-dessus est reportée dans < A.P.Prudnikov. Yu.A.Brychkov.O.1 Marichev, 
Integrals and Series - Volume 1 - Elementary Functions », formule 2.5.46.8 page 468 sous la 
forme : 


EE DEE sl" | (vs) (=) 


+ — co (** | + 
с 


La deuxième intégrale nous intéresse : 


а=2лх с=2л1, Ь=2лу g(x,y) EL Zenit: 
Sinh(c 


| sin = ТЕ 
De plus on a bien Re Tanl 7 = 2) = : = El, =x, y) = I | 
° ol сон — à zl Ул | Cos CS 
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Par le théorème de convolution de Fourier : 
Théorème de convolution de Fourier f * g = F^ {Е{} }х F{g}}= T(x, y) = К^! {F{f}x F{g}} 


f > £6) g — e(l, —x,y) f*g= [ar f(e)s(l =x, y — т) SLE gp teens +2) 


x 


fof) go а(х, у) fag [ar ceo) e = far Ao) Lee Gd 


= px. y)= Ee? | fc e| co z И )re 7 * с" | 


x x 


formule intégrale de Poisson pour la bande de largeur 1, 


Voici une illustration de l'application de cette formule de Poisson sur le graphe de contour suivant, 
avec deux fonctions limites décalées en palier d'amplitude 1 et de largeur 1 : 


0.0 05 1.0 1.5 20 25 
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Prenons par exemple des fonctions limites identiques de part et d'autre de la bande, il vient : 


D D 


f= A=) Ts) ж ы Б ы 


1 CIF 1 


MEL A [p 


x 


s = eo = ZÜ ZE 3 =k 
1+ І І r Ï 1 + Ï Ï 


x х 


M ae en) om E Pa 


£ 


l 


x 


x x 


Comme ad bel H e i Cos [E* bei 20-2) E ad? I sime =) 
L L A L L L 
su pu E=) coa "0 SM 


x 


EE E 


À 


x 


bande : T(x, y)= 1 "m P cos n iz A Se " coa L Ge zl = 


=" со 072) са be Cost 200 Je са 237) 


Voici une autre illustration de l'application de cette formule de Poisson sur le graphe de contour 
suivant, avec deux fonctions limites symétriques en palier d'amplitude 1 et de largeur 1 : 


20L 


оон 


-10t 


0.0 05 1.0 1.5 2.0 25 3.0 
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Problème aux limites de Dirichlet et Neumann sur la bande 


Si maintenant on se place sur une bande [0,lx]x[-eo, +ееј dans le plan, dont on donne еп х=0 la 
fonction limite f), et en x=lx une condition de Neumann homogène. En utilisant la transformée de 
Fourier la solution peut aussi s'écrire : 


OT(x. y) ETE Cosh(2x|A[(1, —x)) ein F T -0лАу 
=0 т = Т(х,у)= [ал F(a ” f()= [4 › 
2r ds: (ку), = fo) Т(у) Í f(a) л ш. (4) f y f(y) 
. CoshQamAx) _ 2x cs Cosh(2z.|AK, — x) "7 Cosh(2z|A|x) „ana 
R Е =L "ELE x 4 =0 = [dâ ES 
E dore Coshl2alal,) ee 5: Cosh(2r ll.) ge») | Cosh(2z AT 
31 x ^ iy l +x-iy 31. + x —iy L—x-iy 
R X X X X 
+2 Cos) an 1 dl 4, | d Al, J- "m ) d Al, ) 
Comme g(x, y) = [da QU = 
3. Cosh(2xll.) 4л 1, 
3], —x—iy L+x+iy 
Zi = =>l-7 = = 
— 1 41, | І 41, кеа 1 Refy(z)-y(1-z,)+w(z,)-w(1-7,)} 
3l, +x—iy LL —x+iy 4л x 
>= 51-2, = 
41, 41, 
Re{Cotan(x z, )+ Cotan(z z,)} _ ! Re an z L +). Zé L => |l 
AL, AL, AL, AL, 


=> gll, —x,y)= pe Tan] 7 = |+ Cotan 5 +? || = 1l me Tan| л 33. |+ Cotan| x = 
AL 41, AL, 41, 41, 41, 


Cosh(2x A (1, = x) ca "t EM -i2n Ày 
Cosh(2z|AL,) je ell, x, v)e 


En utilisant maintenant le théoréme de la transformée de Fourier de la convolution, il vient : 
Théorème de convolution de Fourier f *g = F (Ff ]xFlg]] 9 T(x, у)= F (F(f]x Flgj] 


f > f(y) g > gll, —x, y) f*g- far foel уо) far "lc Be z coran 20-27] 


x 


= T(x.»)- a le f refran z ka + cos 2 =) = Т far Heini: ce), cos 0-2) 


x 


x 


formule intégrale de Poisson pour la bande de largeur 1, 


On peut simplifier encore cette expression : 
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Pour illustrer cette formule intégrale de Poisson pour le problème aux limites sur la bande, voici un 
graphe de contour : 


N 
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Problème aux limites de Dirichlet sur le quart-plan 


Soit maintenant le problème de Dirichlet sur le quart de plan, où les deux fonctions limites sont 
données : 


Transformation de Fourier sinus F, UXA) É B | dy f(y)Sin(2 y) 
л 0 
Transformation de Fourier cosinus F {f (A) (2 Í dy f(y)Cos(2 у) 
л 0 


Convolution sinus 1243. 6 = X [ar rt ell -т|)—в (х+ DI 


Théorème de Convolution (Transformation sinus) F, sl ка) F F, {f X»), {е Y) > за) F {f УЕ, fey) 
[2 — xl e. [2 —| x 
eines fée sina р) = d io e o Eh ey 
Comme un. Des 
2 +00 IT 2 ЕТЕ 
g.(x.y)= ү” Jare Tea y)= = =a > el = Fig Larl 
0 


L(ex£l].-50) бә), =0 


y=0 


TAEA A) е =Flg.(s. y)(a) 


Solution de la forme T (х,у) у)= "deje Sin( (4 y)f, T, (aet 
пл) 2 [at назе rts rs) far riet: 
T ` (x, У, 7,)= HL £,( (= Se 2 ` a 


+Ü = гў х +(у+т) 
Formule intégrale de Poisson pour le quart de plan => T(x, y) =T (x,y, /,)+ T. (y,x, f.) 


J 


reete ey er e e) 
тз) = [ur d - L.) " AG) | 


х +(y—z) Jx? +(+r)) Leef +y (+7) m3 


|у—т))—в.„(х,у+т)) 


Voici un graphe de contour illustratif avec des fonctions limites en palier : 


Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques - p 77 


Si les deux fonctions limites sont identiques, alors la formule de Poisson devient : 


пе н о) тб) B Jae cr note ку ker 


T T T T T T T T T — T T T T T T T T =T Leg pasani 


00} 
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Problème aux limites de Dirichlet/Neumann sur le quart-plan 


Déterminons maintenant une formule intégrale de Poisson pour le probléme de Dirichlet/Neumann 
sur le quart de plan, oü la fonction limite est donnée pour la condition inhomogéne de Dirichlet sur 
l'axe des y et une condition de Neumann homogéne est imposée sur l'axe des x : 


Transformation de Fourier sinus F. YA) EI [dy f(y )Sin(2 y) 
Transformation de Fourier cosinus F ХР) = Толо fo )Cos(2 y) 


Convolution cosinus =l = = je (к — т|)+ g(x + DI 
Théorème de Convolution (Transformation cosinus) Е, PTE FAfYy)F.lgY»y)o че) Fi). eg) 


g. (x, y)= SI [dre "sin(à y)- | = e!" =F fe (y) 


x? + y? 
Comme 
g. х,у) = Nam ""Cos(A у)= (2 => ec Fg x y) 
X ty 


Т(х, у, f)... = f) 7) 


=0 
*i ХА) Е) е =F.fg.(s.y)a) 
Solution de la forme T (x, y) у)= "deje Cos (A yf (Ae ^ 


el g Gy 2) 


TG oo | [ал соңа FAO sea) T Gv D) = A= [a (rye es 


т» =) TEND 


=0 


у=0 


=> 


A x x _ 2x7 x + y? +r? 
T(x, y. f )= T | de CE +(у-—-т)” жеты) л 146 dE күте, 


Pour le problème suivant symétrique en x, y, je donne tout de suite la solution : 


= 0 


y=0 


T(x. y. f), , 7 f(x) aT yf) 


xX +y +t? 


т.з) far fe e(x-cy Xy? "T 
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J'illustre cette construction par la formule intégrale de Poisson avec ces deux graphes de contour 
avec une fonction limite en palier sur l'un des bords : 


Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques - p 80 


Application de la transformée de Fourier pour l'établissement des formules de Poisson dans un 
espace multidimensionnel 


Formule intégrale de Poisson pour le demi plan à deux dimensions puis n dimensions 


Т(х,у) „= £o) 
ту) Jar fe) 


On note que l'intégrale 1 far x] -1, ce qui se déduit de la formule de Poisson mais peut 
л х? +(у-т) 
également se démontrer d'autres façons. 


Le problème analogue à N dimensions donne une formule intégrale avec un noyau de Poisson à N- 


dimensions : 


SE )E R” 
ах" =dx,...dx„ dx" = dx,...dx, , 
T(x, Jo = f(x, Ke? Xn-1 ) = f(x, , ) 


1+” 


zl ) d: 
x,>0 T(x}=— = | ат"! f(t..) 2. ў = z 
nd D л? 


On note par ailleurs que l'intégrale suivante est égale à 1: 


PIT ; 


1 2 
R D + 


n 
2 2 
na^ R D + 


X„;-T 


n-l 


y nd 2 
nr? R (s, + 


Ce qui peut se déduire de la formule de poisson, mais peut également se démontrer d'autres 
façons. 


La solution du problème aux limites et le noyau de Poisson s'écrivent formellement : 


3; 


T(x) = Í ах"! (т, )P(x, , EE Tua ? X, ) 
pri 


n 
2 2 
л? D + 


„| 
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Pour établir cette formule, voyons si l'on peut utiliser la transformation de Fourier 
multidimensionnelle dans la résolution du problème aux limites de Dirichlet : 


KG. )E R” 
A,T(x,)= 0 WO = fe) f&a) 


La solution s'écrit à l'aide de la transformation de Fourier N-1 dimensionnelle de la fonction limite : 


Flt ) == Гах" (х, )e > L: P ws 


+27 i T„1.X 
e 


n- Xn — 


ral zi | > T; => T(x, ) = [e 7G, ) gh is d Bic 
RU! 


Comme résultat admis sans démonstration la transformée de Fourier inverse et N-1 dimensionnelle 
de la fonction exponentielle suivante est proportionnelle au Noyau introduit auparavant, ce qui par 
inversion donne la fonction exponentielle proportionnelle au Noyau 


P(x 


n 

10 

x l- 1 Ens] G+ itna- —1„—2л|х„|[т„а|„+2лїт„1. 

„зз у= п -= (2x) n [e e р |= lot tre Xni — [e e л|х„||т |+ ATERT 
d: pri pri 


-27|xn||Er-1| ` n-l ( ) -2m iT, 4X4 
=> e = dx P. X, H X, e 
RH! 


n 
2 2 
л? (s, + 


„| 


Le théorème de convolution pour la transformée de Fourier N-1 dimensionnelle a lieu en général 
(en supposant de bonnes propriétés des fonctions) : 


f(x...) 2.) (е, 1) fée 7) 


К) Ui = Í dx" f(x, , ) e?" ZE T fe} SR Í dx g(x, ) go Ka 
R"! R" 


Ft} = Í dv" f(v, , ) eg maa Fg) zz far "ele, ) eg maa 
R^ 


Alors. (f*g)(x„.)= ЕЕ, IX Fu gx, a) 


Et par identification de la formule de convolution avec l'expression de la solution du probléme aux 
limites sur le demi-plan supérieur : 


T(x,)- үте tese 79 eoe 
R”! 


gx, à ) = P ? X, ) =” T(x, ) = Í dv Eer )P(x, , = Tua ? X, ) 
к”! 
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Calcul de l'intégrale de Poisson dans le cas à 3 dimensions 


Р(х, y, z) = р 2 | P : s [a dr, e -lzi tx er, +i ( XT„+YT„ 
£ (e + y? +22) 
Plan polaire т, = 7, Cos(r,) T, = Tr, Sin(r,) dt dt, =t,dt,dtg3 x-r Cos(9) y=r Sin(9) 
[ат dr e dale emm) = Í ат, f ат, E 2x z] 
y. e ° (lz|-i (xCos(r, )+ ySin(r, DIR 0 (iz -ir Cos(r, — al (e + y? + z? у; 
1 —|2|\ e ue, +ilxz,+yry 1 Z 
le | | e ( Le Iz 3 


27 (x? + y? +2?) 


L'intégrale est facile à calculer à partir d'une première intégrale donnée раг 
« A.P.Prudnikov.Yu.A.Brychkov.O.l.Marichev, Integrals and Series - Volume 1 - Elementary 
Functions », formule 2.5.17.16 page 421 : 


f йт» 2л ч ат, а* ат, 
d 


z-i( xCos (z,)+ ySin(r, „)) Ё (x? + y? ya } x: d (2-5 (xCos(r, )+ ySin(r, Di du Je —i (xCos(r, Jr ySin(rt g ) H 


sj dt, _ 27 |z| E ат» Е (=1)7' а"! 2л 
` (z-i (xCos(z, )- у5їп(т„)))” (+y +z È » (z—i(xCos(z,)+ ySin(r, )))" (n 1)! du" (x+ y^ +u" 


DER 


Formule intégrale de Poisson pour le quart de plan 


Pour l'espace à 2 dimensions le noyau de Poisson a la forme suivante : 


Т(х,у) 


SE T : T) S T 
T(x.y)= T e cc O Ge) D 


= 7,00) Ty), =) 


La résolution du problème aux limites passe par le principe de superposition de deux problèmes 
aux limites : 
Ty] =O) (0,0) ,=0  x»0 


Thay) 22 [as тя а x SE y^ 

T,(x.y),.=f.(x) By. -0 

= T(x _ 4xy T z T я 
T( y) m [4 — ) 


Chacun de ces deux problèmes trouve aisément sa résolution dans la formule de Poisson dans le 
demi-plan en considérant que la fonction limite est tout simplement impaire : 


Tr]. =) | TG»), =0> f (=y)= (y) 


x Le E T x z Ln ous T Y 
T(x, y) x /( ne = j: et |у suu = ја /( Pon res у 


Puis on fait de méme en intervertissant x et y. 
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Le quart de plan à n Dimensions est donc l'espace situé à l'intersection des deux demi-plan 
supérieure d'un couple quelconque de direction x;>0, x;>0. Partant de la formule de Poisson on va 
également jouer sur la parité de la fonction limite : 


SS n _ n-l _ 
X, xxx, le R Kai ne ЈЕВ. TE). fe, ) 
nl _ dr. dr, » t) ) DES SÉ ) Xi 
dx =Z" x,>0 T(x < | а, F (Tin : 
X; RU 2 2 
X; + Xin Ti, 


Si la fonction limite est impaire en x; alors la formule est également solution du problème aux 


Т\х = X. T X, -0 4 А š 
limites : ( re Wi ) ( À, , et en introduisant de nouvelles notations pour 
distinguer la variable d'indice j, on peut écrire : 


n-2 n-2 dx,...dx 

X; = (x X 15 X; X353; S X,JeR,, dX, >s == m 
i,j,n-2 ` 19009 ie GG EA C GA Ad j—15 "^ j+1 5 Hei i,j ij 

ах,ах, 


Tini = (r, ys Cas fas Tos x,) Fis, eee.) fonction impaire ХР; Le. ))= «fts 1) 
n 
n- x; X. 
T(x,)- ES de, "^ CH dz ,f (s I - i 
n-2 2 2 Es 2 2 о |5 
л? 8. и + X, j,n-2 = Ti jal +(х, = Š P + X, j,n-2 e GE +(х, +r) IE 
En construisant la méme solution par interversion des indices i,j, il vient la solution du probléme 
aux limites : 


n-2 n-2 
X, jn2 = (кызбы ya Xo pase lei R; ; dx;; = 


T(x,). , = f(x...) T(x,). -o =J (х) 


di í йт [ату (т, К x, š 
Wi d 
EE Ee 


D D 
2 2\5 2 2 E 
2 2 
Xi VC 2 Tin] +(x, -т,) | " * Xi j,n-2 =T; j.n-2) +(х, geb 
+00 
X; X; 
n-2 j j 
23 | йт, [arf Сте Я | 7 
R. "2 0 _ 2 2 _ 212 
dé E + |x, j,n-2 mal * (x, т) | s + X, j,n-2 RE * (x, T т) 


DIS 


Voyons ce que cela donne à 3 dimensions, sur le quart-plan x>0, y>0 : 


T(x,y,z) = fz) Т(х,у„г] „= /,(х,2) 
far, far, AGES ëd I Ш ypt 
bu? 9 (e +(z—t, Y +(y- T Jp)? (e +( +(z- т.) Ales Pl 
Т(х,у,а)=у— 
+ Тат, (ату, (т.т. а = > = 
œ 0 (y *(z-v,y +(x ny (v + (z t.) exe, Y)? | 
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Formule intégrale de Poisson pour le huitième de plan en 3D 


Les conditions aux limites sont maintenant : 


T(x, y. z) ved yox) Thx, y. z) „= f,(x,z) T(x, y, z) 


SEH = f(x,y,z) 


En jouant sur la parité des fonctions limites à l'égard de la variable z, on obtient facilement la 
formule suivante : 


x x 
e к (x2+(z-_r.) 6-50) еу re | 
(een Gs)? | EMI 
У y 
EE + far. fanis) (е0) Ааа) Я 
ba) а. 

KAP ОООО е 

| (pud exiens] 


Formule intégrale de Poisson pour une bande plan de largeur Ix 


T(x,y) 2 A6) = TG»), = £O) 
Pede = ја | f (ne cons cg) fee ras => =) 


x 


ЕСО ош 
"ec cost L ud C QU cost L eed 
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Pour le problème multidimensionnel, la solution peut s'écrire sous la forme d'une transformée de 
Fourier (ou son inverse) : 


Ts.) sl TG) 


x„ =l; E h (x, ) 


al Sinh( T К =X )) = Sinh[ т, |X ) ; 
T = d n-l п=1 (Ах п 1|n +2літ, 1-Х, 
&J- [ cie, E, | ‚йе, EE 
Ais, [as f(x emm Zi ^ dx"! f (x, e "nte 
R” 
Remarque Lim Sinh a d Lim ea) =Q = Lim Sinhl2z E ER Lim eK) — 0 
nal Sinh(2z|e, p [eys] >+ [a +œ З т, |. e 


Sinh( 


Së n-1 X, ils) +2літ, 1.X„.4 n-l 
g(x,,.x,)- [ dt e sl x ct "p ат 
gu Tua "e = Geh 


Cas à 3 dimensions, l'intégrale à calculer est la suivante : 


„i 24 т X, ) ein LG ЖИГ dE 
1) 


Tri 


Plan polaire т, = t,Cos(r,) f£; — Sin(r,) dr, dr,-T,dvt,dt, x=r Cos(9) y-rSin(9) 


Sinh( z|) Miles) = Sinh( dÄ. 
Í = dë зноя) ” I=I-I= Jar, [arr T. оры) 2 s(2z r,(xCos(r,)+ y Stats, II 
2л +00 
=> Í ap |49 Гаа и PE GE Sin(r, ))) 
inh(ul,) 
— Calculons [az z, SR t, ) = 2 far, Sinh(a т T.) Sinb т,) 
d Sinh(c c, ) db} ' Sinh(cr,) 


с 
bz ал 
1+ Cosh C 
Nn Sinh(a v, ) л? ui | с | a c | 1 (a+c+ib (К a+c-ib 
= [атт MA Cos(b tr, )= ;|v +w 
inhic z, 2c 2c 


On obtient une intégrale simple pas forcément trés facile à calculer sous la forme suivante : 


Plan polaire т, = т,Соѕ(т,) T, = т,Ѕіп(т,) dt dt, =т,ӣт,йту х= Cos(9) y-r Sin(9) 


14 соңа — Ü d 


D 


x 


Sinh(2z "E iani ) = uds [ar xCos(r, )+ y чы + са) |z d 
9 
l 


dr = 
|е E (2л z, 1 .) dh 


x 


` g sicot) s), (i etlech? = 


21 


x 


È 


A ce stade je n'obtiens de formule explicite pour la fonction g qui permettait une formule de type 
Poisson en faisant jouer le théorème de convolution. 
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Solution de l'équation de Laplace à deux dimensions en coordonnées cylindrique (r 9) ou (r,z) 


avec conditions aux limites Dirichlet, Neumann ou Robin par la méthode de séparation des 
variables 


Le systéme cylindrique est défini par le changement de variable : 


х = т Соѕ(9) 

y = r Sin(9) 

z 

Typiquement nous recherchons la solution du problème aux limites intérieur en coordonnées 
sphérique sur un espace limité dans l'espace des valeurs ” © [0.1] ; 6eelo.2z] : ze 


L'équation de Laplace en coordonnées cylindrique 3D (r,O,z) s'écrit sous la forme 
' 2 2: 
1 ð „оТ (7,0,2) | 1 G^T(r,0,z) Я д TUS) 8 


2 2 0 
r òr Or r 00 Oz 
Le Jacobien en coordonnées cylindriques et la métrique sont : 
à x, 
26 Ge Cos(0) -rSin(0) 0 T MA 
J=det Z L 0 |= del Sin(0) rCos(9) 0|= 0 5 -= HCos?*9+Sin*9)=r 
or 00 ^ Š Í òr 00 O0 dr 
0 0 1 


ds? = ар? + r!d9? + dz? 
Ѕіпаи/агіёё du système de coordonnées cylindriques 2D-3D et régularité du gradient 


Autrement dit, le Jacobien s'annule en r=0. Dans le cas où le domaine d'étude de la solution de 
l'équation de Laplace comporte l'origine r=0, alors le système de coordonnées cylindrique 
comporte une singularité. Cela se traduit par l'annulation du Jacobien en r=0 et la singularité 
"apparente" du Gradient de toute fonction solution de l'équation de Laplace. La solution doit donc 


neutraliser la singularité apparente du gradient en assurant une valeur bornée à ce dernier soit : 
(7r,0,z)eQ/r=0 


Grad(T(r,0,2))= (2066.2 y, " 1 (222), (£682). 
r 


Or 00 ez 

2 2 2 

Grad(T(r,0,2)f =| Z Git | eo E 
Or r 00 ez 


Singularité du Jacobien d =r =0 condition Grad(T(r,0,z)) reste bornée Vr,0,z 


1 E 


> Cste 


r? 00 


r=0 
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Cette contrainte n'est en général pas évoquée car elle est de facto remplie par la contrainte de 
finitude de la solution: Exemple : les solutions admissibles par séparation des variables (voir plus 
loin) sont : 

T(r,0,z) = R(r)O(0)O (z) 

R(r)= ar^ +br ^ > finitudeenr =0= R(r) =r" = R'(r) =A r^? 

@(0)= e Cos(10)+ f Sin(10) => O'(0) = —eA Sin(20)+ fA Cos(20) 

Le respect de la finitude de la solution pour les valeurs r=0 suffit pour assurer la condition de 
finitude du gradient soit : 


2 
LOAD od œ r?*(q Sin A0) + B Cos(A0)) œ r”? 
r 00 
2 
l je) > Cste si22>l 
r? 00 E 


Dans le cas d'un disque complet les valeurs de À sont égales à n 21, pour le cas d'un secteur de 

disque on А, = THE иеаз A >n>lr Lorsque l'angle d'ouverture du secteur de disque est plat 
a 

(а=л), les valeurs propres également égale à n. Il reste le cas où l'angle d'ouverture du disque est 

plus grand que n. 


Dans ce cas, on peut faire appel au raisonnement impliquant la construction de la solution sur le 
problème complémentaire du secteur au disque entier. Appelons problème n°1, le problème sur le 
secteur d'angle obtus. Il y a existence et unicité de la solution du problème de Dirichlet intérieur 
construite par les fonctions de l'angle obtus. Cette solution s'étend mathématiquement à tout le 
disque, donc sur les frontières du secteur d'angle aigu. Cela induit donc un problème aux limites de 
Dirichlet n°2 sur le secteur d'angle aigu. Ce dernier se construit également à l'aide des fonctions 
propres qui possède la régularité du gradient requise sur le système de fonctions et valeurs propres 


В=2л-а å, = azn,B<n=A,>n>1. Et par existence et unicité de la solution induise 


le respect de la régularité du gradient de la solution dans son domaine complémentaire qui doit 
être identique à la solution du problème n°1. Ce type de raisonnement s'étend à d'autres types de 
configuration avec des domaines complémentaires. 


Dans les configurations de coordonnées rz, la singularité du gradient n'est pas présente car les 
coordonnées ne dépendent plus de l'angle 0. 
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Équation de Laplace en coordonnées cylindrique polaire (r,9) 


Lorsque l'on réduit le probléme à deux dimensions , les deux systémes étudiés sont le systéme 
polaire (r,O) et le système cylindrique axial (rz). En coordonnées polaire (r,0) l'équation devient : 
10 ÔT(r,0) l10O7(r0) OT(r,0) 1 ôT(r,Q) 1 ô*T(r,0) 
r +> — = SE L ;—-0 
r òr Or r 00 Or г Or r 00 
Et en coordonnées cylindriques (rz) également dénommées axiale ou axiale-symétriques: 
10 OT'(rz) O’T(r,z) OT(rO80) 1ôT(r,0) OT(r,z) 
r + 2 = 2 + + Э = 
r òr Or Oz Or r дг Oz 
La séparation des variables en coordonnées polaires (r,O) donne : 
OT(r,0) l0T(r,O) 1 O°T(r,0) 
2 + + 2 2 
Or г дг r 00 


= R"(@(0)+ Á moe) + — R(r)@"(9) =0 
r r 


0 


-0 T(r,0)- R(r)G(0) 


= rR'(r)G(0)*r R'(r)@(0)+ R(r)@''(0)= 0 
rR'(r)rR'(r) 2 @''(0) 


R(r) Ө(0) 
rR'(r)r R'(r) _ @"(0) 242 
R(r) |. 0e) … 


PR (r)+r R'(r) e £A R() = r (r R'(r))-+ ÆR(r) = 0 

Ө"(0) + 20(0) = 0 

La constante peut être soit positive, soit négative, soit nulle. Les valeurs doivent être choisit еп 
cohérence avec la configuration géométrique. 


Dans le cas A , les solutions sont : 
r (r R'(r))-A R(r) 202 R(r) = ar^ «br ^ 
Q"(0) + А?Ө(Ө) = 0 = @(0) = e Cos(A0) + f Sin(A0) 
Dans le cas — Ё, les solutions sont 
r (r R'())V'HVR(r)=0= R(r) = ar? + br ? 
" 2 ма _ 20 -10 
Ө"(0)-10(0)=0= Ө(0)= ae" tbe . Mais dans ce cas la solution en variable Ө n'est pas 
périodique, car la conformation géométrique impose que la solution soit périodique 


Ө(0)=0(0+27) lorsque le problème est posé sur tout l'intervalle des valeurs де Ө, soit sur le 
disque complet. Pour un probléme sur un secteur restreint du disque g Io. 0,]c [0,27], оп 


pourrait penser а priori, qu'il n'en est pas de même. Les solutions de l'équation différentielle 
ordinaire en г : 


R(r) = ar^ «br * 
R(r) = aCos(2 Log(r)) + bSin(A Log(r)) 


Présentent un inconvénient majeure, celui de comporter une oscillation rapide lorsque r->0 rendant 
instable toute solution théorique, totalement impropre à représenter un phénoméne physique. On 


2 
écarte donc les valeurs ^ À . 
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Dans le cas de la valeur propre nulle, les solutions sont simples et pourront toujours être associées 
avec l'un ou l'autre des cas de valeur positive ou négative. 


г (г R'(r) 20 R(r)=a+bln(r) 
Q'"(0)20— @(0)=c+d0=c car O(0+27)=O(0) 


L'application des conditions aux limites permet de trouver les valeurs du paramètres À admissibles, 
soit une suite en général infini, ce qui donne la solution à rechercher sous la forme. 


T(r,0)=a+bln(r)+ D R(r,¿)9(0,2A,) 


HEURE? 


Lorsque le problème aux limites conduit aux solutions fonctions propres O(0, 4. y, ces fonctions 


sont orthogonales sur l'intervalle o = [0,,0, ]c [0,27 |; Soit : 


N 10 (8,2, (0,2, ) 0 si n # m 

O(0,2,)=Sin(2,0) ои @(0,2,)=Co0s(2,0) 

1-Cos(20) 1 ГӨ Sin(20)|"^ 
À | | à 


x 2 9, s: 9, 
OH =Í do Sin (4,0) = |, 10 ; 2 A 


JCos(4, 6) = N d0 Cosi. 0) = f. 40 


1+Cos(20) _ 1 É , seo 
2 A2 4 


Am 0, 


Ici le choix du signe de la constante se fait sans trop d'ambiquïté conduisant à un probléme de 


Sturm-Liouville en coordonnée Ө 


Le gradient dans ce systéme de coordonnées prenant la forme : 


Grad (rir. oy) =( LQ y, (25650, 
Yr r 


00 
Les conditions mixtes de Robin sur une iso-courbe r=Cste (cercle) ou ü=Cste (cône) prennent alors 
la forme : 
Iso — courbe r = Cste Iso — courbe 9 = Cste 
OT(r,9 l10T(r,9 
C.L. a => = UM -BT(r,9) = (9) C.L. de S ) -BT(r,9) = foo (r) 
a dQ 


Si les conditions sont homogènes, et dans le cadre de la méthode de séparation des variables, il 


vient pour une courbe O-Cste (cône) 
Séparation des variables => T(r,9)= R(r)8G(8)) 


1 дӨ(9) 
a = t BO) 


=0 


9-9, 


C.L. 


La variabilité radiale de la condition aux limites introduit donc une difficulté de résolution 
supplémentaire puisque le systéme de valeur propre et fonction propre angulaire dépendrait alors 
de la position en r. Dans ce cas de figure on simplifie le probléme en supposant, de maniére tout à 
fait académique et certainement irréaliste, que la condition aux limites mixte posséde des 
coefficients indépendant de l'espace. Ce peut étre réaliste pour une configuration de section 
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annulaire creux où l'épaisseur radiale est faible et donc le facteur r variant peu. Mais certainement 
faux pour un cône complet. Car dans la limite des rayons faibles, la condition aux limites se 
comporte comme une condition de Neumann, tandis que dans les rayons grand, c'est plutôt une 
condition de Dirichlet. En conclusion les solutions sont exhibées volontairement avec des conditions 
mixtes fixes, comme suit : 

Iso — courbe 9 = Cste 


a SO) + BO(9) = 0 
og 9-9, 
ce qui permet d'ailleurs de condenser en une seule formule, les deux cas limites Dirichlet et 


Neumann. En revanche les problémes avec condition de Neumann ou de Dirichlet sont eux 
susceptible d'étre résolues entiérement par la méthode de séparation des variables. 


C.L. 


Équation de Laplace en coordonnées cylindrique axiale (r,z) 


La séparation des variables en coordonnées cylindrique axiale (rz) donne : 
OT(r,z) lOT(r,oz) O'T(r,z) 
+ + 
or? r Cr oz? 


20 T(r,z)= R()Z(z) 


Во) (о) ВР) (с) RAZ") = 0 
r 


= rR'(r)Z(z)- R'(r)Z(z)+ rR(r)Z''(z)= 0 
rR''(r)+ R'(r) ne Z"(z) 


rR(r) Z(z) 
rR''(r)+ R'(r) SR ZH) 242 
rR(r) Z 


Ce qui conduit au système de deux équations aux dérivées partielles. 


rR''(r)+ R'(r) = EArR(r) 2 (rR'(r))-X XrR(r) = 0 
Z'(2)+2?Z(z2)=0 


La constante peut être soit positive, soit négative, soit nulle. Les valeurs doivent être choisit en 
cohérence avec la configuration géométrique. 


2 
Dans le cas +4 , les solutions sont : 


(rR'(r))- X rR(r) 202 R(r)=al,(Ar)+bK, (Ar) 
Z'"(z) - AZ(z) = 0 = Z(z) = e Cos(Az) + f Sin(Az) 


où DANK) 


modifiée. 


sont les fonctions de Bessel modifiées de première et deuxième espèce 
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Dans le cas — Ё, les solutions sont : 
(rR'(r))- X rR(r) 202 R(r) = aJ,(Ar)+bY, (Ar) 
7''(2)— AZ(z) = 0 Z(z) = e Cosh(Az) + f Sinh(Az) 


Oü J (Ar) Y (Ar) sont les fonctions de Bessel de premiére et deuxiéme espéce respectivement. 


Dans le cas de la valeur propre nulle , les solutions sont simples et pourront toujours étre associées 
avec l'un ou l'autre des cas de valeur positive ou négative. 


(rR'(r) 202 R(r)=a+bln(r) 
Z'(z)202 Z(z)=c+d z 


L'application des conditions aux limites permet de trouver les valeurs du paramètres À admissibles, 
soit une suite en général infinie, ce qui donne la solution à rechercher sous la forme. 


T(r,z) = (a t bln(r)Xc +d z)+ > R(r, ¿,)Z(z,1,) 


À,#0,n=1,00 


Lorsque le problème aux limites est homogène soit dans la direction r, soit dans la direction z les 
valeurs propres À sont solutions du problème de Sturm- Пову. C'est l'homogénéité des 
conditions aux limites qui dirige le choix des fonctions et valeurs propres de l'équation séparable. 
Pour la direction r le probléme de Sturm-Liouville conduit aux fonctions de Bessel de premiére et 
deuxiéme espéce. 


ptr)=r;wr)=r;s(r) = dira zr PO zr Ro) =0 refr,r,] 


C.L: 

a R(n.) * bR'(n) =0 

a;R(r,)+b,R'(r,)=0 

Admet comme solution sur [r1,r2] les fonctions propres P „(x) qui sont toutes orthogonales entre- 
elles sur l'intervalle [r1,r2] pour le produit scalaire de mesure r 


[arr(M)O,MO„()=0 sinzm 
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Pour se convaincre de l'orthogonalité on utilise une démonstration souvent utilisée dans la 
formulation des problémes de Sturm-Liouville. Dans notre cas les fonctions propres sont solutions 
du probléme aux limites, et supposons deux fonctions propres de valeurs propres différentes : 


(o) + rA D,(r)=0—=®,(r)= = (0,0) 
r 


n 


(но, (5) +522 Ф, (r)=0= 9,0) = -— (ro) 
r m 


C.L. 


' ' a ' ' a 
а,Ф (қ) * b, (n) = 0 = o, (1) == c ds) a5G, (r,) * 5,05, (r) = 0 = D, (7) =; RES 


1 2 


i , a a 
aP „(.) * jb, (n) 2 0— Ф,„ (л) = = 9.) aP, (r;) c b, (5) =0>> D „ (12) = = Фо) 
1 2 


Par une double intégration par partie 

= [arr 8„(r)O „(r) = -5f dr (тФ(г))Ф„() 
1 ' > m ' ' 

ООСО) 


n 


-- => (67 6,08 -[o, c) (7, o + [ 49, GE, 2) 


Š re. Ф, c] - Io. e Act -g f aree) 
or Ho, tc cl -[ro, c) 9, cT А 


a, а us uy SCH = 
É н S n n 9.69.6) | SA Pat) tn.) 0) 0 


2 1 2 1 


2 D 
= ['arr®,0®,0- af dr à, (ry, (> ("dr D,(M)O„(r)=0 sih, 2, 


Ce qui démontre l'orthogonalité. 
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Une autre démonstration est équivalente et consiste en ce qui suit : 


M (rO0,())+H22 Ф, (0) =0 ХФ, (Р) 
(2) (ro, (r) )+r22, D,(r)=0 x®,(r) 
C.L. 


' ' a ' ' a 
aD, (rn) +b,P (7) = 0 = o, (1) =< Don) a5, (r,) * bo, (r,) A 0 = D, (r,) = = Din) 


1 2 


I ' a ' ' a 
aP, (n) * b, (n) -0— $,6) 7 -,- D, (ri) a5, (r) t b, () - 0— Deler 0,0) 
1 2 


e (н, (cht. Ae) +12 0,00, 0) CMO LAO - 72, Ф, 0)9, (7) 20 
er(Z-2 ),00,0«69,0)9,0)-(9,(0)9,0) 20 
orle -2 y, (90.0) «(0)», 0)-9,0)9,0))) - 0 
intégrant sur l'interval |г„т,]= (a2 -A, fv r D, (r)b, (r) + BO (r)D,(r)-®,(r)®, (J: =0 
Le 09,0) - 6,09, (el 
(22 =X, ) 
et comme précédemment BO (r)®,(r)-®D,(r)®, (JE =0 


soit [ar r $,(r)b,(r)- 


=> (2-2 f атт Ф„()Ф„()=0 


Ce qui démontre également l'orthogonalité. 


Le problème du calcul de la norme des fonctions propres n'est toutefois à ce stade pas résolu. La 
norme des fonctions est calculée pour chaque problème de Sturm-Liouville sur le produit scalaire 


Jo, c = arro, (0) 


de mesure r, . Avec les fonctions de Bessel de premiére et deuxiéme 


PUIS ar оч РИ) , alors les normes se calculent à l'aide de la formule 


kË D", (r)? +®, eJ 


2 


espéce, 


Ip cy = Ñ dr r D „(r)? = 
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On peut retrouver cette formule de la manière suivante : 
posons À, = А, + АА dans l'équation = À - X, = AA(A, + А„) 
Le. vc Ae) Ach, (ell 

(a -An ) 


substituons Ф (r) = D (r)* (22m et en passant à la limite АА — 0 


[ar ro, (п) (т) = 


en considérant les fonctions dépendant de r et А = = g — = 
r r 
oo r = r 
D, (r) D D, (r) = Da 
2 
ЕБС) ÔD „(r) | SS 06, (r) Айу AD 0) a°, (r) m 
е 04 J ô ôr rô 
2 
x АЛ ôO „ (z) 0, (r) -o (r) PauU) 
дА Or OrOA 
2 n 
(20 O A vi 270 
n y дА Or ÔroÀ 
d D 2 T 
[ar rto, => 


lorsque À,, — À, 


[ SE E "um | 


Га (о ( o дА Or Oro 

М rr(®,(r)) = 2A. 

Avec les fonctions de Bessel de premiére et deuxiéme espéce Фу) =a] (Ar) c bY sr) nous 
trouvons : 

P) _ —r(aJ(A,r) -bY(A,r)) et CL (à (aJ Ar) BEP) 

— 3 909,00. 00,0) 


QA Or 


BE? д ФФ, (r) = Ô б 2.0) =-rX Ф (r) équation de Bessel 
ôr дА Or Or 
д'Ф„() = —rÀ, b, (r) 
Oro? 
2 E D : | 
(бөө е] (ре) een] 
[ dz r (B, (r)) = 24. = 22, 
r? (207 +D (r)? | 

(ei We il | 

=> [4 r CAO) = 2 | 


Ce qui est bien le résultat attendu. 
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Calculons ce que serait la norme si les fonctions propres étaient de la forme: 
b (r)=aJ (A„r)+bY (Ar) 


, où p est un réel quelconque et À, valeur propre tel que les fonctions 


oo 
pe "cm 2. (ar? —р?)Ф,, (r)=0 et 
r 


retrouvons une formule analogue avec les calculs suivants, en procédant par le même passage à la 
limite : 


propres soient solutions de l'équation de Bessel el 
дг 


On obtient 
М 0o, (r) Ф, „(r) LD. mu 0*6, (r) 
" дА ôr T OrOA 
[arr (o, (0) = F 2 Ф,„ (r) - aJ, (Ar) - bY,(À,r) 
ОФ ОФ ob 
> А, P e) y —? e) Equation de Bessel >r Cl, np e + (22 _ p°)®, (r)= 0 
0A ôr ôr ôr d 
= À Ô 0o, (r) = Ô 99.50) E SC dw (f) (ё "ele, (и) 
"ôr дА дг òr r OrOA ¿r 


n n 


" ob. (r) 2 | (27° —р?) n 
Tate 


[ "dr r(o, (r) = =: 
2 ó, , ; 2(42,.2 2 Í 2 oo, , | 2 : j 
l E +Ф, (o) ir - р ) [ (S2) *9o, ,(r) 9) 
= [4 r (o, vc = TE Dto. ; n 


Ce résultat coincide bien avec une formule connue des intégrales indéfinies sur les fonctions de 
Bessel : : 


[Га rU, r)} = г, an CCD) 


2 
AE AD À 222 Ae 
[arr an} = ou EX Dar =p Jj rA Jur) = pJ Aur) And, (Aur) 
Клад, а) o aor sf 
i 24? 
bitt, оо) 2р0 Ale nU VD. 
Е ой? 
I er, An) +J „u AABI, UD) 2 PAC) — e 
242 ecurrence pil Vn = CR p Vn p Vn 
Ракаў = LG, EE -p yr S LATIO) CHAP DE "m 


n 
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On a démontré que pour un système de fonctions propres orthogonales satisfaisant à l'équation de 
Bessel, la norme se calculait comme suit : 


eo 
Equation de Bessel > r - E = 2) +( er? - p)o, , (r) 20 
r 


Or 


OG, (r) Í p? i 
ДЕЕ 


2 


Роиг ип ie orthogonal de fonctions propres de Bessel modifiée et d'ordre purement 


Jo, e = [ar ro, 0) = 


imaginaire, Qu) ed PA) DE G dar): l'équation de Bessel modifiée devient : 
On note o, (r) 7 D„; (ir) 


Equation de Bessel modifiée d'ordre imaginaire 


GI 6ob Gir) 
dee Ë Se Lie no, =o 


0| OO, (r) 


SCA -(p?- 222 D = 0 
= r Es | (р 272) "mU 


> dëi. "cke — 2), (r) = 0 


L'orthogonalité se démontre comme suit : 

(1) r(r$, Hlp = 2°), „ (r) =Ü x o, (r) 
D r$, hp 2) Фо) 0 хф, (о) 
C.L. 


oa, (lr b ()70 8,0) 26 (n) aj0,,(5)* b, (n) 202 0, ,'() => Zil 


n,p n,p 


o, (n) + b, 


s, '()=0= Cm ph) = Er a6, 0) 55,05) 0 8, (= 28, (2) 
S иф 2 KOC OLII -2⁄2)@ (08, -reS 0), 0) - (p -Ar )ф, FOUR MOL 0 
er SE — A LO (ld, (r) +r[r@, On, (r) -r(, cl, (r) =0 
S r (22-22 Y, (Ф, 0) (Ф, 0, 0-9, (Ф, 0) -0 
intégration sur [гт] = (42 – X fa r O, (r)5, , (r)+ [r ©, AD: p- D, NOU „© =0 
He. (Ф, 00-9, Ф, O 
(2-2) 


et comme précédemment I5, 06, (r)- ФФ, cl = 


soit [а r Ф,, (Dn, (r) = 


=> (42, — A fa r Ф,, (D „ , (r) =0 
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La formule de la norme peut s'établir par passage à la limite, comme réalisé précédemment : 
posons А, = А, + AÀ dans l'équation > А, =â? AAL. + А„) 
Ho. 6,,0-6,,08,,0)] 
(2-2) 
ôD „ „(r) 
ОА 


farr Ф, (Ф, Ar: 


substituons o, (r) x D (r) 1 ju et en passant à la limite ЛА — 0 


ER ; : d Ô 
en considérant les fonctions dépendant de r et А => — > = 


r Or 


Ф od оф, (2) oð _ ë 
Ф, ac A o 22200 SC m. (z) ÔD „ (r) P (59, 
AE e ôA а PY ra, 
r ôD „ „(z) ôD „ (r) - Du. пб Ф") 
8.) ôr | EEN 
dr D r GR D 
Í, Á Яш ) (2, " n) 
lorsque À, — À, 
r Ô, (r) CO, (r) = OO, p (r) OO, (r) 
n P ) d Oro дА ôr 
d D = п 
[arr €, „(0 2 
Si d. (r) = aJ, (UA) + bY, (ir) 
ER ap ôd _ 
> À, — у) =r— 9) Equation de Bessel modifiée — r 2 r _ P (r) + (p? = Rr), (r) =0 
дА ôr ol ôr > 
n Or QA Or Or r ESCH Ar 


Е (22 — p?) $ ob, (r) ali 
: (е0) Ar Al & 
=> e r (@,, (ө) => à, + 
| jJ (= 2 p 0o, (r) A. 
o etes] 
=> [ar r (@,, (9) SE 1 


п 
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Reportons cette expression avec une fonctions propre de la forme J;,(iA,r) pour confirmer une autre 
expression connue de la méme intégrale indéfinie , soit : 


LG aJ =J; GA DI, GADJ 
rf = " 


à démontrer Í "dr r Fe (i À, 
т 


2 
Ф, Achs J,, GA) 
So? ~ 2 
OD, r) ASi Ar) =. v. OO, (r) ge 
An | Z | mo ADF 
2 
Ë о, б! = m | EMT 222 
[arr GA) = x 
Dérivée  irA Jip (iA,r) =i pJ; Ar) -CiA rJ; GA, r) => rA J, GAP) pJiGA,r) - ALPJ; А) 
| p p il 
di (d > (i А9! "n rp? | + É J;, G À,r) E Jipa (i Z) 
E 2 
n 2 2 e 
dl JG UL 3 => | + v (J GA) (аА) 22Р, GA), a G an) 
H 2 
| | 2 | | І 
EUX Af (0а - ZF J. (AJ pui 2) 
8 2 


2 
Récurrence J; , (iA, r) = т Ji GA, r) - Ji 4GA,r) 


kb. eA Dy Ae, GA 


=> "drr (J. AD) = ; 


c.q.f .d 
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Exemple : Cylindre plein (r,z) soumis à des conditions aux limites de Dirichlet homogènes en z 


Sans restreindre le problème on peut fixer le rayon et la hauteur du cylindre à 1. 


2 2 
д T(r,z) | 1 OT(r.z) | O^T(r,z) E 


or? r Or ez? i 
T(r,2)|_, fini 
Тар, =] 
T(r,z) ,-0 
T(r,z) ,-0 


La condition homogéne en z fixe le choix des fonctions propres sinusoidal en z. La solution en r est 
une combinaison linéaire des fonctions de Bessel modifiée. Seule la fonction I est admissible car sa 
valeur est finie en r=0, la fonction K, quant à elle, diverge en r=0. 


Il vient en prenant en compte toute les conditions aux limites énoncées plus haut : 
2,(2) = A,Cos(A,z) * B,Sin(A,z) 

R. (r) = (Ar) 

Z,(0)20— 4,20 

Z.) =0= 1, = пл 

T(r,z)= Y` B,Sin(nz z)I, (nza r) 


À, #0,0 
Т(1,2) =1= > B,Sin(nz 2)1,(пл) 
А, 50,00 


L'orthogonalité des fonctions propres appliquée sur les conditions aux limites еп г et notamment la 
dernière, conduit à la détermination des coefficients : 


e Sin(nz z) I 21-(-1") 


nz Iy(nz) 


B = 


n 


1 
I (nz) f dz Sin’ (nz z) 
0 


_4 : Iy((2n *1)z r) 
T(r, z) = = D Sin((2n+1)x z) Qn + DI (Qn + Dz) 
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Pour le problème plus général 


2 2 
Ô T(r.z) | 1 OT(r:2) | O'T(r,z) 8 


er? r =ò dz? 9 
T(r,z) , fini 
T(r,z),, =f,(@) 
T(r,2)., = 0 
T(r,z), , =0 


La solution s'écrit sous la forme 


L 
А, = e je Sin (A, )=-= 
ГА 
2 Í dz f. (z) Sin(À, z) 
=" 0 
m LI) 


C, = [dz f. Sin(A, z) 


2 C 
T(r, z) = r" TA r) Sin(À, z 
( 2 l, Z Ia d í ) И 
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Exemple : Cylindre creux (r,z) soumis à des conditions aux limites de Dirichlet homogènes en z 


Sans restreindre le problème on peut fixer le rayon intérieur à 1 et extérieur à k et la hauteur du 
cylindre à 1. 


2 2 
д T(r.z) | 1 OT(r:2) | Ô T(r,z) Е 


or? r =ò дг? | 
Tz) ., =0 
tech, =] 
T(r,z), , =0 
T(r,z), , = 0 


La condition homogène en z fixe le choix des fonctions propres sinusoïdal en z. La solution en r est 
une combinaison linéaire des fonctions de Bessel modifiée. Cette fois-ci les deux fonctions | et K 
sont admissibles puisque le cylindre est creux. 


Z (z) = A,Cos(A,z)  B,Sin(A, 2) 
R (т) = С L(A,r)+ D K, (A,r) 
Z,(0)=0—= 4, =0 
Z,(1) 202 4, = пл 
L) 


R(1202CA(AL)-DK,(4,))202 D=-C——— 
„ (1) "i „) 0 ( „) K,(,) 


T(r,z) = > B, Sin(nz 200 . | 


T(kz)-1- У, B,Sin(nz 109 Ge 


L'orthogonalité des fonctions propres appliquée sur les conditions aux limites en r et notamment la 
dernière, conduit à la détermination des coefficients : 


i үй Sin(nz z) Š (- cay) 


B, 7 = 
Е (2,0) Ehl | [ae Séiss 2) nz E GA KE) | 
DA) КФ һы) Ko) 


А, = (20а + л 


(222 s K, (Ar) | 

I Sin(A, z) VI) KG) 

T(r, z) m e À, I, (A, Kk) = K, (A, K) 
LA) K,@,) 
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Pour le problème plus général 


2 2 
Ô T(r.z), 1 OT(r.z) | д uo. -0 


or? r дг 02 
T(r,z) "s 0 
Tez), =/,(® 
T(r,2)., = 0 
T(r,z), , =0 


La solution s'écrit sous la forme : 


1, 
[az Sin° (A, z)= 2 
d à 


[а 7.02) Sin(À, z) 


2( dz f. (2) Sin(À, z) 


B 0 — 0 
н Gell S rd E Gel 
0 


1,0411) Ko) 
C, = [dz f. Sin(, z) 


LAr) 


I, (4,1,1) K, (4,1,1) 


T(r,z) Е > C, Sin(à, 2) pas 


z n=1l,+0 


Kr) 
KA) 


À, E (4,1,5) 
I, (4,1,1) 


К, (А) 
K, (4,1,1) 
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Exemple : Cylindre creux (r,z) soumis à des conditions aux limites de Dirichlet homogènes en z 


Sans restreindre le problème on peut fixer le rayon intérieur à 1 et extérieur à k et la hauteur du 
cylindre à 1. 


2 2 
д T(r.z) | 1 OT(r:2) | Ô T(r,z) Е 


or? r =ò дг? | 
T(r,2)|_ =1 
tesch = U 
T(r,z), , = 0 
T(r,z), , = 0 


La solution est calculée de la même manière que précédemment et donne le résultat suivant : 


А, =(2n+Dr 


|200 њал 
y Sin(A, 2) 1(4,0) | KE) 
Oe. À. ` (DO) Ka) 

1,0,0) Ky) 


Pour le probléme plus général 


2 2 
Ô T(r,z), 1 SD O T(r,z) а 


or? r дг oz? : 
Т, = f.) 
Trz), =0 
T(r,z), , = 0 
T(r,z) _, =0 
La solution s'écrit sous la forme : 
of dz f.(z) Sin(A, z) А 
dt qeu C, = [dz f.) Sin(A, z) 
l dE 
1,041) К) 


Eu 
т.222. У с, Sin(A, z) 104,2) КМ) 
L Z A, (LA) EEN 
Hd E 


n'r2 
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Exemple : Cylindre plein (r,z) soumis à des conditions aux limites mixte de Dirichlet et Neumann 
homogènes en z 


Sans restreindre le problème on peut fixer le rayon et la hauteur du cylindre à 1. 


2 2 
д T(r,z) | 1 OT(r.z) | O'T(r,z) S 


or? r Or дг? š 
T(r,2)|_, fini 
T(r,z)| _, =] 
T(r,z)| _, =0 
T;(r,z)) =0 


La condition homogène en z fixe le choix des fonctions propres sinusoïdal en z. La solution en r est 
une combinaison linéaire des fonctions de Bessel modifiée. Seule la fonction I est admissible car sa 
valeur est finie en r=0, la fonction K quant à elle diverge en r=0 : 


Z (z) = A,Cos(A,z)  B,Sin(A, 2) 
R, (r) =T (A,r) 
Z,(0)20— 4,=0 


Z (202 Cos(4,) 202 À, _ DE 
PE Y B sin 2" +l) x ze +1)7 5 
2,#0,00 2 
T(Lz)-1- Y B,Sin(A, z1,(À,) 
А, 0,00 


L'orthogonalité des fonctions propres appliquée sur les conditions aux limites en r et notamment la 
derniére, conduit à la détermination des coefficients : 


1 
dg сэл [а Sin(à, 2) = Zu — Соз(А„)) Cos(A,)20 
0 n 
1 
[az Sin(A, 2) 
d 2 
ÉD | ALG) 
1Ь(%„)| SQ, z) ^" ° 
0 
T(r,2)=2 $ Sin(ÀA, 2) 1,(A, r) 


n=0,+00 À 1, (А, ) 


n 
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Si les conditions aux limites sont inversées sur z 


2 2 
Ô T(r.z) | 1 OT(r.z) | O'T(r,z) 8 


0 
or? r o Oz? 
T(r,z) , fini 
T(r,z) , =] 


T.(r,z)) „=0 
T(r,2)., = 0 


Il est aisée de voir, moyennant un changement de variable z en 1-z, que la solution est : 


, _ Qn+Dz 


n 


ze l-z 


mai У SE z)) Le : 


n 


T(r,z) 22 Y (=1)' CEST = — 


On retrouve cette solution en appliquant les conditions aux limites et l'orthogonalité comme ceci 


Z,(z)= A,Cos(A,z) + B,Sin(A,z) 
R,(r)=1,(,r) 

Z',(0)=0= B, =0 

Z ()=0= Cos(4,) -0— À, SEMI 
T(r,z) - > A Cos(À, 2)1,(, r) 


2,#0,00 


T(Lz)-1- Y A Cos(À, 2)1(А„) 


À #0,00 
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L'orthogonalité des fonctions propres appliquée sur les conditions aux limites en r et notamment la 
dernière, conduit à la détermination des coefficients : 


1 | | 
А, = CQn+l)a => [а Cos(A, z) = 2 [sina,2)], = Sin(A, ) 2 (1) 
2 0 À À. 1 
1 
[ 4 Сой, 2) 
А = 0 2(-1) 
| | 2, T (A 
WEE Cos" (A, Z) n ok 2) 
0 
ке аи 
2 
(—1)” Соз(А„ z) I(À, r) 
Т H =2 n b 
зл” 2) z. à. TE 


Si les deux conditions aux limites de Neumann axiales Z sont nulles (adiabatiques dans le cas de 
conduction thermique, voir plus loin pour le détail) : 


2 2 
Ô T(r.z) | 1 OT(r.z) | O^T(r,z) = 


or? r Or дг? j 
T(r,2)|_, fini 
T(r,z)| =] 
T(r,z) „=0 
T;(r.z)), | -0 


On trouve une solution particuliére immédiate du probléme en utilisant la valeur propre nulle de 


А . ino AE (Руш И : 
l'équation séparée : 13202) puisque le terme en In(r) donne une divergence en r=0. Avec de 


tels conditions aux limites homogènes еп z, si le profil de la condition aux limites en r-1 est une 
fonction quelconque, il y a lieu d'envisager toutes les fonctions propres y compris celle de valeur 
nulle (soit une constante en l'occurrence puisque la condition de finitude nous 


Pour les problémes plus généraux : 


2 2 2 2 
Ô T(r,z), 1 ôT(r.z) | O T(r,z) _ д T(r,z), 1 OT(r,z) O'T(r,z) Е 


0 0 
др? r дг oz? ôr? r дг oz? 
Trz) , fini Trz) fini 
T(r.d),.; =f.@ Тг] | = 7.0) 
T(r,2).., = 0 T; (r,z) „E 0 


GE = 0 Trz) =0 
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Les deux solutions s'écrivent sous la forme, pour la première : 
L 
quet z [az Sin? (Q, 2) = а 
21, o 2 
L 
2 | dz f. (z) Sin(À, 2) š 
B, = —— >C, = | dz f.(z) Sin(À, z 
E Se „ = | f.) Sin(À, 2) 
Tiras Y € NON) ыд 
L n=0 ,+00 PA (4, L.) 
Et pour la deuxième : 
L 
A, IT fa cos (a, 2) 6 
21, d 2 
L 
2 [az f.(2) Cos(à, z) 1 
0 


B = = С = |а Cos(A 
; БЕ; „ = | dz f.G) Сов, 2) 


Fast У С, LATE a 
L n=0 ,+00 I; (А, n) 


Pour le probléme plus généraux avec section haute et basse "adiabatique" (dérivée nulle) 


2 2 
д T(r.z) | 1 OT(r.z) | O'T(r,z) d 


0 
or? r дг oz? 
T(r,z) , fini 
Тг] =f.@) 
Тг] _, =0 


Cl, =a) 
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II faut ajouter la solution de valeur propre nulle 

2002) = С2+р 

Z (z) = A,Cos(A,z) + B,Sin(A,z) 

R (r)= Aln(r) +B => T(r,z) fini en r 0 R, (r) = B 
R. (r) =r) 

Z', (0) =0= B, =0= Z,(z)= A,Cos(A,z) 

Z' (1) =0= Sin(A,L)=0= À, = CS 
Z',(0)=0æ7,()=0—=C=0 
T(r,2)=D+ Y A Cos(A, 2)1,(À, r) 


А„=0,® 
T(L,z) - f.(2)=D+ > ACos(À, 2)1,(4,,) 
À, 30,00 


orthogonalité des fonctions propres, n > 0 => 


Z Í dz f. (z) Cos(A, z) 


[œ us 


D = et À = 


l " L 1002,1.) 


2 


Les valeurs de C et D dépendent du profil de la fonction aux limites f.(z) : 


f.@)=1 


L 
= D=1 et А, с |а Соз(А„ 2) < Sin(4, L) -0 Ty (r, z) =1 
0 
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Pour une autre fonction limite : 


f2)=1-0(z-k) à,-77 
2 f 2|Sinz); _ 2 Sin(A,k) _ 2 Sin(A,k) 
A, 2 —  — [dz Соз(А, z)- = 1 = = 
L IL) ALIA) ALLG) nz IA) 


k k 2 Sin(2,k 
D= "m Tum uU > < Se m Cos(A, 2)1,(А„ r) 


Si la fonction limite est f.(z)-z 


faz ie 


Z 


1, 

[а z Cos, 2) = = [z Sin(z)+ Cos) |^ = => (Соз(А/.) 1)= n = (Sin) 
0 Я Я п п п 
À, = — — | z Cos(A, 2) 
l, Iy(A,,) 0 


"m zdz 
> 4. 0,5. 177-5 „z et D=2 её 
T (2n + 1) I, (ul) L 2 
m (2n 41) x 
і. 
L 41 1 


Jar, zl < > > 


Cos(À, z)1,(4 
л? Ze, (2n+1) L.A) os(A, 2)Ь(А„ r) 
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Si la fonction limite est f.(z)=z° 
f(nez A= = AL-nmz Sin(A,1.)=0 Cos(A,1.)=(=1)" 


2 


Í dz z Sin(z) = [Sin(z) -z Cos(z)] 
[а z Cos(z)= z Sin(z)+ Cos(z)] 
Vé z^ Cos(z) = (z E Sin(z) +2 z Cos(z)] 


j dz z? Cos(À, 2-3 le — 2) Sin(z)+2 z Cos(z)["* = 2 LG cata) das QT 
A == 2 77 z z? Cos(À, z) 
" LIS 
[zaz к , И 
Е Е Е = 
L 3 LA, HAL) n° x LG) 
E (-1} ns 
T(r, z) = SE 2 LD Cos(i, z)I,(A,r) А, = Л 
Si la fonction limite est f.(z)-z/k sur [0,k] et }.(2)=(12-2)/(12-К) sur [k,lz] 
2 (L-z) nz | п 
== nuire D'E = =(-1 
f.(z) elt ei f.(2) m eel А, 1 Sin(A,L)=1 Cos(A,L) = (-1) 
fa a (=a) Соѕ(А, z) = 
J k -k) 
; Ak L _ 1 d 4,1, 
=p lF Sin(z) + Cos(z)];" + DEE > — [sinok A-B ` [z Sin(z) Cos) 


=> (A,k Sin(A,k)+ Cos(A,k) -1) Z Sin(2,k) 


A А, (L -k) n 2 o -k) ( 1)” Ak Sin(A, k) = Cos(2,k)) 


| п 
-zg 0560. k)-1,+k-k(-1")= zd a (0090.0 - D -Cn ) 
n 
A, = ГІ leo Со5(А, z) = -OL nan CAD- ) 
21, I ep 
NETT C WA ien) 
© (L -2) £ (2 , А L-k 
"m E a. ded lc LL | 
I L i L i L B 2 
A = Dt: 
n = L 
_1, 2 ( 1.(Cos(a,k)-1)+k(1-(-1")) 
V< 2 “HE - k) л? 2. n? LA) Cos(2, 2)1(А„ r) 


Exemple : Cylindre creux (r,z) soumis à des conditions aux limites de Dirichlet homogènes en r 


Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques - p 111 


Sans restreindre le problème on peut fixer le rayon intérieur à 1, extérieur à k et la hauteur du 
cylindre à 1. 


2 2 
Ô T(r,z) 1 ÔT(r.z) | Ô T(r,z) а 


0 
or? r дг дг? 
T(r,z), , =0 
T(r,z)| , =0 
T(r,z), , =0 
T(r,z), , 71 


La condition homogéne en r fixe le choix des fonctions propres de Bessel en r. La solution en z est 
une combinaison linéaire des fonctions hyperboliques sinusoidales (exponentielles). Les deux 
fonctions de Bessel de premiére et deuxiéme espéce sont admissibles car il n'y a pas de conditions 
de finitude de la solution en г=0: 


R.() =C„J,(A,r)+D,Y,(A,r) 

Rf 202 €,4,0,) DEI 

AO LO 
DA) FO 
JA LG 

LY RAJ 


R (k) 0 C,J Gk) = -D,Y (Ak) > 


valeurs propres À, solution de 


Z (0)=0= A, =0 
Z. (z) = B,Sinh(A,z) 


Yr) 


T(r,z)- У, EE n avec O,() e TR See ou D, = 27 


T(rD=1= > B„Sinh(2,)D (r) 


2, #0,00 


Y,G,) 
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L'orthogonalité des fonctions propres Ф (г) appliquée sur les conditions aux limites en z et 
notamment la dernière, conduit à la détermination des coefficients : 


JA) XO Dë , Tote. ten jT 
AT AD nO) SO el Ce di e E J 
„бы? т) _ gut) X) ЛО) 
JA XA А, 0,0) JA 
OA) XQ JA 9,0 YQ) ZA) 

Ж. ` TUR JD A LUA AO 


n 


b, (| _ Dr uo 104 оор _ á E _ Gei | Y(A) A4) ) 
: " dh AA) JolA,k)) 20104,6) JA 


P,(r)= 


Jy G)- Jo) K'E) = C) 


D (1) 2 6,(k) =0 


` 1 
| rdr J,(A,„r) = 
` 1 
| ЕСЫ 
k 
[rar Ф (r) = 1 kJ (Ak) - J,(4,) _ KY (A,k) - Y(A,) 
I 4, Ј,(4,0) Y,(A,k) 


(Сла Y (AA) | X) 

Alan ah) (AGO. XO 

2(( 40,0. EG) (ZA) EQ 
Jrdr e, 0 А, СЕ 24 p KH 


o. sinna | [| AGE) оо (ла) AT 
"U (LA) OD (OD. LO 


Л - (QD) 


n 


BIELESER 


1 


= Sinh(A,) 


n 


B 2 


TEE E WOR || 
LA AH) (JUD. KA) 


Sinh(À, d'A Ju 
T(r,z)=2 > Jk) Y,(A,k) 


à Sinh( Ze JO _ dul JG) _ Y) ) 
LAN X9) (AO. XO 
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дп) KAM) WAD Jä. AG.) É A,r) Hd 
Remarquant que : 1 Jo (An) Ү,(4,) Ji (Ak) Y, (1„k) Jo (A, Kk) J (24) Ү,(4,) 
les fonctions propres peuvent également être sous cette nouvelle forme : 
k 2 D 2 À 
a, iq 00) „АЮ р of- TEE Ë D Ф, ex 
Jo (А, ) Y, (А, ) 1 2 À, 1 
Jo (A,r) Y, (À,r) => Ф', (r) - Y(A,r) J, (A,r) 
Jo (4,) Y, (4) À, Ү,(4,) Jo (4,) 
$',(D Al AA) Ф", (к) GEI AE 
À, Y, (4„) Jo (4,) À, Y Jo (,) 


СЕ Gel E say] 
oA) A) Ү (А) JAJ 


Ô (r) = 


J=- 100) = Bir 


$,(1) - 6,(k) -0 


k 


[rar 44,0) = => [19 (r) = = (GD - J (4,)) 

[rar 1,0, = Д = AK A) 

le ES Sn ed 

Al A) LO. 

(лоо) XO) (ло) XQ.) 

AUTO): d) (IO: TA) 
` 21 [лоо KA) (JA) XQ) 
pem ALES mE "o Ta) 


Bé ee r (R = 
n Si h À 2 2 2 
A, Jo, n) E Gel (262 a) | 


2 
¿Lema 200 Gre Jä Sall 
JA)  X(A,) IA) XA.) 


Sinh(À, {бә Had 
ES A) KA) 


“re asiman HAGR кеа 09 Sall 
JU. nA) GQ) 3,05) 
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Pour le problème plus général : 
OT(r,2) , Lëtz), Ô T(r,z) _ 


or? r Or дг? : 
T(r,z) xm 0 
T(r,z) E 0 
T(r,z)| _, = Ü 
Тг] = 0) 


Ji ye VU) 
J (ALD Y (ALl) 
J G¿,r) | Y,G,r) 
TA) XL) 


À, valeurs propres solution de 


Ar) _ Yr) 
Jo (4,1, ) Ү, (4,1, ) 


tee X: вто) 


А, 0,00 


li 
a (r) = JD n) аар coe] 
AO TAL) ` 2 
AHA) А) 
HG) UD 
9.0) YO JALI 9,0) XQ). JL) 
ы AL) CAL À HA) 4012 


n n 


2 2 
Фр A) _ ASI) L(XQL) AQ) 
2 LO) 40,3). 2404. HOL) 


| ауес Ф (r) = 


D'„ (r) = Ју) =-J;(r) KE) = =Y (r) 


D, (la) - 6,(,,) - 0 


im П 


1 
r2 r2 B B 
B, = |rdr f (r) J (Ar) В, = |rdr f (r) Y, (A„r) |rdr f.(r)®,(r)= e D 
„=J | i | à | FA) LOL) 
k 
B B 
rdr f,(r)®,(r) Ú » 
= | => J (Alo) Y, (2„l,.) 
n ] 2 2 2 
Sinh ALB, C swa Y 2 [4642 пола) (Ala) Ka 
TR 2 J, (4,/,) Ү,(04,/,,) 2 Jo (4,/.,) Y(A.L5) 
| В» By, | 
B = Jo (¿,1.,) Ү,(4,/,,) 


п 


2 2 
Sinh(A,L) E Tad d JUL) mE 
Jo (A15) Y (Al.,) Jo EM Y (A15) 


t Se z 2 Jun, “| aan | 
T(r,z)=2 Kä 004,2) Yolla) 004,12) YO) 


2 2 
WEE 


Jo (А) Y, (A,1,,) Ji LA) Ү (Al...) 


Jo (4,1,) =) Y, (4,1,) 
Jo (4,1,1) Y (Al) 


avec А, valeurs propres solution de 


Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques - p 115 


Lorsque la condition aux limites est constante alors : 
š 
r2 1 1 
B, = frar Jo (A,r) = А? LA (п) SS ER Rach (4,1,5) = I, (4,1,1 )) 


l4 n n 


1,2 


1 1 
By, = | rdr KAD ett = COLE Aln) 
Li n n 
I p (ry = Le Uu La Qu) _ CH Cubes) YJ (A.D) 
Li | À, Jo (4,1, ) Y, (4,1, ) 


_1(, лал) XO) 4 (лалә) OU 
EE TALI VIAL) YL) 


(Ла) 11) | ( AQ) GL) 
A. Jal. BAL) LS) FALS) 


2 2 
sica) (4642 EG) E Y (A, y 


n 


Jo (4,1,) Ү,(4,/,) Jo (4,/,) Ү,(4,/,) 
2 


В, = 
à sma, [1 (A5) (4,1,2 y JG) YU) J 
JS) 0,1,5) J (AL) UE) 
зин, = 7 ET 
Ehe o( ula) al DIE 


n=1,+00 


JA) sa), у 0) 


À, Sinh(A,L.) 2 
Jo (4,1,) Ү,(4,1,) Jo (4,/,) ON 
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Exemple : Cylindre creux (r,z) soumis à des conditions aux limites de Dirichlet homogènes en r 


Sans restreindre le problème on peut fixer le rayon intérieur à 1, extérieur à k et la hauteur du 
cylindre à 1. 


2 2 
Ô T(r,z) loT(nz), O'T(r,z) E 


дг? r Or дг? | 
rej. =0 
T(r,z), , =0 
T2): = 0 
Lech =] 


La condition homogène en r fixe le choix des fonctions propres de Bessel en r. La solution en z est 
une combinaison linéaire des fonctions hyperboliques sinusoïdales (exponentielles). Les deux 
fonctions de Bessel de première et deuxième espèce sont admissibles car il n'y a pas de condition 
de finitude de la solution en г=0: 


К, (7) = CJ, sr) * D,Y Ar) 

R,' (120 C,J,(4,) 2 -D,Y (À,) 

JE) YA) 
JA) Y(A) 
Jk) XA) 

JA) X) 


R (k) 202 4,C, J (4,k) = -À,D, Y, (À,k) > 


valeurs propres À, solution de 


Z,(0)=0= 4 =0 

2,(2) = B,Sinh(A,z) 

T(,2)= > вто) Gr) SN avec Ф (r) = J Ur) A 
Йе JE) Ү,(А,К) JE)  Y,G,K) 

T(rÀ) -1- > B„Sinh(2,)0 (r) 


An #0,00 
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L'orthogonalité des fonctions propres Ф (г) appliquée sur les conditions aux limites en z et 


notamment la dernière, conduit à la détermination des coefficients : 


à i D LOS ko, f=] war S| D LJ 


ЛА) Jä. 2( 4 
JO (r) kär 
Jk) 04,0) 
AY Anr) _ 4„J sr) 
Y,U,k) ` J Uk) 


o (0-0 9500. LAH JOD o m=- d Y.(À,) SEH 
À, (4,0) JAk) Y (Ak) JAP 


n 


P,(r)= 


pue 


Ју) = (и) 000) = (0) 


HA. RA) 
OJQA YA, 


EE Goal 1 Ju) na) | 


c, (1) 


Y, (U,k) (4,0) Jk) 104,0) 


[г 4,0, = ESA 
[rar Y (A„r) = => [ro = = (КҮ, (4,k)- Y, (2,)) 


[rar o, c) - 1210200) 0,0 TG) 
СА AGO Y, (4k) 


_ 11 049 YG) | ү Лб) _ Т 
A LAG) XOU) LAG 10,0) 


[raro VE 
САЛ, XO) 


[raro J (205 mA 
NE. JA) Y, 
` Sinh(À,) 


Ф g sa [er ар) | A). X) J 


n 


Y, U,k) (4,0) Jk) YE) 


Sinh(A en TUE 
T(2)-2k Y JA GD (АЕ) XE 


"бле ima, [er Goal Es SH 
XA AD) (JA) YO) 
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Pour le cas plus général : 
OT(r,z) 1ÔT(r,z) O'T(r,z) 
2 + + 2 
Or r дг Oz 
Т2) = 


= 0 


T(r,2),., ef 
Tr, z) = 
REO 


La solution s'écrit : 


O, (|° = [rar Ф (r) = ВЕ но) 


AO) От) Aua. go) rer Ô 9200 HAN _ 40D 
HAE) KE) `” š i A YA.) AGL) 
oun Palis) YG) ЛО) 9,0) YQ) HAL) _ 
SC д HAL) SULLO. À Tolla) Aide 


JO) Y) 1 (yl - ea) A042] EA XQ) TO) 
(А) PU) Z GU FAL) Pl ids 


nra 
Ву, 
J "m 13 Y (A15) 


RCE SH 
J 104,11) Y (A, La) 


D (r) = 


Ф„(1)= Ké 


В, = aer conan B, = [arr tto s 


Vua asa (r) - Sé , re 


sein J) LI HA) л) Ж EE SCR 
Yi) ЈА) 2 (HA) KA) 


2 É Ee ium a. ë Sinh(A, d EE К Wd 
T(r,z) = 25 OX^"n^r2 0 "n° r2 0 "n° r2 On" r2 


нао лол) „(мә н) 


п 


Ү,(4,/,) Jo.) J (A, L, ) Y, (Al) 
J, (À, 12) _ Y, (À, L3) 


À, tq 
J, (4, La) Y (Z, la) 
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Exemple : Cylindre creux (r,z) soumis à des conditions aux limites de Dirichlet homogènes en r 


Sans restreindre le problème on peut fixer le rayon intérieur à 1, extérieur à k et la hauteur du 
cylindre à 1. 


2 2 
д T(r,z) | 1 OT(r:2) | Ô д -0 


др? r дг Oz 
Tür zy a =0 
T.(r,z) 470 
T(r,z) ,-0 
T(r,z), ,-1 


La condition homogéne en r fixe le choix des fonctions propres de Bessel en r. La solution en z est 
une combinaison linéaire des fonctions hyperboliques sinusoidales (exponentielles). Les deux 
fonctions de Bessel de premiére et deuxiéme espéce sont admissibles car il n'y a pas de condition 
de finitude de la solution en г=0: 


Ries C, Ja (A,r) * D,Y Ar) 

R.) =0= C,J,(4,) = -D,Y,(4,) 

JE) Ү(А&) 
TA) X) 
AG, XQ, 

JA) X) 


R (к) 302 À,C,J,(A,k) = -A,D,Y (A, k) > 


valeurs propres À, solution de 


Z,(0)20— A, =0 
Z,(z) = B,Sinh(A,z) 

g | JG К,Р) JA) Xr) 
MAS à. взто,2{ LA) X) | a £ UE 
T(r1)=1= > B,Sinh(A,)O, (r) 


An #0,00 


Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques 


— p 120 


L'orthogonalité des fonctions propres Ф (г) appliquée sur les conditions aux limites en z et 


notamment la derniére, conduit à la détermination des coefficients : 


LA) (AO tero. d k cou] 
JG) Ou A К" 
JD X) 
JQ) KA) 
d, ti XQ) Ji 
A KA) Ji 
A (k) XO 
JQ) X) 
00 Y) AG) PMB HI AQ). 
2, YA) JA) À FU Cu) 


n n 


оү 
А, 


D, (r) = 


AGE TO AGE Hir 


$,0)-0 ®,(4)= 


Le e 


LUE Dr A JO) XY 
JA) g 2\ J 775) 
L (J (4,5) AG) 


Í rdr J (A„r)  — [rI nf" = 


A 


1 
2 e 
` 1 
rdr Y (À r)  — 
| o Gr) 2 


[rY @) = l(a,0-Y(,) 


À: 


јего = kJ (A,k)-J,(4,) KY (A,K) -Y.(A,) 
À, Jo.) LO.) 


(лао EG) (ZA) Y) 
ak EE YG) J] UG) SCH 
[rar ©, = EH 
| A, JI XO) 


k JA) Y() 
e D, (r) |299 223 


жы [0 HORT 1050). xà 
PON КЕ X09) S 


Sinh (A, ТЕ AC) [262 d 


Y (AK) 


T(r,z)- Y. YA) JU) X JU) KO) 
ЕВЕ Goal (29) Sail 


Ju) YQ) ) 244405). Y,@Q,) 


Pour le cas plus général : 


t J, (A, Kk) = 


Jo) 


Ү,(4,) 
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2 2 
O T(r,z) " 1 ôT(r, z) " O'T(r,z) 


-0 
or? r дг oz? 
T(r,z) , =0 
L(nz) , -0 
T(r,z), , = 


тз = f) 


La solution s'écrit : 


J (å, 12) _ Y (A, L) Ф'()” А 1,2 
À, tq LAL) XQ КОЈЕ = [rar o,e) = i, { — +0 (7) J 


IAr Xr) I'M At yop Ea Y(Ar) JA) 
JE Dp ROAS Se | Ao ` TUL) JOL) 
LG) 35055) 

J (À l4) Y. 0 (2, l4) 

QUU MO) hd) ` D'AL). X ay J 1) 

À Y) AG) | A Al) JO) 


EE SE 2 | b "US > N AE L) YO, Е | 
2 А, 2 J (Q, la) Y O (A, La) Jo (А, la) Y O (A, 11) 
1,2 la 1,2 В, В 


B„ = [rar FOJIA) В, = |на f) Xr). [drr 0)Ф, (0) = 


La 


Ò, (r) = 


$,(,)-0 0,0;)- 


= Yn 


Jo (4,11) Y (la) 


J (1) An) 


{arr 0,0 d DE? ) 
sica М E EA BEE m) 


J Ul) Y GL) J (1) 0,11) 


Sinh(A, d В, В,, | (Аг) | XA) | 
T(r,z) -2 У J ia I DADALO FA) 


vt тм Е a) E Ë AE L) _ YG, " | 
Jo (А, la) Y, G, La) J, (Z, La) Y, o (A, l4) 


J, (4, luo) _ Y (A, L2) 
Joada) Al) 
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Exemple : Cylindre plein (r,z) soumis à des conditions aux limites de Dirichlet homogènes en r 


Sans restreindre la généralité du problème, on peut fixer le rayon à 1 et la hauteur du cylindre à 1. 


2 2 
д T(r,z) | 1 OT(r:2) | O T(r,z) E 


or? r Or dz? i 
T(r,2)|_, fini 
Т(т,2)| _ =0 
T(r,z) ,-0 
T(r,z), ,-1 


La condition homogéne en r fixe le choix des fonctions propres de Bessel en r. La solution en z est 
une combinaison linéaire des fonctions hyperboliques sinusoidales (exponentielles). La fonction de 
Bessel de deuxième espèce n'est pas admissible la solution en r=0 doit être finie : 


R,(r) CA) + DK sr) 

R,(0) est borné > D, =0 

R,(D=0= J,(4,)=0 

valeurs propres À, solution de |. J4(A,) = 0 
Z,(0)20— 4, =0 

Z,(z) = B,Sinh(A,z) 

T(r,z)- > B Sinh(2,z)J,(A,r) avec D, (r)= J (A,r) 


2,#0,00 


T(r1)=1= > B„Sinh(2,)0 (r) 


An #0,00 


L'orthogonalité des fonctions propres Ф (r) appliquée sur les conditions aux limites en z et 


notamment la dernière, conduit à la détermination des coefficients : 


À tq J (2,)=0 ЕЕЕ 2 eeu) om 


0 


n 


D, (r) =J,(QU,„r) PD ла) Ф (1) =0 aD - A.) lo, (nf = J, v 


1 1 
[rar 9,0) = [rar J (2,7) = 1 pup 108) 
9 0 À, À, 
1 
[rar o, (r) 
a -* 9 _ 2J,(2,) Е 2 
Sinha E, — A, Sinh), (A,Y 2,Sinh(A,)2,(À,) 


T(r,z)-2 У Sinh(À, z2)J (A,r) 


2, tq J,(2,)=0 
„03 А,9іпА(А, )Ј7,(2А,) „ tq Ju) 
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Pour le problème plus général : 
OT(r,z) 1 ôT(r,z) OT(r,z) 
SE + 5 = 

Or r дг 02 

Trz) fini 

T(r.2)| ; =0 

z=0 = 0 

7e) 


La solution s'écrit : 


0 


T(r,z) 
T(r,z) 


1°ф' 1 2 LJA 
А, tq J, (2,l,) = 0 lo, (M) = F L4 „) = Cr -< „) 


nr 


B 


Jn 


` Sinh(A.L) 


nz 


B, = [dr r f.) 9,07 [dr r f; (0) 0,0) В, А 
0 0 D, (ғ) 


ni y Ba Sinh(À, DIAN) 
l, n#0,+00 Sinh(A,l.)J, (GA) 


VE 
posons r'= Ат 


А, tq Jo (4,1,) = 0 


d 1 ^t ' ' ' L, 
si f.()=1= В, EE Jur r л) 1 40, L) 


pct У Sinh(?, DIr) 
LE, À Sinh(A,1 


nz 


À, tq J, (2,1.)=0 
)J,(4,1.) n q al Së 
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Exemple : Cylindre plein (r,z) soumis à des conditions aux limites mixtes de Dirichlet et Neumann 
homogènes en r 


Sans restreindre le problème on peut fixer le rayon à 1 et la hauteur du cylindre à 1. 


2 2 
O'T(r,z) R 1 ÔT(r,z) B O'T(r,z) 


or? r Or дг? EE 
T(r,z)| , fini 
T(r.2)| _, = 
T(r,z)| _, = 
T;(r,z)) _ = 


La condition homogène en r fixe le choix des fonctions propres de Bessel en r. La solution en z est 
une combinaison linéaire des fonctions hyperboliques sinusoïdales (exponentielles). La fonction de 
Bessel de deuxième espèce n'est pas admissible la solution en r=0 doit être finie : 


R (r) = C,J (4„r)+ D,Y (sr) 

R,(0) estborné > D, =0 

R,()=0= J,(4,)=0 

valeurs propres À, solution de Ј,(А,) = 0 

Z,(0)=0= 4,=0 

Z,(z)= B,Sinh(A,z) 

T(r,z) = DE B,Sinh(A,z)J,(A,r) ауесФ (r)= J, (A,r) 


À,#0,00 
T.(rD=1= > 1,B,Cosh(2,)®,(r) 
À,#0,0 


L'orthogonalité des fonctions propres Ф (г) appliquée sur les conditions aux limites en z et 
notamment la dernière, conduit à la détermination des coefficients : 


A tq J (2,)=0 |o, GY = | war JEE 2 ° 
pue Aen Ò 50, 0,0-0 LO = 79.) buet = 28) 
[ar o, (r)= [rar Jr) => [rd (r) = = 

frar bd = 2J,(4,) _ 2 


A, Cosh(A,)J, (A, = A Cosh(4,)J,(À,) 


Sinh(2, z)J,(2„r) 
T(r,z)=2 
GE Z À, "Cosh(A, )J,(2,) 


À, tq Jo) = 0 
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Si le flux est imposé sur la section basse en z=0, alors la solution est obtenue par changement de 
variable de z en 1-z 


rr y =2 y, SG. U-DI) 
і n=0,+00 A, Cosh(A, PA (А, ) 


Pour le problème plus général : 
O°T(r,2) , Lëtz), Ô T(r,z) _ 


2 2 
Ô T(r.z) | 1 OT(r.z) | O'T(r,z) z 


0 0 
or? r дг дг? or? r дг дг? 
Trz) fini Trz) _, fini 
T(r,2).., = 0 T(r,2).., = 0 
T(r,z), , = 0 ELE „=/,@) 


Ciel, =f.() 


Les solutions du premier et deuxième problèmes s'écrivent : 


Tr, z) =0 


ou le flux en section basse 


2 r2 
A, tq 4,0,,)-0. | (0)? =D 
1, 1, 
B, = [arr 7,0) 8 „(r) = [ar r LORAN 
0 0 
B ES B, EH 2B, 
" ACosh(A,L), (o) Ayl, Cosh(A,L)J Qs) 


Probléme 1 
T(r, z) = —D > B, Sinh(A, Cr) 
1 nono A, Cosh(A,L)JV(A,l.) 
Probléme 2 z«€»1-z 
rr.) = 2 y, Визит, 0, = DII) 
Lace ¿Ü СОЗА l, ) J (A,L.) 


nr 


À, tq J GL.) m 0 


nz 
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Exemple : Cylindre plein (r,z) soumis à des conditions aux limites mixtes de Dirichlet et Robin 
homogènes en r 


Sans restreindre le problème on peut fixer le rayon à 1 et la hauteur du cylindre à 1. 


OT(r,z) 1 ôT(r,z) O'T(r,z) 
+ + = 
or? r =ò dz? 
T(r,z) , fini 
T. (r,z)* hT(r, z) Жы 0 
T(r,2) 


0 


z=0 =1 


Te. ,=0 


La condition homogène en r fixe le choix des fonctions propres de Bessel en r. La solution en z est 
une combinaison linéaire des fonctions hyperboliques sinusoïdales (exponentielles). La fonction de 
Bessel de deuxième espèce n'est pas admissible car la solution en r=0 doit être finie : 


R, (r) = C, Ja(2,r) DX sr) 
R,(0) est borné = D, =0 
R (1) - AR, (1) =0 > —4,J, (4,)+ AJ (4,)=0 
valeurs propres À, solution de À,J,(2,)=hJ,(2,) 
A, Sinh(A, ) 
` Cosh(À,) 
Z (z) = A,Cosh(A,z) + B, Sinh(A,z) = A, (Cosh(A, z)Cosh(A, ) — Sinh(À, z) Sinh(A.,)) 
Z,(z) = A,Cosh(A, (z —1)) = Cosh(A, (1— z)) 
T(r,z)- > B,Cosh(AÀ,0—-z)J,(A,r) avec Ó, (r) = J (Ar) 


DEE? 


T(r,) 21-2 > B„Cosh(2,)8 (ғ) 


А, 0,00 


Z. (1) -0— A,Sinh(2,)+ B,Cosh(À,) - 0— В, = 
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L'orthogonalité des fonctions propres Ф (г) appliquée sur les conditions aux limites en z et 


notamment la dernière, conduit à la détermination des coefficients : 


à 2K A 
ФП) 
À 


A AJ.) M.) [OO = frar o, = ВЕ © ro |] 


0 


o. (= (А) =-J (4r) 9,0) =J,(2,) 


(patr) жа 
ELE =-Л(4,) 


n n 


o Af Sek E dr il le ea?) 


n 


[rar o0) [rar а) = Zei (E = SC 
= аи z 2J,(,) 
Cosh NO А, Cosh AU 0,44.) 
Bn, O DELS » 


" A COSA DIE AN +24) Cosh(a,)J(A,)h2 A) | Cosh(A MA. IA 
À, tq À „J, (¿,) = hJ, (А,) 
T(, 2) = 2h Y Cosh(À, (1 -D)e Gu 

„71. Cosh(A, J (A, V? + А, 
Pour le probléme plus général : 
OT(r,z) lOT(r,az) O'T(r,z) 
+ + = 
or? r or Oz? 

T(r,z) , fini 


T.(r,z)* hT(r, z) Е 0 


T(r,z), „= f, 
5 


A, tq ÀJi (A,)= hJ),(A,) 


0 


À, tq ÀJ (4,/.) = hJ, LA) | 


dl = | rdr © (r) EK Oig er] 


i 2 2 
out Mike Al 


D) - J ur) Te 2-4) | 


n n 


F 


в, = [drr LOD, [dr r f) 40,7) 
| dr r f.r) Ф, (о) ids 


n 


= Совт, (A). Соз) 14 UL +) 


dg B, À. Cosh(A, 0, - z)J Aur) 
"TU, Cosh(A], Jy QUI A) 


À, tq À „J (2 l )=AJ (A l ) 


nr nr 
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Exemple : Cylindre plein (r,z) soumis à des conditions aux limites mixtes de Dirichlet et Neumann 
homogènes en r 


Sans restreindre le problème on peut fixer le rayon à 1 et la hauteur du cylindre à 1. 


2 2 
д T(r,z) | 1 OT(r.z) | O'T(r,z) S 


er? r =ò ez? : 
T(r,z), , fini 
T (а = 0 
Tr) =0 
T(r,z), , =] 


La condition homogène en r fixe le choix des fonctions propres de Bessel en r. La solution en z est 
une combinaison linéaire des fonctions hyperboliques sinusoidales (exponentielles). La fonction de 
Bessel de deuxième espèce n'est pas admissible la solution en r=0 doit être finie : 

R (т) = C,J(A,r) - D, HAr) 

R,(0) est borné > D, =0 

R,'()=0=> J,(1,)=0 valeurs propres À, solution de Ј,(А,) = 0 

Z,(0)=0— A4, =0= Z,(z) = B,Sinh(A,z) 

T(r,z)- > B„Sinh(2,2)J (A,r) avec d, (r) = Jär) 


2,#0,00 


T(rD=1= > B,Sinh(4,)®,(r) 


An #0,00 


L'orthogonalité des fonctions propres Ф (г) appliquée sur les conditions aux limites en z et 


notamment la dernière, conduit à la détermination des coefficients : 


A J,2,)=0 Jo, = j rar d BE e e] _ o. a 
POS EAN 920-0 [фу == —- 
Jo, c = ^02) 
fra D, (r) = frar J (A) = Ze IOP = чыз Se 
[rar d (r) 


0 


B,=_ >> T 0 
Cosh(A, ) b, (ғ) 
T(r,z)=0 
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Il est clair que ce ne peut être la solution recherchée tout simplement parce qu'il faut rechercher 
une de ses composantes également dans les solutions admissibles de l'équation séparée à valeur 
propre nulle. Soit tout simplement 10-2) = 2 


On aboutirait à la méme solution 1*2) =Z pour le probléme : 


2 2 
д T(r.z) | 1 OT(r.z) | д ы) E 


er? r Or 02 
T(r,z) , fini 
T, (rz), , =0 
Tr us = 0 
T,'(nz) 4 =1 


Ceci indique que dans la recherche d'une solution, il ne faut pas toujours d'emblée s'engouffrer 
dans une recherche calculatoire des coefficients des fonctions propres du problèmes aux limites. Il 
faut également considérer que ce type de probléme présente un cylindre isolé sur sa surface 
extérieure hormis sur les sections hautes et basses. Dans ce cas, la solution ne peut dépendre du 
paramètre r lorsque les conditions aux limites en z ne dépendent pas du rayon r. La solution ne 
peut qu'étre purement axiale (z). C'est le sens du choix de la valeur propre nulle de l'équation 
séparée. 
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Exemple : Cylindre plein (r,z) soumis à des conditions aux limites mixtes de Dirichlet et Neumann 
homogènes en r 


Reprenons le même exemple que précédemment et sans restreindre le problème on peut fixer le 
rayon à 1 et la hauteur du cylindre à 1. 


2 2 
Ô T(r,z), 1 ôT(r.z) | O T(r,z) 


=0 
ôr? г Or z 
Jini 
T (r, Z) 39 
T(r,z), „= 
T(r, z) =1,re[0,k]ke]0,1[ =l 


T(r,z) 


z-l,re[Kk,1]ke]0, vu" 


La condition homogène en r fixe le choix des fonctions propres de Bessel en r. La solution en z est 
une combinaison linéaire des fonctions hyperboliques sinusoïdales (exponentielles). La fonction de 
Bessel de deuxième espèce n'est pas admissible la solution en г=0 doit être finie. De plus il faut 
ajouter la solution de valeur propre nulle. En effet il faut considérer le problème symétrique en r 


O^T(r,z) 1 ôT(r,z) OT(r,z) 
2 LA + 2 
or r дг Oz 
Trz) , fini 


= 0 


T(r, Se = U 


=1 


z=1,re[0,k]ke]0,1[ 


T(r,z) 


z=1,re[k,1]ke]0,1[ 


Or si nous superposons les deux problèmes précédents nous avons la solution du problème suivant 
qui est triviale, et est à rechercher dans les solutions à valeur propre nulle de l'équation séparée : 
d T(r,z) " 1 OT(r,z) u OT(r,z) 


or? r Qr ez? Se 
T(r,z) , fini 
T (r, 2) 470 
T(r.z) _, =0 
T(r,z) _, =1 
=> T(r,z) =z 


II faut donc introduire un terme A, z+B, la condition en z=0 donne B, =0, il nous reste donc 
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R, (r) = C, JAN + D,Y, (A,r) 

R,(0) est borné = D, =0 

R,'(120J,(4,)20 valeurs propres À, solution de J (¿,)= 0 
Z,(0)20— A, 20— Z, (2) = B,Sinh(A,z) 

T(r,z)- Az > B,Sinh(A,z)J (A,r) avec D, (r) = Jy (À,r) 


An #0,00 


T(r,l) = À, + > B,Sinh(1, )®,(r)=1-9(r-k) 9(7 - EN fonction de heaviside 


А, #0,00 


L'orthogonalité des fonctions propres Ф (г) appliquée sur les conditions aux limites en z et 


notamment la derniére, conduit à la détermination des coefficients. On remarque également que 


la fonction propre de valeur nulle qui est égale ici à la fonction identité : guod est également 
orthogonale avec les fonctions propres de valeur Ф (r) , comme on le vérifie dans le calcul suivant 


1 
n —0,À, = 0 = intégrons avec la fonction Jy (Aor) => [rar J (A,r) = We = 0 
0 n 
1 k 
=> orthogonalité de Ф (r) avec Ф (r) > A [rar = [rar > A, =k? 
0 0 
n> 1,¿, + 0 intégrons avec la fonction J,(A„r) => 
1 
1 2 ' 2 
â, tq J (A,)=0 uf = [rar o, oy = EC a a |] 
0 n 0 
Ф' (r Ф' (1 
9,0)-44, 93€). gut Zell A0 
1 k 
lo, GY = E © - 4 ) 2 г [1-90 kp, (rar = [rar J, Qr) T Ды | = LU 
0 А, п 
k 
[rar Ф (r) 
7 __ 2kJ(4,k) 


À,Sinh(A,)J, (A,) 


Ti(r, zz kz +2k У J, л; эд) 
Hor A, Sinh(A, )J, (À,) 


À, tq J,(4,) = 0 
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Regardons maintenant le problème symétrique en r 


2 2 
д T(r.z) | 1 OT(r:2) | Ô T(r,z) 


=0 
ôr? г Or _ 
fini 
T als = 0 
T(r,2)., = 
T(r, z) =1,re[0,k]ke]0, I 


T(r,z) =1 


z=l,rek,1]ke]0,I[ 


Dans ce cas la solution recherchée est de la forme : 
T(r,2)= Az * > B,Sinh(A,z)J (A,r) avec D, (r) = J, (À,r) 
А, 0,00 
T(rA)-2 А, + > B„Sinh(2, )b,(r)29(r—Kk) 9 fonction de heaviside 
À, #0,00 


et la prise en compte de la dernière condition aux limites donne : 


1 
п= 0,4, = 0— int égrons avec la fonction J (Agr) => [rar J,(,r)= Z 
0 


1 1 
A, [rar = [rar = À =1-k? 
0 k 


n> 1,¿, + 0 intégrons avec la fonction J (A,r) => 
Ф 2 Ech 


lo, = : 
мт PAG pe 
ES o, (r) 
d 2 2k J,(A,K) 


n 


À,Sinh(A,)J, (À,) 
ES J (A,k) Sinh(A, 24, (Ar) 
no А„5їпА(А„)Л, (A,) 

T(r,z)+T„(r,z) =z 


T,(r,z) =(1— k°)z À, tq J(A,) 20 


Nous voyons bien que 


Remarque : on voit que l'on peut suspecter la présence d'une fonction propre à valeur nulle dans le 
développement de la solution : 

- si la fonctions propre à valeur nulle est orthogonale aux autres fonctions propres à valeurs 
positives 

- lorsque la superposition des deux problèmes aux limites complémentaires conduits à une solution 
triviale 

- cela permet enfin dans certains cas de résoudre les solutions de problèmes aux limites 
discontinues, comme c'est le cas ici. 
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Si le problème aux limites est le suivant : 
OT(r,z) lOT(rz) OT(r,z) 
2 + + 2 
Or r =ò Oz 
T(r,z), , fini 


а = 
T(r,z), , ЕС 
T(r,z) 


=0 


z=l,re[0,k]ke]0,1[ = Т, 


T(r,z) 


z=l,rek,1]ke]0,I[ =T, 
Alors la solution selon le principe de superposition est : 


T(r.z) =T, tr, zit T,(r,z) 


T(r,z) = k?z T, E (1— &?)z Т, +2k(T, -T7) > J (1„k) Sinh(À, SE 
n=0,+00 A, Sinh(A, Li (А,) 

À, tq Л(А,) = 0 

Regardons maintenant le problème suivant еп г 


OT(r,z) 1ôT(r,z) O'T(r,z) 
+ + 
or? r =ò ez? 
Trz) fini 


T'(r,2).., = 
Та). = 
T(r,z) 


=0 


=1 


z=l,rek,k, ]k,k,e]0,1[ 


T(r,z) =1 


z=l,re[ki,ks ] ks €]O,1[ 


T(r,z)- Az > B„Sinh(2,2)J (A„r) avec D,(r)=J,(â„r) 
А, 0,00 


T(r1])-2A + > B,Sinh(A, )b,(r)2 9(r—-k,)-9(r—k,) 9 fonction de heaviside 


А, #0,00 
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et la prise en compte de la dernière condition aux limites donne : 


A _, 


1 
п= 0,4, = 0 = int égrons ауес la fonction J (A, r) => [rar J (A,r) = 
0 n 
І D 
A, [rar - [rar = A, =k; -k; 
0 D 


n> 1,1, + 0 = int égrons avec la fonction J (år) => 


o.) A. 
Je 


le, c = 


k, u 
[rar 48,7 => lr Jp Ech б: kJ (4k) 


k n n 


1 


| ue E 2(k,J, (,4)-kJ,(2,k)) 


Ф (o) A Sinh(A, VO) 
T2) - (d ien? Y. (kJ (Aka) TAC) cht, 2)J,(Q,r) 
n20,+œ A,Sinh(A,)J, (A,) 
Pour le problème plus général : 
2 2 
O^T(r,z) Е 1 ÔT(r,z) u д ы) -0 


^  Sinh(A,) 


À, tq Л(А,) = 0 


or? r дг Oz 
T(r,z)| , fini 
T '(r, D = U 
T(r,z) ,-0 


TE, 2), = f, 0) 


La solution est de la forme : 

1, 1, 2 1, 1, 

n=0= A, [rar = [arr f.) => 4 =F [dr 0) = B,=[drr ti 
0 0 r 0 0 


n> 1,¿, z 0 = intégrons avec la fonction J (A,r) => 
1o, (1) 5 LJ, A) 


D, (el = 


В» c 2B,, 
Sinh(A,L)|Ib, (| 1, Sinh(A,L) J Q,1.) 


T(r,z) = St 2+ > BAC z)J (A,r) 


k, 
В, = [arr 50) Ar) В, = 
0 


„55. Sinh(A,L)Jy (AL) 


r 


| À, tq J, (4,1,) = 0 
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Exemple : Cylindre plein (r,z) soumis à des conditions aux limites mixtes de Robin homogènes en 
r 


Pour le problème général : 
OT(r,z) lOT(rz) GOT(r,z) 
„+ + 5 —— = 
Or r Cr Oz 
T(r, 2) „fini 


œT (rz) B,T(r, 2) _ = 0 


0 


a;T,(r,z)+ B,T(r,2)| =0 


oT. (rz) + BT z) = f.) 

Pour assurer la condition aux limites radiale homogéne de Robin, il vient : 
А, tq 05À,J (A...) = BJ Qi) 

Et l'expression des normes des fonctions propres devient : 


Jo, c = f rar oo = EO + ver) 


0 


Ф", (r) 2 E 2 2 
Ф) Лт) = AN EVA) + AL) 
Comme  Q,A,J,(ÀA,1.) = B.J, (2„1.) ilvient 
2 2 2 2 12 2 
rf 259211. P | ou bien |O (mj? = 22h), Аг 
2 2 В, 
Pour le respect de la condition aux limites en 2=0, il vient : 
Terme axiale x à,À, Cosh(A,z) — D,Sinh(A, 2) 
La solution se développe donc sous la forme d'une série : 
T(r,z)- >B, (z,À,,Cosh(ÀA, z) - B,Sinh(À,z))J, (Ar) 2,4a,2J(Al)=B,J,(A,1.) 


п=1,+о0 


2, 
п 2 


La dernière condition en z-l, donne : 
/,@)= УВ, (oA (0A, Sinh(A,L)— B;Cosh(2„1,))+ B. (oA, Cosh(A,L) — B,Sinh(A,L))J, (Ar) 


n=1,+00 
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Ce qui donne pour les coefficients du développement en série, la solution suivante : 
C,=(a,2,(a,,Sinh(2,L.)- B,Cosh(2,L))+ B, (o, A, Cosh(A,L ) — B,Sinh(A,L))) 


nz nz 


C, = A,Cosh(À,l. XB, -a,B,)+ Sinh(A.1.)la,a A — B.B.) 


nz nz 


2 2 2 27,2 2 
but =-= D, L- oubien Wu! = Ur di 
À, Qt B 


2 ; 2 S 
L 
[arr f 440,7) 1 
В, =" — = B, = [drr f) JAN) 
C, Ф (r)| 0 
2 B, | 
T(r,z)- T У: n (aA, Cosh(A,z) - B,Sinh(A,z) J (Ar) 
ir c Je a+ A J 
n 0, 
ou bien 
2 B, | 
T(r, z) = E »3 à 04A, Cosh(A,z) — B,Sinh(A,z) J (Ar) 


FN 
dE слона) 
В, 
ауес À, tq aA, (4,/,) = Bd; (4,1,) 
Cette expression permet d'envisager tous les cas limites possibles simulant les conditions de 
Dirichlet ou de Neumann, homogènes ou inhomogènes : 


a, =0,1 
B, =1,0 
a; =0,1 
В, =1,0 
a, =0,1 


B, = 1,0 
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Exemple : Cylindre creux (r,z) soumis à des conditions aux limites mixtes de Robin homogènes en 
r 


Pour le problème général : 
OT(r,z) 1ôT(r,z) O'T(r,z) 
+ + = 
or? r ОР дг? 
aT, (r,z)* BT(rz) , =0 


0 


eT.(r,z)* B,T(,z) , =0 
oT, (7,2) + B,T(r,z) = 0 
a,T(r..) + Bitze, = f 0) 


Pour assurer la premiére condition aux limites radiale homogéne de Robin, il vient : 

J (Ur) E (А) 
BJ Anla) - 0A JO) | BX Ul) - eA, Y (LA) : 
Pour assurer la deuxiéme condition aux limites radiale homogéne de Robin, il vient : 
В.Ј,(4,1,,)-а,4.1(4,,) _ B5 Y (9.5) - 5A, Y Cl) 
BI la) oA IOS) BA) - 0A, Y (uu) | 
Et l'expression des normes des fonctions propres devient : 


ЕЕ) 


Fonctions propres radiales ос 


À, valeurs propres solution de 


P 24 A A 
gie JA) Y (A) 
| BJ, L4) = 044,7, (AL) BY G4) - 0A, Y (Al) 
pu) = Yr) Jr) 
À, Ё BY, (4,04) a AY (A) BJ4O,L4) G A.J (ua) 
ae JUL) _ Y ln) 
BJ (Ala) E ad, (4,1,1) BY, (4,11) -aA,Y, (4,11) 
bu МАШ Y L:) 
BA) oA JIOU BYGL)- (AY OU) 
Ф'„(Ь) _ Y (Ll) J OA) 
À, BY GL) a A4,Y, A la) Bio lAl) 04,4, 0L) 
Ф'„(%) _ Y (А) Ji UL) 


À, BX, (Aulus) - aA, Y (Al) Bd (4,11) -AAJ (4,1) 
» LD) L. Сау 
=o f = (8:627 eur (9:6 ао (ау 
Pour le respect de la condition aux limites en z=0, il vient : 
Terme axiale x œ„A,„Cosh(2,z)- B,Sinh(A, z) 
La solution se développe donc sous la forme d'une série : 


n=l,+œ 55 (A14) x QA, J, (A14) BY, (Anla) EN 0,4,1, (4, l 


nrl 


T(r,z)= Y. B,(a,4,Cosh(2,2) - B,Sinh(A,z)) Jour) Yr) ] 
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La dernière condition en 2=1, donne : 
/,@)= УВ, (oA, (oA, Sinh(A,L) — B,Cosh(A,1.))+ B, (oA, Cosh(A,L) — B,Sinh(A,L Up. 


n=1,+00 
Ce qui donne pour les coefficients du développement en série, la solution suivante : 
C, = (0,4, (aA, Sinh(2,1)- B;Cosh(2,1))+ B(0,2,Cosh(2,1)- B;Sinh(2,1.))) 
C, =2,Cosh(2,L)e,8,-a,p,)+Sinh(2,)a,a,2, — B.B.) 


Lz 
Lä, o | JA) Yr) ) 
B _ l4 BJ, OA) 7A, J OU) BY, (L4) = oA, Y (Al) 
: C, |o, c) 


1,5 15, 
=> B„ = аЛ) В, = (dr r FONAN 
l4 l4 


В _ By, 
Bd (Anla) E QA, J, (Anla) Pi (CAM) - a Y, (CAM) (0,4, Соѕ/(А, 2) D B,Sinh(A,z))x 
T(r,2)= У, C, |o, Cr)| 
Te ( J (A,r) D Y Ar) ) 
BJ, (4,11) = KE (4,1,1) BY, (Ala) St QA, Y, (Ala) 


BJ (А ) = 0,4,7, (4,1,,) A BY, (2,1,42) = a,^„Y, (AL; ) 
BJ (A, ) ж QA, J, (4,1,1) BY (Anla) um 0 A, Y, (Anla) 
Cette expression même sous cette forme la plus indigeste possible permet d'envisager tous les cas 


limites possibles simulant les conditions de Dirichlet ou de Neumann, homogènes ou 
inhomogènes : 


avec À, valeurs propres solution de 


a, =0,1 
p, =1,0 
a, =0,1 
B, =1,0 
a; =0,1 
В, =1,0 
a, =0,1 


B, = 1,0 
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Exemple: Cylindre plein (r,z) soumis à des conditions aux limites mixtes de Robin homogènes en z 


Pour le problème général : 
OT(r,z) 1ôT(r,z) O'T(r,z) 
+ + = 
др? r Or дг? 
T(r,2) fini 
aT, (r,z)+ BT(r,2) = f.) 


оТ, (rz) + B,T(r,z) | =0 


0 


eT (r,z)* B,T(r,z)| _ =0 

Pour assurer la premiére condition aux limites axiale homogéne de Robin, il vient : 

Fonctions propres axiales > Ф (2) = 04A, Cos(A,z) — B,Sin(A,z) 

Pour assurer la deuxiéme condition aux limites axiale homogéne de Robin, il vient : 

Ф (2) = 04A, Cos(A,z) — B,Sin(A, z) 

Ф„'() = -A, (oA, Sin(A,z) + B,Cos(4,z)) 

oO, (.) + B®, (l) =-a,A„(a;4,Sin(A„1.) + В,Соз(д,1.))+ B(a;2,Cos(2,1.)- B;Sin(2,1.))= 0 
À, valeurs propres solution de Sin(A,L Zoe, + BB.) = A„Cos(2,„1.)(B,a,- Bo, 

Dans le cas des fonctions sinusoidales, les normes se calculent ainsi : 


|o, c = f ao, cy = Joc EH cr -+0,00', el 


n n 


1 


Ce qui donne dans notre cas : 


z,=0 z,-l Ф, (2) =о,А,Соѕ(А,2) – B,Sin(A,z) – Tc = Q5 À, Sin(A,z) + B4Cos(A,z) 


1 


|o, c = ioun j SE e») P (2)9", e 


0 


Ф (0) = а,4, -———- 


D, 0.) = a 2, Co) PSG) -EC aA Sin.) Сох.) 


n 


Га л,Соз(2,.)– Ват, )) +(a;A,Sin(A,1.) + B,Cos(A,1.)) ) 
Эа = (a,A,Cos(2,1)- B,Sin(2,L)\a,2,Sin(A,L)+B,Cos(A,L))-a,B, 


n 


le, of = 


La dernière condition en г=1, donne : 
f.) = > `B,L (Q, LY(054,Cos(A,z) - B,Sin(A,z)) 


n=1,+00 


Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques - p 140 


Ce qui donne pour les coefficients du développement en série, la solution suivante : 


Í dz f (z) (a@3A,Cos(24,2)-B,Sin(2,z)) 
B = 0 


n 


1,(,1,)|e, (2) 


L ГА 
C, = | а f,(z) Cos(A, z) S, = | dz f.(e) Sin(A, z) 
0 0 
_ 1 
"KALEO 
Et la solution s'écrit 
1 ((@,4,Cos(2,1,)- B,Sin(A,1.)) +(@,2,Sin(a,1.)+ B,Cos(A,1.) ) 
| 


(a,2,C, - 8,5,„) 


— 


D, (GI = 


HEEL 


n 


À, valeurs propres solution de Sin(A,L YA, a.a, t 8,8, )= À,Cos(2,l (8,0, — Ж 


C, = [dz f. Cos, 2) S, = [dz 7,02) Sin(A, z) 


TO ye. (a;A„C, – B3S„) 10, r) 


a,À,Cos(À,z)- B,Sin(À 
Ио 
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Exemple:Cylindre creux (r,z) soumis à des conditions aux limites mixtes de Robin homogènes en z 


Pour le problème général : 
OT(r,z) lOT(roz) O'T(r,z) 
+ + = 
or? r Or дг? 
aT, (r,z)+ B,T(r, z) NE 


oT, rz) B,TG.z)| = 7.0) 


oT, (rz) + B,T(r,z) | -0 


0 


оТ (rz) B,T(r, 28 -0 


Pour assurer la premiére condition aux limites axiale homogéne de Robin, il vient : 
Fonctions propres axiales > Ф (2) = 04A, Cos(A, z) — B,Sin(A,z) 


Pour assurer la deuxiéme condition aux limites axiale homogéne de Robin, il vient : 

Ф (2) = 04A, Cos(A,z) — B,Sin(A,z) 

Ф„'() = -A, (eA, Sin(A,z) + B,Cos(4,z)) 

oO, (.)- Bb, (l) =-a,A4,(a;4„Sin(2,1.) + B.Cos 4,1.) В, (nA, Cos (4,1) - B;Sin(2,1.))=0 


nz 


À, valeurs propres solution de  Sin(A,l Ya, Zoe, + B,B,]- À,Cos(A,L Ba, - Bo, 


nz nz 


Dans le cas des fonctions sinusoidales, les normes se calculent ainsi : 


lo, c = f” dz e, ey = Joc EH e) - 40,00", el 


n n 


1 


Ce qui donne dans notre cas : 


z,=0 z,-l Ф, (2) =о,А,Соѕ(А,2) – BSin(A,z) — m = Q5 À, Sin(A,z) + B4Cos(A,z) 


n 


L 
2 1 Lbs 1 I i 
|o, C = С ta o) "bc, ol 
o,'(0 
P,(0) = «2, — р, 
| © (L) ; 
D, (L) = 0,4, Cos(A,L) - B,Sin(A,L) — E E = 05A, Sin(A,l.) + В,Со5(А,1.) 


(le A. Cos (4,1.) - B,SinCA, LY + (аА, Sin(A,L) + B,Cos(A,1.)) ) 


1 
Ie, c =>) 1 | | 
2|+ T (0,4,С05(4,1.) ~ B,Sin(2,L)\a;2,Sin(A,L)+B,Cos(A,L))-a,B, 


п 
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La première condition en r-l,, donne une partie radiale sous la forme : 
Fonctions radiales 


R(r)= AI(Qur) * B K (A,r) 


= -AL(,r)- B K (A,r) 


a, (4 L (2,14) - B KA, ))+ B, (4 I(A,14) +B Ky (A, )) =0 

Ala, T, ula) + BU = Bla KU) BEA 

Fonctions radiales R(r)- (e,K,(A,1,) BK (4,1 D) A,r) +(0,1,(2,1,4)+ BIG ))K ur) 
La deuxième condition en r-l,; donne : 

A, = (eK Qi) - Ks (A1) 

B, = (1,14) + BAAL 

R(r)- A„I (A,r) - BK, (Ar) 

R'(r)= A (4,1, (A,r) - BK, (A,r)) 

a Rr) В,К(т) _ SA A OUS) B.K AA le BA, I O12) B Ky ula) 
Posons D. = G5 A (A, 1,1.) - B, K s) BA, 1,1.) * B, Kol.) 


f.(@)= Y B,D, (aA, Cos(A,z) - B,Sin(A,z)) 


n=1,+00 


Ce qui donne pour les coefficients du développement en série, la solution suivante : 


Í dz f. (z) (e,À,Cos(A,z)— B,Sin(A, z)) 
B = 0 


n 


rn 


Dt 


C, = [а f.) Cos(â, г) S, = dz f.(2) Sin(A, z) 


1 
B = pm. -BS 
n D, sc) 3 "nn P, „) 
Et la solution s'écrit 


(le, A Coat A) — B,SinA,L)Y +(a;A,Sin(2,1.)+ B;Cos(A,1.)P ) 
| 


dE =] 1 
2062 eA Co) Sin. Uer Sin.) B;Cos(A,1.))- ap. 


n 


À, valeurs propres solution de Sin(,L YA, Zoe, t 8,8, )= А„Сов(А.1, X8,o., Е Bol 
1. L 

C, = | dz f.(2) Cos(A, z) S, = | dz f.(2) Sin(A, z) 
0 0 


4, T (aK, (4,11) CS BK; Al) 
B, = (al, (4,1,1) + ВІ, (A,1,)) 

D,, — G-A, (4,4, (4,1,5) я B,K, (4,1, ))+ B, (А „1, (GJ. + B, K, (,1,,)) 
T(r, z) = У (0,А,С, - B.S,) (4,1, (4,0) RE B.A, (4,r)) 


A7. ff D. (2,4, Cos(4,z) - B,Sin(4,z)) 
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Exemple : Soit en coordonnées cylindrique (r,z), la solution de l'équation de Laplace sur un cylindre 
plein infiniment long, de rayon l, partant de l'origine z=0 (demi axe) avec des conditions aux limites 
de Dirichlet. 


2 2 

д AUZ m 1 OT(r,z) 4 Ô Aue =0 T(r,z) 
ôr r r Oz 

T(r, 2), = 0 


T(r,z), , =f,(r) 


„o fini T(r,z), _„ fini 


La condition homogéne en r fixe le choix des fonctions propres de Bessel en r. La solution en z est 
une combinaison linéaire des fonctions exponentielles. La fonction de Bessel de deuxiéme ordre 
n'est pas admissible car la solution en r=0 doit être finie : 


ET) S Geh eC 


Ar Ar 
n " j 


R,(0) est borné => D, =0 
№,(1,)=0= J (2,)=0 


valeurs propres À, solution de J,(¿,)= 0 


2,2 
T(r,z) , fini— 2,(2) = Be" 
Anz 
T(r,z)= Y Be* AC 
À, 40,00 r 


On rappelle la formule pour le calcul de la norme des fonctions propres : 


kË p. (ry +D, (ry ) 


Ar 
; - n-0 nl, d) JC 


r 


Ip cy = [ dr r b (ry = 


o, (,)-0 9) - 740.) 


1, 
A La (r)? + Ф (r)? play (1)? | 1? ТА рд) 
A n n Я r 2 n Nr r 2 1? 1 п J? (À ) 


r 


D "s = = = T 
= | Bisi 2 2 2 2 
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La solution se présente sous la Torme d'une série de fonctions cylindriques : 


Ј (A, 
A, tq 700,)=0 |Ф„()] = ^0) 2 
8, - [d Bj, 
"but 
Se 
Bë cy 5”) 
T(r,z) -2 —— A, tq J (à,)=0 
(r Z) z. TOA q al ) 
si f.(r) =1, pere 
1, l 24, 2 1° 
в, =] dr r JC ES 544] [arr 407 (z) PAGE = 209 
_À,z 
TE Ағу 


= T(r,z) -21/ y - 


=== — À і Y. P, =0 
n#0,+œ AJ (A) n tq oC 4) 


Si la solution se décline à l'aide des fonctions propres : 
valeurs propres À, solution de | J4(A,1,) = 0 


T(r,z) = Y Ве Lr) 


An 40,00 


D (r) = Jy (Ar) 


Alors la solution s'écrit 


À 
2, tq Jo(A,l.)=0 |o, DG) 
B, -Í drr f.) 2,0,0) B, ы P= 

° KAQ) 


2 B, e" J G.) 
T(r,2)= Z D à, tq 44(4,)20 
L z, JUL) Í x 


si f.(r) = T, r=4r>r=— 
ү 1 ү E Ï 
б! ' ' va Hal, * 
в, пама. DEELEN Ae" = To ZA AL) 


n 


27, S е (Аһ) 
>T V „Z = 0 0 n À [1 J À =0 
CD 2. AAOD C490 


r 
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Exemple : Soit en coordonnées cylindrique (r,z), la solution de l'équation de Laplace sur un cylindre 
plein infiniment long, de rayon l, partant de l'origine z=0 (demi axe) avec des conditions aux limites 
de Dirichlet. 


2 2 
д T(r,z) + 1 OT(r,z) + д T(r,2) = (0) T(r,z). , fini T(r,z2)|._, Jini 


or? r дг Oz? 
T(r,z r=l, = Г Т(у, 2 r=l, =0 
( Js 9 ( la 
T(r, z) z=0 = T, 


Le problème s'inspire en partie du problème précédent et représente la superposition des deux 
problèmes suivants : 
О =fr.,z/r e [0,1,],ze[0,»]]=Q,UQ0, О, =fr,z/r e[0,1,],ze[0,4,1) О, ={r,z/r e [0,1,],z [1,00] 


То, (7,2) 


reo] = Т, То, (7,2) rel04,1 = To, (22) rel0,4,] 


To, (r, z) 


r-l, = 0 


2є[1, ,+о] 


reto] a To, (r,2) re[0,1,1 To, (r. z) 


La solution du deuxième problème est obtenu à partir de l'exemple précédent : 
2 C) 
А, ee Т, (rz) - > B, =Z m > =*> 
| L, n#0,+00 Ji (¿,L.) 
Bn étant des coefficients déterminés d'après les conditions aux limites de surface en z=lz. 


Pour ce qui est de la solution du premier domaine, on peut facilement la ramener à un problème 
homogène en considérant la solution triviale du demi-cylindre porté uniformément à la 
température TO, et celle d'un cylindre de hauteur z aux conditions homogène de Dirichlet sur la 
surface latérale et la section basse et une fonction limite donnée sur la face supérieure. Il vient 
donc pour le premier problème : 
EE 

| lee Ji (A,L.) Sinh(A,£) 
An étant également des coefficients déterminés d'après les conditions aux limites de surface en 
z=lz. 


On peut remplacer le terme constant par la partition de l'unité obtenue avec un problème aux 
limites normalisé à une valeur de face supérieure égal à 1, et toutes autres faces portées à O, soit : 
2 J (A 2 J, (A,r) LJ (¿I 
-2 у AUDI E у II A) лал а,)=0 
L, AJ AL) ее (1) А 


nr nr 


n 


Ce qui donne la solution du premier probléme : 


2 y Jr) ЛО) ү, _ 4 Sinh, 2) 
pou TAIL `À Sinh(À, 1.) 


n 


Ta (r,z)= J À, tq J (A,„1.)=0 
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Pour déterminer la solution il suffit maintenant d'assurer les conditions de continuité С, et C, sur la 
section supérieure du tube cylindrique : 


LJ (AL) 
о: 


п 


Co) __,ь — д COSA e — LG: _ al а) 


r Sinh(AL) КЕ” " Sinh(A,L) `" А, ° Sinh(A,L) 
= Á, = те Sinh(A,L.) 
= B, = A) тес Cosh(A,L.) 


n 


Ce qui donne pour le problème : 
eT(r,z) n 1 ôT(r,z) = eT(r,z) 


=0 T(r,z ini T(r,z ini 
or? r дг д2? Gas À ( se fi 
То) SCH s Seen E 

T(r,z), , = T 


la solution suivante : 
А, tq J4(A,1,) = 0 


SE 


2 Jor) rr 
Io (r,.z)=—T Cr N-e ^ Sinh(A z 
d 1, pa 21205 i ) 


2 Jo (Ar) E 
T. (r,z)=—T "O Cosh(A,L )e ^ 
o, (^, z) 7 »» AJ. (4,1) 


r nr 
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Exemple : Soit en coordonnées cylindrique (rz), la solution de l'équation de Laplace sur un demi 
anneau cylindrique infiniment long, plein de rayon intérieur l,, et extérieur 1, partant de l'origine 
z=0 avec des conditions aux limites homogènes de Dirichlet sur les rayons intérieur et extérieur. 
eT(r,z) n 1 OT(r,z) E O^T(r,z) _ 

or? r Or д2? 
T(r,z) , = 


0 


T(rz) , = 

T(r,z), „= f. (r) 

T(r,z) , fini 

En tenant compte de toutes les conditions aux limites, la solution prend la forme : 


70,11) = Y (A14) et D (r) = J (À,,r) Y, (4„r) 
Jo (4,15) Y,G,1.,) " Jo (4,15) Y (Anl) 


1,2 2 2 
"m =| LL a viel - 1 (Hd) л) наа) AG 
" 2 A 2 Y (A1) Jo (GJ. 2 Y) Jo (4,1,2) 


À, tel que 


l4 


12 12 Vë 
B, = [ar f) О) Br = fr fA). [ar Ф, Z = Lt 
l4 l4 la 


Jo (4,15) Ү,(4,/,) 
| В By, | e^ { Ji (A,r) Y (A,r) | 
Jo (LA) Y, (41.2) Jo (4,1) Ү,(4,1,) 


2 2 
2 [ { J, (Al) AURA) | l { J, (4,1,1) Y (l, | | 
r2 rl 


T(r,z)= 2 


Jo (4,15) Y (A1) Jo (4,1,2) Y, (Anla) 


Si la fonction limite f,(r)=T, alors : 


fra f. (r) D (r) = T, La Jl )- TL, Ji (4,1,1) 1, Y (A5) - l4 Y (Anla) 
la | ; À, Jo (4,1,5) Y, (4,1, ) 


А, Jo (AL) Y L3) Jo LJ. Y (A,L5) 


o.) =! JS) О) | | ЛА) YO) ||, 
À Pr) 4013) VIU UT) 


x 1 J, (4,15) Y (AL) + l J, (4,1,1) Y (Ala) 
PUA) KA) JU YO) 


nra 


|, гә - 0429), (342. Ha) 


2 


eaa n) 
T(r,z) = 21, > 7,0,1) 04,2) 


1 
NEE iM 
С Jo (GJ. Y (AL) s Jo (4,15) Y (A1) 
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Exemple : Soit en coordonnées cylindriques (rz), la solution de l'équation de Laplace sur un demi 
anneau cylindrique infiniment long, plein de rayon intérieur l,, et extérieur 1, partant de l'origine 
z=0 avec des conditions aux limites homogènes de Dirichlet sur le rayon extérieur et de Neumann 
sur le rayon intérieur 


OT(r,z) lOT(roz) O'T(r,z) 
—————À + р 
or? r Ог д2? 

T(r,z)) =0 


0 


T(r, 2), = 0 
T(r,z). , = 5.0) 


E: < P dE La solution s'écrit : 
(А) = Ү,(4,,) D (r) = Jo (A,r) HAr) 
| J, (4,1, ) Y, (Al) É Jo (4,1,,) Y, LJ. 


12 2 2 
Ф О = p^ Ф', (° + (r)? = LÉ Ү(А„1„) J, (LJ. 1 Ji (4,1,1) Y (AM 
I 2 2 Y (41,52) Jo (GJ. 2 Jo (4,1,,) Y, (4,1,52) 


B m By, 
Jo (LJ. Ү,(4,,,) 


B„=[drr f,O) Jr) В, = |а уа) [drr Ф, 0) = 


| B, B, je Jur) aan) 
AGREE) Yo) J, GL.) FA) At Jo.) _ D,a) 
2 2 n 
dee [ (ie gel dE Sall Aula) EU) 
r2 rl 


Y (A15) Jo GJ. Jo LJ. Y (Al) 


Si la fonction limite f.(r)= To alors : 


[^ f (r) D (r) = To Li JU) = L, Ji (A11) _ La HOA, Y (A1) 
La ( í À, Jo (4,1) Y (L4) 


m T, l J, EM Y(A,l,:) 

ü À, x Jo (А) Y, (4„1,2) 
JA) YO) | A) А) |, 
Jo (4,1,) Ү,(4,,,) Jo LJ) JA) 


1 
CR 2 2 
es l À XL) J, dala) | l 1 J (Al...) d 
' Ү,(4,1,,) Jo (4,1,,) ' Jo LJ. Ү,(4,1,,) 


T(r,z)= 21l. 


Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques - p 149 


Exemple : Soit en coordonnées cylindriques (rz), la solution de l'équation de Laplace sur un demi 
anneau cylindrique infiniment long, plein de rayon intérieur l,, et extérieur 1, partant de l'origine 
z=0 avec des conditions aux limites homogènes de Neumann sur le rayon extérieur et de Dirichlet 


sur le rayon intérieur 


2 2 
CT(r,2) | 10T(2) d TEZ 4 
др? r Qr дг? 
T(r, 2) = 
Т, (Р, 2), 2 а 
Т(г,2) „= 5.0) 
fini 

А, tq Ја, LJ Y (À,1.,) Ф (r) = J, (2.4 r) HA, r) 

Jo (А, Li) Y, 0(4, la) Jo (А, La) Y, 0(4,/ л) 


J (Al a) Y (ALa) 


J, (4,1 rl ) 


1,2 
1| Ф' (° 1 
PRU! ear EL 


„= fur eu n, = [ote f,O KAN fersa, (r)= 


B 


Jn 


B J (A,r) 


Yn 


Y, (1„r) 


Y, (А. oi 


| g TI) 


i ( al LA) 


| 


E 


В 


J (2, la) Y, o (Ål a) 


B 


Yn 


T(r,z)= 2 


"i la fonction limite f.(r)= To alors : 


m0, oo 1 2 J, (A, 1) Y, GI AT l 2 J (Al Ze TY GET 
FACE) XOU AGE) YOL) 


| rdr f (r) D (r)= 2 ал (Al) 


1, Y (A „l ua) Li Y (4, = 
Y, (A, la) 


J Ai Ze Y GI KS 


SE l T, J (A, Al Y (4, La) 
CORSI S 

[ш a) IA D JA) LA) ) 
T(r,z) =2T,l„ V. (Аһ) Q, ы) ЈА) KA) 


J (4,1 м) 


AM 2 J (Al „2 ) YA ms Lp? 1 Y (A, Lo 
LIAL) LOIS) LOGE) QU) 


| 


| 
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Exemple : Soit en coordonnées cylindriques (rz), la solution de l'équation de Laplace sur un demi 

anneau cylindrique infiniment long, plein de rayon intérieur l,, et extérieur 1, partant de l'origine 

z=0 avec des conditions aux limites homogènes de Neumann sur le rayon extérieur et sur le rayon 

intérieur 

e^T(r,z) , Lëtz), O°T(r,z) 
др? г Cr Oz? 

Tez) = 


-0 T(rz) , fini 


Tz) = 
TÉLÉC 


Le développement en série de la solution donne les calculs suivants : 
Ф (r) SS J (A, r) Y, ACA r) et А, t J (4l 12) _ Y Q, 1,5) 
J, (4,1 л) Y (A, 1) J, (4,1 л) Y (4, la) 


eo =Z DES 97) Б ^ (2622 DE AU Lala E (29 L) X, satay) 
2 À, l 2 J (À, Li) Y (A, La) JU м) Y (A, La) 
12 В 


B, = (rdr f) Jr). В, = [rar sona, r) ferr tivo, (r de 5 E ^ 


sch B,, By, | AQ) X, à 
E Ala) МО) А) 0,1) 


“э, "É AS) X4. ш) j E nce Si 
rl 

J, (А, La) Y (А a) J, (4,1 м) Y (А Al 
Mais Il faut rajouter à la solution précédente le terme de valeur propre nulle. Du respect des 
conditions aux limites homogénes radiales de Neumann, il s'agit d'une fonction radiale de la forme 
A+B Log(r) qui ne peut être qu'une constante. Soit la solution finale suivante : 
a (r) A DAND aa q AO YA) 
J.Q, by Y (A, la) J DER Li) Y Q, LA 


1,2 


B, = [rar Go 440,0. By ан в, [rr fO 


Pu B, A J LAN) | Y, 2| 
Tio 2A 42 IAL) KUIDAI A) Kala) 


+2), 
e =li n#0,+00 [ E lj) GER l E La) Y, (2, si d 
r2 rl 


J.Q, Li) Y Q, La) J.Q, Li) Y Q, la) 
Si la fonction limite f.(r)= To alors : 


| rar f.) Ф, (= ЕЕ A Ou) ia Gua LY. a» : 


la 


lochy 


в,= [rdc ROET, 7> =T(r.z)=T, 
l4 
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Exemple: Disque creux (r,O) soumis à des conditions aux limites mixtes de Dirichlet et Neumann 
homogénes en r 


On rappelle que les solutions de l'équation séparée en (r,O) respecte la périodicité en Ө 


T(r,0)=e+ f In(r)+ V en +dr” \a,Cos(n0) + b, Sin(n0)) 


n-l,oo 


On fixe les conditions aux limites de maniére à ce qu'elles soient paires selon l'axe des abscisses 
cartésien, cela induit que la solution suit cette propriété également. La combinaison des deux 
conditions aux limites extérieure et intérieure, et leur symétrie sur l'axe des abscisses implique que 
la solution du probléme posséde elle-méme cette symétrie qui se traduit par la propriété suivante : 


T(r.0) = T(r -0). Pour que cette propriété soit effective il suffit que la constante bn soit nulle. Les 
solution s'écrivent donc 


T(r,0) =е+ f In(r)+ Y (en +dr” Jcos(n6) 


n=1,00 


Problème n°1 

Sans restreindre le problème on peut fixer le rayon intérieur du disque à k et le rayon extérieur à 1, 
l'angle 0 variant de O à 2 x. Dans un deuxième temps on fixera les conventions suivantes rayon 
intérieur du disque à 1 et le rayon extérieur à k 


10 ÔT(r,0) 1 OT(r,0) OT(r,0) 10T(r,0) 1 OT(r,0) 

d e Ces 27 WIS 2 `= 
r òr Or r 00 Or r Or r 00 
T.(r,0)-hT(r,0) _ =0 


| Cos(0) 0e[- 5,7, 
T.(r,0)+RT(r,0) _ = + 


Cos(0) Oe E 


л 3x 
0 0e[—,— 
SE 54 


La condition aux limites en r-k, donne pour toutes les valeurs de n la condition suivante : 


posons С,(0) = 


n(c, e" а") сна") 202 а," +h) = с", п) 
> d =c k?" (n- hk) 

" “ (n+hk) 

et pour la valeur propre n=0 cela conduit à : 


X I d _ = ES e 
п he— f h In(k) = 05 f (1- hk In(k)) = hke = f LETTO) 
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La solution s'écrit donc : 


т-а s sn рн UE AE "e s(n0) 


i) 
(1- — (k)) my (n+ hk) ` 


Si maintenant on applique la condition aux limites en r=1, cela donne 


а ШИ а + д2" " h, ) » (n— k) 
(E Go WEE 1-4 "HD losen = C,(0) – 5e — h, Ye Ge LD ce 


n=1,00 
L'orthogonalité des fonctions propres du problème séparée donne les relations suivantes sur [0,л] 


Б 
Í Cos(n0)Cos(m0) = =J = > si n=m + 0 
7 si n =m =0 


En intégrant sur [0,л] l'équation multipliée раг la fonction propre Cos(n8 ) (en changeant dans la 


série n par m) cela devient : 
п= 0 


hk л/2 
eh -= f Cos(0)40 – her 
(I—hkln(k)) 5 


_ (1— Ak In(&)) 
(hk  h, (1 — hk In(k))) 


= 
Zefi k? ne). cor tno ie [veio CRD Ï Cos?(0)d0 = | Сои ET 
á (Ab) À 2 (+0) | Je 4 
1 
G 
Al E= e +h, 1+ p e ш) 
(I+ hl) (I+ hl) 
n>l 


л/2 
Len 1-k” GRO [Cos(6)Cos(n6)40 — 2 с, xt (=) 
2 (n+hk) A 2 (n+ hk) 


avec Í Cos(0)Cos(n0)d0 = cl "Z ) /(t- ny 


M 2Cos(17/) 


n ka KL, dru =) 
(п+ hk) (n4 hk) 
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La solution complète est donc 


T(.8)- (1— AK In(K)) j. Ak 
`^ Gm(hkeh,ü-hkln(k))A ` (1—AKIn()) 


h + 


1 (1- AK) y 


@ 5 
= 2) (I+ hk) 
(1+ hk) (I+ AK) 


am + 


Cos(nz E 
2 > oss) @ » hD 2 
uice [> caa), h, c e en] (n-- hk) 
(oi Ak) (oi hk) 
1 r 
TO = en en OR LUE 
i Ë Ë ch a 
d 0-40 Jen [rei mJ (1+ А) 
TEENS (I+ hk) 
Cos(nz А 
2 > oss) @ pe DA 2 
He DR e eA). ie ph ean] (п+ hk) 
GERE) (ҮЛ) 


Problème n°2 
Passons au second cas où l'on fixe le rayon intérieur du disque à 1 et le rayon extérieur à k, l'angle 
variant de O à 2 n. 


10 ôT(r,Q) 1 O0T(r,0) OT(r,0) 1ÔT(r,0) 1 O’T(r,0) 

r + 2 2 = 2 + + 2 2 = 
r òr Or r 00 or r © r 00 
T.(r,0)-hT(r,0) _ =0 


0 


| Cos(0) Get. 
T'(r.0)+hT(r,0) = 2 2 


Cos(0) Get. 
posons С,(0) = 2 2 


л Зл 
0 0e[—,— 
SE e 


On rappelle que La solution est de la forme : 


T(r,0)=e+ f In(r) + Y ler" +d r” XCos(n0) 


n=1,00 


Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques - p 154 


La condition aux limites en r-1, donne pour toutes les valeurs de n la condition suivante : 


a snos sd de PER) 
(n+h) 


n(c, - 


n 


et pour la valeur propre n=0 cela conduit à f-hez0-f- he La solution s'écrit donc : 


T(r, 0) = e(1+ h, In(r))+ > c |+ + о r Jostn0) 


n-l,o 


Si maintenant on applique la condition aux limites en r-k, cela donne 


А Уе, de = CA) p Jcostn0) eh +h, Ye = = EE C,(0) 


L = des fonctions propres du problème séparée donne les relations suivantes sur [0,л] 
O0sinzm 


Í Costnó)Costme) d0 = z5 n=m #0 


T sin=m =Ò 


n=l,œ n=1,00 
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En intégrant sur [0,л] l'équation multipliée раг la fonction propre Соѕ(пӨ ) (en changeant dans la 
série n par m) cela devient : 


п= 0 
л/2 л/2 
ex = [Соѕ(0)аө — ел [Соѕ(0)аө =] 
0 0 
D LM 
z(h + h,k) 
n-l 
x/2 л/2 z/2 
T el WA = [Cos*(0)d0 =" h c, p ЕА) f Соз°(Ө)аө = | CRC 
2 k(n+h)) 3 2 k(n+h)) 5 d 2 4 
E EN А Ы 
Dos k'(n+h)) ^4  k(n+h))) 4 
"m 1 
z 
Al MEA) en 1 97,242 
(n+h) : (n+h) 
n>l 
zí2 
Zen prds T) ы. [ Cos(0)Cos(n0)d0 – 2 с, p ХМ) үз 
2 (n+h) d 2 (n À) 


avec Ï Cos(0)Cos(n0)d0 = Cos "7 ) /(t- п) 


Zela pri D'A el ‚| p DA E I о") fn 
2 (n А) (nh) 2 

2CosV TX 
n оз(плу) 


| a =h) a + vi (п) i || 
(n+h) (nh) 


Problèmes aux limites de Laplace en coordonnées cartésiennes, polaires et axi-cylindriques - p 156 


La solution complète est donc 


k 


Є z(h + h,k) 


(1+ bh In(r))+ 
I Ë ад 
dÉ (n—h) kal Mel (n+ |) 
(n+h) (n+h) 
Cos 7. _ 
2 osl A , GO h) a 


> x= 
л, Lech, (n-h) e. h, efis (n-h) 20) (nh) 


(n №) (nh) 


210 + 


an 


Exemple : Disque creux (r,0) soumis à des conditions aux limites de Dirichlet en r 


Nous allons illustrer un problème avec des conditions aux limites présentant une symétrie par 
rapport à l'axe cartésien des ordonnées (symétrie par la transformation Ө-> л - Ө). 


On rappelle que les solutions de l'équation séparée en (r,O) respecte la périodicité 2rt en Ө 


T(r,0)=e+ f In(r)+ D +dr” Ya,Cos(n0) 4 b, Sin(n0)) 


n=1,00 


Si la condition aux limites présente une symétrie cette transformation : Ө-> л - Ө cela induit que la 
solution suit cette symétrie également. En effet, la combinaison des deux conditions aux limites 
extérieure et intérieure, et leur symétrie sur l'axe des ordonnées implique que la solution du 
problème possède elle-même cette symétrie qui se traduit par la propriété suivante : 


T(r,0)=T(r,7 - 9. Pour que cette propriété soit effective, reportons là dans l'expression de la 
fonction: 

a, Cos(n0) +b Sin(n0) = a, Cos(n(z —0)) & b, Sin(n(z —0)) 

= a,Cos(n0)  b,Sin(n0) = (-1)' (a,Cos(n0) - b, Sin(n0)) 

—a,-(-M'a, b, --(-1'5, 

=> a,, # 0 et Б», # 0 


T(r,0)=e+ f In(r)+ V (c. r?" +d,,r °" JCos(2n6) + > C +d, т" )Sin((2n +1)0) 


n-l,oo n=0,0 
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Sans restreindre le problème on peut fixer le rayon intérieur du disque à 1 et le rayon extérieur à k, 
l'angle Ө variant de 0 à 2 л. 


180 ÔT(r,0) 1 Ô°T(r,0) OT(r,0) l0T(r0) 1 Ô°T(r,0) 

d F = ITERS 7. ГА dg z =0 
r Or Or r 00 Or r Or r< 00 
T(r,0)| -0 


ms 0 € [0, л] 
10 8 e[x2z] 


T(r,0) 


r-k 
1 0e€[0,z] 
0 elei) 


La condition aux limites en r=1, donne pour la valeur n=0 : € = 0 et pour toutes les valeurs de n : 


posons С,(0) = | 


C, +d, =c, +d, .. =0 . GEN 
2n 2n 2n+1 2nd . La solution s'écrit donc : 


T(r,0)= f In(r) + py D =" Jcos(2n0) + b» = š )Sin((2n + 1)9) 


n-l,o n=0,00 


Si maintenant on applique la condition aux limites en r=k, cela donne 


f In(k) + Y c, (k?" =k 2"JCos(2n0)+ Se, (k?"! =k ?"3)Sin((2n+1)0)=C,(0) 


n-l,oo n=0,0 
L'orthogonalité des fonctions propres du problème séparé donne les relations suivantes sur [0,2л] 
2л 
[ Cos(n9)Cos(n'0) dO =0 si n z n' 
0 
2л 
| Cos(n0)Sin(n' 0) dO =0 si n z n' 
0 


2л 
[ Соз(2лӨ)Соз(2л'®) 40 =0 si n z n' 
0 


2л 
| Cos(2n0)Sin((2n'+1)0) d0 =0 si n z n' 
0 


2л 


In docs 


0 


2л 

Í Cos? (2n0) аб = 
V 

f Cos? (ne) 40 =m 
kd 2л 

[sin (2n +1)0)40 = 2x — f Cos? (ол +1)0)40 = z 
0 0 


2л 
Í Sin*(no) dO = л 
0 
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En multipliant par l'une des fonctions propres et en intégrant sur [0,2rt], il vient : 


п= 0 
1 
21n(k) 


2 f In(k)z = je 040 = [40-55 f = 
n>0 


2л л 
x c, (k? =k" )= Í C (0)Cos(2n0)d@ = | Cos(2n0)d0 = 0 > c,, =0 
0 0 


m 6,4 7 Ek fa (9)Sin((2n + 1)0)40 = j Sin(2n +1)0)d6 


_ | Cos((2n + 2 S (1 = Cos((2n + 1)л)) (1 Е Cos(x)) 2 


(2n+1) Е 
Рета 2 
2n+ — л(2п +1 pn uper) 


2n+1 2n+1  2n+1 


0 


La solution s'écrit donc : 


lr) 2 Lo pt Jsin((2n + 16) 
T.(r,0) = 
Bis ) 21n(k) iE T 2, (2n + (enn E 


Pour le probléme symétrique à l'axe cartésien des abscisses: 
1 ə г OT (r,0) " 1 e?T(r,0) _ OT(r,0) " 1 OT(r,0) Í 1 eT(r,0) _ 


= =0 
r òr Or r^ op др? r дг r^ og 
T(r,9)|.., =0 
0 Өє[0,л] 
T(r,0)|_, = 
= 1 Өє[т,2л] 
0 Oeo, 
posons С,(0) = E102] 
1 Өє[л,2л] 
La solution se calcule pareillement : 
T,(r,0)= Inr) 2 28 _ г iz 00) 
2In(k) e Zei, (2n+1)(k?-k2") 
Et pour le problème : 
10 OT(r0) 1 0°T(r,0) OT(r,0) 1 ôT(r,Q) 1 &T(r,0) 
r F= = >= = =+ EG „ £20 
r òr Or r 00 Or r дг r 00 
T(r,0)| , EI 
T(r,0) _, =0 
La solution triviale de valeur propre nulle est 
T,(r,0)=1- In(r) 


Int? 
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Par le principe de superposition, la solution complète pour le problème : 
1 ô OT(r0), 1 Ə°T(r,0) _ O°T(r,0) ,l8T(.0), 1 Ə°T(r,0) _ 


= -0 
r òr or r^ og? or? r Or r^ og 
T(r,0)| , =Ó 
0, Өє[0, 
ks 2 e[0, 7z] 
r= |ô, Өє[л,2л] 


T(r,0)=0,7,(r,0)+0,T(r,0)+0,T,(r,0) 


T(r,0) = 0, 1— In(r) À 0, In(r) р 2 Län =r?" )Sin((2n + 1) 
In(k) 21n(k) Л „0 (2n + 1 2n+1 -k 2n 1) 


+6, | Inr) 2 |! ад 


эщ m, Оле =k?) 


(M) In) РИИ (r -r )Sin((2n +1)0) 
T(r,0) - 0, [ 2). 2In() (д, ш д,)+ = (6, 2» (2n "m 2п+1 d 


Pour un probléme plus général de la forme : 
10 OÔT(r,0) 1 oOT(r8) OT(r8) lOT(r0) 1 ô*T(r,0) 

r + = + + = 
r Or Or г? 020? or? r дг r^ 00? 
T(r,0), ,. -0 T(r,0), ,. = f,(0) 


Supposons dans une premier temps que la fonction fə(8) présente la symétrie par rapport à l'axe 
des ordonnées f9(0)=fo(n - Ө) alors la solution est à rechercher sous la forme du développement 
suivant (on notera X ou Y pour désigner une fonction paire/impaire par rapport à l'axe X ou Y, et + 
et - pour désigner une même fonction paire/impaire) : 

Fonction limite ` f, (m -0) = f,(0) 

Т; (7,0) =e + f In(r) + Y her” Td, r ?" cos(2n0) + Y (ca, r?" + d, r 2" )Sin((2n +1)0) 


n=l,œ n=0,00 


0 


n=0 conditionr=l, =e+fin({,)=0=e=-fIln(.,)=x f i 


rl 


condition r =l, > E = Í dO f,(0) = [de f, (0) Í dO f,(0) 


rl 


=> Cam Lg = /,C0)> f = 2] dO f,(6) 
= Casn°2 f,(0)=-f,(-0)= f =0 


partie valeur propre nulle => 
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Fonction limite f, (z -0)= f,(0) 
T; (r,9) - e* f In(r) + > (c. r?" +d, r” Jcos(2n0) + У (c, rm +d, ar )Sin((2n +1)0) 


n=1,0 n=0,00 


n>0 


2n -2n 
Se n -2n r r 
condition K= La => e + d, „l | = 0 => F d ü p | 


2nd -(2п+1) 
= Е И + T RE = 0 о Jl Ë p | 
rl rl 


condition r = l, => 


Z5 42) |. јао f, (0)Cos(2n0) 


ri 


Í, = [do f, (0)Cos(2n0) = fao f, (0)Cos(2n0)+ fao f, (0)Cos(2n0) = 0" =n -0 
1, = [49 f,(0)Cos(2n0)- | dO' f(x -0')Cos(2n(r -0") 
Toc fao f, (0)Cos(2n0) - [аө f, (9')Cos(2n0') = fao 7,00)Соѕ(2п0) + fao fot—-0)Cos(2n0) 


=> Cas n°1 f,(0)= f,(-0)= L = 2] dô f9(0)Cos(2n0) 


—Casm2 f,(0)=-f,(-0)= I, =0 


2п+1 -(2n«1) 2n 
(8) 42) |- [40 f,(0)Sin((2n+1)0) 


I 1 


rl rl 


I fao fa (0)Sin((2n +1)0)= fao f, (0)Sin((2n +1)0)- [do f, (x — 0)Sin((2n + 1)(z -0')) 


L. = fao f, (@)Sin((2n + 1)- [do f, (0)Sin((2n +1)0)= fao f (0)Sin((2n+1)0)- [дө f,(=0)Sin((2n+1)0) 
=> Саѕп°1 = f (0) = f,(-8)= 1, =0 


= Cas mi f,(08)--f,(-0)2 I, = 2f d0 f,(0)Sin((2n +1)0) 
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En conséquence la solution s'écrit pour les deux sous-cas du problème 


Jo(0)= f, (m -0) 7,00) = f$ -0) 
2[аө f,(0)Cos(2n0) 


f= [do LI) с, = — z — fao f, (0)Cos(2n8) 
0 I l 0 
«(87 5) 


2n -2n 
r E WEEN 
LS 2 L ) É ) | 
rl 
Toetab j= De, = Cos 
l 


Jo(0)= f,(0-0) f,(0)=-f,(-0) 


2 | dO /,(Ө)5їп((2п +1)9) 
Сы — p SES Cj 2nd -fao Jo (0)Sin((2n + 16) 
(ЕЧ Т 


2п+1 —(2п+1) 
Ka m is 
L, J Li 
Ty, x. (7,0) = => C; 2nd 


E Ë li fe | ET 
І 
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Et pour le problème général suivant symétrique sur l'axe radiale r : 

10 ÔT(r,0) 1 0T(r,0) OT(r,0) 10T(r0) 1 T(r,0) 
r + = + + 

r òr Or r^ op др? r P r^ ӨӨ? 

T -f$4(0) tq f,(0)= (7-90) 


=0 


avec l'hypothèse de symétrie par rapport à l'axe des ordonnées fe(0)-fo(rt - Ө), on retrouve 
également les deux sous-cas : 


Jo(0)= f (x -0) f,(0)= f,(-0) 


К 2[аө f,(@)Cos(2n 0) 
f =[d0 RO) о, == 
( 


2n -2n 
RO 
r2 1, 
2п -2n 
r SS NE Ed 
^w 2 ( Е p | 
Ty. xy, (r,0)= == = > C jn ER Cos 
tri =0, 


Jo(0)= f,(-0) f,(0)=-f,(-0) 


2 | dO f,(0)Sin((2n Se 
Ci, 2п+1 = | 40 7,00)51п(( (Qn + 10) 


2n+1 —(2п+1) 
l4 


2n+1 -(2n41) 
PF r 
i" e 


Ty, x- (r, 9) = E -Y С j 2nd 2nd (2n41) Sin((2n + 10) 
n=0,00 Ï P La 
(5) 49) 


= c, = |40 f,(6)Cos(2n0) 
0 


Connu = 
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On peut également prendre en compte des conditions aux limites impaires de la forme : fe(rt -Ө)=- 
fe(0) et fo(-0)=-fo(0) transférant cette propriété sur la solution Т(г,Ө). 


Pour que la propriété impaire T(r r -0)--T(r,O) soit effective , reportons là dans l'expression de la 
fonction: 

n-0— 

e+ f In(r)=-e- f lIn(r) Vr=ez=f=0 

n»0— 

a, Cos(n0) * b, Sin(n0) = —a,(Cos(n(z —0))- b,Sin(n(z —0))) 

— a, Cos(n0) + b, Sin(n0)  -a, (C1) Cos(n0)-b, (-1)" Sin(n0)) 

=> a,Cos(n0)+b, Sin(n0) = —(- 1) (a,Cos(n0) — b,Sin(n0)) 

—a,--(-lfa, b,-(-1yb 


=> а, =Oet b, , =OVet du *£0etb,, +0 


2n 


sT n) Y (c, rn +d, ee XCos((2n +1)0)+ Y (c. r?" +d,„r ?" )Sin(2no) 


n=0,co n=1,00 


2n+1 
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II vient pour une condition aux limites impaire par rapport à l'axe des ordonnées : 
Fonction limite  f,(x-0)=-f,(0) 


T, (r,0)= SS l Pd ps Jcos((2n +1)9) + > (c. r?" +d, r” Jsin(2n0) 


n=0,0 n-l,oo 


2n —2n 
ee + Hm = 0 о Jl {7 | 
rl rl 


2n+1 —(2n+1) 
2n+1 —(2л+1) r r 
Cond + dala m 0 => Œ Солы 5 É p 
rl rl 
condition r = | , => 


(8) 42) | јао f, (0)Sin(2n0) 


ri 


n>0 


condition r = 1, => 


І, = Í 10 f, (0)Sin(2n0) = | dO f (8)Sin(2n0) + Í dO f,(0)Sin(2n0) 


1, = | 40 f,(0)Sin(2n0) - fao f(x -0")Sin(2n(x -0')) = 0' = x -0 
L, 


Il 


ja 
fao F (0)Sin(2n0) — fao f,(0')Sin(2n0") = fao F (0)Sin2n0)— fao f, -0)Sin(2n0) 
—Casml f,(0)= f,(-0)= 0 

= Casn2 f,(0)=-f,(-0)= I, = »[do f, (8)Sin(2n0) 


2п+1 -(2n«1) эл 
ell 42) |- [do f,(6)Cos((2n + 1) 


І I 


rl rl 


1,= fao 7,00)С0»((2п+1)0)= fao f,(@)Cos((2n +1)0)- fao f(x -0')Cos((2n+1)(x — 0") 


L fao É (0)Cos((2n+1)0)- [ав f, (0')Cos((2n+1)0') = fao fa(@)Cos((2n + 1)0)+ | dO f, (-0)Cos((2n+1)0. 


=> Cas n*1 f,(0)= f,(-0) > I, = jj d0 f,(0)Cos((2n--1)8) 


—Casm2 f,0)--f,(-0)2 I, =0 
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En conséquence la solution s'écrit suivant les deux sous-cas : 


fO (0)=-f (x -0) fo (0) = fo (-9) 


2 | dO ne : 
= [40 f,(6)Cos((2n +19) 


Const 5 | 2п+1 | —(2п+1) С опт 
ame. 2153 | 
La La 
| 2п+1 —(2п+1) 
2 La La 
T, x. (1>0) == > Cjan Ze. чола TV N nÀ 
T n=0,0 La E 12 
La la 


Jal) = —fo(r -0) (0) =— fa C0) 


2] dO f,(0)Sin(2n 0) E 
a = с, = [40 f,(0)Sin(2n0) 


A II Sin(2n 0) 
III 
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Dans le cas des problèmes symétriques sur l'axe radial cela donne, avec les mêmes deux sous-cas : 


10 ôT(r,Q) 1 OT(r,0) OT(r,0) 1ÔT(r,0) 1 O’T(r,0) 
r + = + + 
r òr or r^ og? or? r © r: og 
rz = f,(0) tq f,(8)- —fy(x -0) 


Т(ғ,0)|_, 5 = 


=0 


f4(0) =-fo-0) f,(8)-— f,(-0) 
2 Í dO f,(0)Cos((2n Um А 


Са 7 — um 63,47 | dO f, (0)cos((2n--10) 
1, l, 


2 
T, x. (7,0) = > Con 


r2 
Cos((2n +1 
Ws ^ "m Ges T -(2n41) ( y) 
Ls L5 


fo (9)=-fo (x -0) fo (0) =— fa C0) 


2 i d0 f, ee 


с = „=| dO f,(0)Sin(2n0) 


ат “ig 38 
У 


Sin(2n 0) 
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Enfin si la fonction limite ne possèdent aucune symétrie, on peut toujours construire une 
décomposition en partie paire et impaire suivant l'axe des x et paire et impaire suivant l'axe des y 
et impaire : 


Jo(0) = Jo x (0)+ Lors x- (9)+ foy- y. (9)+ fs y. X. (0) sachantque },(0)= f,(27 +0) 
Jo(0)+ fo(x —0)+ f,(-9) + /,00-л) 


fes 
fax 0) се =. 
faya (0) = BONE O+ CO) HO п) 
fosa (09) 0 0) f CO) f, (8x) 


On vérifie aisément que les fonctions construites vérifient bien les propriétés de symétrie : 
Jan, (0) = Soy. x. (mT -0) et Joys (0) = foy x, (70) 

ЛДу+,х-(Ө)= Joy. x. (Œ -0) et foy, x. (0) = foy. x. 70) 

foy x. (O) = — foy x. (F -0) et foy x (9) = foy. x, C09) 

Soy-,x-(O)=-Soy. x. (R -0) et Soy. (9) = ух C0) 
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Et la solution s'écrit : 

10 ÔT(r,0) 1 OT(r,0) ô*T(r,Q) 1l0T(r90) 1 &T(r,0) 
r + = + + = 

r Or Or r: og? or? r Or r^ og? 

T(r,9) , =0 


TCO) = f, (0) 


f,(0) = foy. x. (0) Jar. (0)+ Joy- x+ (0)+ Jan х (0) sachantque  f,(0)= f,(2z +0) 


fo Í (0) = POHLED CY 0л) 
foa (0) LID fn 2 Al 0)- f (0-7) 


0 


у (9) = LE) ur 70) CO 00 -л) 


4 
far 0) RO HE D 1,70) f 07) 
EE j^ ИВ" ja foy, (0)Cos(2n0) cy, х_ = fao fy x-(O)Sin((2n +1)0) 
Сузу, = fao for „x. Cos((2n+1)9) с, х = fao for „x. (0) Sin(2n0) 


2n -2n 
r nes 
fraxi di 2 (E) E | 
T(r,0) = I : + 2 Cy, x+ 2n -2n 
T nl) с. la) _|Ьь 
La l l 


Cos(2n0) 


Sin(2n0) 


2 ri 
n+ X- n -n 
e) œ) ` (Gt 
rl rl l, l 
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Le développement pour la solution en C.L. symétrique radiale : 

10 ÔT(r,0) 1 OT(r,0) ô*T(r,Q) 1ôT(r,0) 1 &T(r,0) 
r + = + + = 

r òr Or r^ op др? r дг r^ og 

T(r,0),. = f,(0) 


F0), = 0 


700) = foy x (0)+ fo. r (0) Jan x. (9) + foy. (0), sachantque f,(0)= f, (2x +0) 


sag = 100 fe Ai OU. 9(-0)+ Lë —z) 


0 


f (0) _ OR 0 Al 0)- f, (0-7) 
f, Á (@0y= 80—80) Nn 4-0) - 6-7) 
far 0) RO HE D 1,70) f 07) 


2n -2n 
r BE 
( I, | É | | 
+ Cy, x+ 2n SES 
Zx di Л пә mur 
La [> L, 


Cos(2n0) 


Sin(2n0) 


| 
D 
3 
d 
Sun == 
Q 
© 
= 
N 
S 
+ 
= 
S 
— 
+ 
| 
су 
T 
i 
LEES 
= 


r2 
DRE 
1, I, 
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Exemple : Disque plein (r,O) soumis à des conditions aux limites de Dirichlet en r 


Voilà encore un probléme avec des conditions aux limites présentant une symétrie par rapport à 
l'axe cartésien des ordonnées. Sans restreindre le probléme on peut fixer le rayon du disque à 1, 
l'angle Ө variant de 0 à 2 л. 

10 " OT (r,0) " 1 e?T(r,0) E e^T(r,0) y 1 OT(r,0) d 1 OT(r,0) = 


0 
r òr Or r 00? or? r дг r^ 00 
T(r,0)| _, fini 
1 Өє[0,л] 
T(r,0) = 
"1 10 Өє[л,2л] 
1 del, 
posons С,(0) = elo] 
0 elei) 


La solution s'écrit alors selon la symétrie recherchée (voir exemple 23): 
T(r,0)=e+ f In(r) + У; (c. r?" +d,,r ?" \Cos(2n0) + > e, CH CT ат" )Sin((2n + 10) 


n-l,oo n=0,0 


Tout d'abord la condition de finitude en O implique que f =d,,=d4,, = 0 : 


T(r,0)=e+ У с," Cos(2n0) 4 Y Car Sin((2n +1)0) 


n-l,oo n=0,0 


Si maintenant on applique la condition aux limites en г=1, cela donne 
e+ Y. c, CosQn0)* Vc, ,Sin((2n41))- C,(0) 


n=1,00 n=0,00 
En multipliant par l'une des fonctions propres la condition aux limites en r=1 et en intégrant sur 
[0,2л], il vient : 
n-0 


2л л 1 
2ел = (C,(0)d0 = [403 > e= — 
ел | (0) | T e 2 
n>0 


2л л 
лс, = Í C,(0)Cos(2n0)d0 = Í Cos(2n0)d0 =0 > c,, =0 
0 0 


T Gau = Í C,(0)Sin((2n +1)0)40 = | Sin(@n +1)9)d9 


É Cos((2n + 1)0) i » (1- Cos((2n + 1}r)) _ (1- Cos(x)) PE 

-| (2n+1) E 2n+1 = 2n4l 2n+l 
DES 

TURIS zx (2n 41) 

La solution est donc 


1 2 — r"sin((2n41)) 
T = 
36) 2 x 2, (2n +1) 
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Pour le problème 
T(r,0),  , fini 


_]0 Өє[0,л] 1 2 — r?" Sin((2n +1)0) 
ТӨ = | 0 ef[r,2x] T (r.0) = 2 > (2n+1) 


on a la solution suivante : 20:6 
Par principe de superposition pour le problème 


T(r,0)|. , fini 
T, Өє([0, 

Hej = 0 €[0, 7] 

"1 in Өє[л,2л] 
La solution est 
T(r,9) =T, (7,0) * T, Т,(7,0) 

2n+1 Ç; 2n+1 ©: 
-n 1,2 y^ Sin((2n +1)0) ат 1 2 yi Sin((2n +1)0) 
2. жу (2n +1) 2 Z s (2n+1) 
T, +T, 2 r?" Sin((2n +1)0) 
T(r,0 T, — T, 
(r,0)= +(T, ER Сағ) 

Pour les catégories plus générales de problèmes : 
1 8 , T(r.,9) | 1 eT(r, 0) _ eT(r, 0), 1 ôT(r.9) | 1 Ə°T(r,0) 


=0 
r Or Or r^ o9? or? 1 or r^ 00? 
T(r,0)|. , fini 
T(r,O), , = f(0) 


f, (0) = foy. x. (8) + foy.x- (8) Jox- x+ (8) + Jazz (0), sachantque  f,(0)= f,(2z +0) 


DEOR 1,00) * fo (x -0) CO) XL => 
DOE 7,00) + f, (x — 0) LA 0)-f,(6-7) 

DEOR fo (0) fe zt _ COS (9л) 

E LO- LE-O- 49-0) 

fro = аө foy (0) Crx = fao for.x:(0)Cos(2n0) cy, y = fao Sors x (0)Sin((2n+1)0) 
Cy. ys = [аө for x. Cos((2n+1)9) сх. - аө foy x. (0) Sin(2n0) 


n=1,00 n=0,00 


2n 2п+1 
T(r,0)= Js у 2 У сү. 06 бану” У Én Sin((2n 4 10) 
T T z 
2n+1 2n 
É > ж 06 Соѕ((2п +1)0) y2 buc <) Sin(2n0) 
T 


Л n=0,œ0 n=1,00 
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Pour le cas où les conditions aux limites sont impaires par rapport à l'axe des x, la valeur de la 
solution à l'angle л est nulle, il suffit donc de prendre en considération les conditions aux limites de 
symétrie : Ү+,Х- et Ү-,Х- qui sont les deux conditions aux limites pour lesquelles la solution 
s'annulent sur la base х=0, soit : 


AT(r,0)-0 
AT(r,0) = 0 кра 
Т(г,Ө)]„Лпї Siro- 16) 
T(r,0),., = f,(0) avec f,(8)  — f5(-0) =, 


T(r,0), , =T(r,0),. =0 
f,(0) = foysx- (0) Sox-x- (0) sachantque  f,(0)= f,(27 +0) 
SE (6) = 18) 0) - f,(-0)- 7,00 T) _ ],00)+ f,(x — 0) 


4 2 
fag s (0) LOG 0) L Аа оа мө мл 


Cy, = [40 foy, x-(O)Sin((2n+1)0) c, = [40 fay x.(0) Sin2n0) 
0 0 


iega У 2 Sin(2n 19). 2 > «d Sin(2n0) 


T n=0,00 n=1,0 


Cela permet de calculer dans le même temps la solution du problème aux limites sur un demi 
disque pour lequel la base plate est portée à la valeur nulle. 
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Exemple : Disque plein (r,0) soumis à des conditions aux limites de Robin enr 


Voilà encore un problème avec des conditions aux limites présentant une symétrie par rapport à 
l'axe cartésien des ordonnées. Sans restreindre le problème on peut fixer le rayon intérieur du 
disque à 1 et le rayon extérieur à k, l'angle Ө variant de O à 2 л. 


2 
18 óT(n8), 1 076,6) _ 


0 
r òr Or r^ op 
T(r,0)|. , fini 
| 1 Өє[0,л] 1 Өє[0,л] 
T'(r,0)+hT(r,0)| = роѕопѕ С,(0) = 
= (0 Өє[2,2л] 0 Өє[л,2л] 


La solution s'écrit alors selon la symétrie recherchée (voir exemples 23, 24 et 25): 
T(r,0)=e+ f In(r) + > (c, +d,,r °" Jcos(2n0) t > (c, nn +а, ат" Jsin((2n t 10) 


n=1,00 n=0,0 


= d 


2n+1 


Tout d'abord la condition de finitude en O implique que =й, =0 : 


T(r,0)=e+ У: с," Cos(2n0) KN Sin((2n +1)0) 


n=1,00 n=0,0 
Si maintenant on applique la condition aux limites en г=1, cela donne 
> c,,2nCos(2n0) + > c, (Cn+1)Sin((2n+1)0)+he+h D c,, Cos(2n0) - h > c, Sin((2n+1)9)=C,(0) 


n=1,00 n-0,oo n=1,00 n=0,00 

En multipliant par l'une des fonctions propres la condition aux limites en r=1 et en intégrant sur 
[0,2л], il vient : 

n-0 


2л л 
1 
2her = | C;(0)d0 = | dO =z => e = — 
| 100) | T 
n>0 


2л л 
л (2n c,, +hc,,)= | C;(0)Cos(2n0)d@ = | Cos(2n0)d0 = 0 = c,, =0 
0 0 


л ((2n 1)c,, + he, )= c,  ((2n+1)+h)= Í C,(0)Sin((2n -1)8)a0 = j Sin((2n +1)0)40 


Е | Cos((2n + 2 _ (1 ú Cos((2n + 1x )) (1 _ Cos(x)) 2 


(2л +1) E 
et 2 
2 z (2n Y((2n 1) 5) 


_ 2 r?" Sin((2n +1)0) 
Ө Z Z Ga n D о) 


2n+1 2n+1 = 2n+1 


La solution est donc : 
Pour le problème : 


2 
10 TC) 1 TCA) 


= 0 
r Or Or г? 00° 
T(r,0)| _, fini 
0 Oe, 0 Өє[0, 
T'(r,0)+hT(r,0)| = SE posons С,(0) = SE 
ri |1 Өє[л,2л] 1 del Ze) 
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T. (r,0) = 1 2 Y r?" Sin((2n-- 1)0) 
on a la solution suivante : 2h m Ze (2n +1)(2n+1) e n) 
Par principe de superposition pour le probléme : 
2 
1 Oo ,T(r.8) , А д Fe) E 
r Or Or r 00 
T(r,O)| _, fini 


| hT, Gelz 
T'(r,0)+hT(r,0) _ zx hT. Oe [л 2л | 
1 H 


La solution est 
T(r,0) = РТ, Ti (r,0O) - hT, T)(r,0) 


Е zl Ë 2 y E, У ( NES, y E, 


2h x,%.(2n+1\(2n+1)+h 2h л„&(Оп+1\(2л+1)+л) 


Eat _ m y2h r?" Sin((2n+1)0) 
RES Ee P CE 


Pour le probléme général : 

2 
1 ð „T(r.0) | 1 OT(r,0) ET 
r òr Or r^ 00? 
T(r,0)| , fini 
T.(r,O)* hT(r,8) = f,(0) 


La prise en compte des conditions aux limites mixtes donnent pour une fonction propre soit de 


valeur propre nulle, soit non nulle : 
n=0 R(r)-1 


T.(r,0)+ RT(r,0) , > Coefficient diviseur h 
n>0 R(r)-r^ 


T.(r 0) hT(r,0) | => nl," + hl," = É + 72 => Coefficient diviseur É + ) 


r r 
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On reprend les mêmes termes d'analyse des conditions de symétrie des fonctions aux limites : 
1 0 ,.0T(r,0), 1 O°T(r,0) 9) _ e'T(r,9) | 1 ôT(r.,9) | 1 O°T(r,0) 
r òr Or Doo др? g= Or "e EFT ` 
T(r,09)| | fini 


T,(r,0)+hT(r,0),, = 7,09) 


700) = foy x (0) + ух (0) + foy y (9) * foy. (0), sachantque  f,(0)= f,(27 +0) 
fars. xa (0) = f,(0)+ f, (z I Ja (0 =z) 
fo(0)+ fox -0)- f,(-0) - f (0 — x) 
4 
fau (0) = EP C 0л) 
foy a (0) LOAA D Ө)+ f, (87) 


-0 


Jar, (9) = 


= [аө foy x Cos((2n+1)0) c, х. -Í dO fay x, (0) Sin(2n0) 


0 


2n 2n+1 
тө) 4.2. Уу E У жос mn 
ha ЫЕ) (ан, i l, 
l, l 


2 >. SEN 8 Cos((2n + 1)9)+ 2 >. ee 8 Sin(2n0) 
Ë D | 


л n=1,00 + ) r 


r 
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Avec une condition aux limites mixtes de forme encore plus générale, on obtiendrait : 


10 ÔT(r,0) 1 ƏT(r,0) OT(r,0) 10T(r0) 1 &T(r,0) 
r + = + + = 

r Or or r^ 00? or? r Or r^ op? 

T(r,O)| | fini 

aT, (r.0)- BT(,O)J | = 7,00) 

SaO) = fo у. x (0) Jan, (0) foy x (0) + foy (IL sachantque f, (0) — f,(2x +0) 


0 


ТД) E LEE) S 
pe г (s) 

Ру (0) = fy foie - 0) IED) - f.(0 - x) 

0) ROAD FEES 

frag. = fao Sov+x+(0) Crx = fao fay x+ (0)Cos(2n0) cy, x = fao Sors x -(0)Sin((2n+1)0) 
es m [40 fos x. Cosa). o, = [49 Jay x Ө) Sin(2n0) 


2n 2n+1 
T(r,0) = fux |2 s 8 Cos OD) S Lud B Sin((2n +1)0)+ 
px T n=1,00 ae L, n=0,00 a be L, 
І 


r r 


+В 


(a 


2n+1 2n 
2 Cy. y. r 2 Cy. x- r š 
; Cos((2n +1 : Sin(2n0 
За = jo os((2n + КЫ 2 d in(2n0) 
o а 
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Exemple : Disque creux (r,O) soumis à des conditions aux limites de Dirichlet en r 


Sans restreindre le probléme on peut fixer le rayon intérieur du disque à 1 et le rayon extérieur à k, 
l'angle Ө variant de 0 à 2 л. 


2 
10 , T(r.0) | E ô 00) 20 
r Or Or r 00 
T(r,6),, =T, 


(7.0), =T, 


La solution est à rechercher avec la valeur nulle de l'équation séparée : 


di 
La 
La 


Et pour le problème général suivant : 
10 ÔT(r,0) 1 @T(r,0) 

r + 2 2 = 
r òr Or r 00 


T(r,0),, = fi (0) 
r=l„, zi Jao (0) 


La solution a été donnée en exemple 23, que l'on compléte par principe de superposition : 
f, (0) = Луз, х+ (0) Jior x- (0) + Jar z- (0) + Jar x- (0) sachantque  f,(0)= f,,(27 +0) 


fig (0)+ fi 9 (T —0)+ f19(C0) + fi (0 — z) fig (0)+ f (r - 8) - fi (C9) - fi (0 — z) 


T, (1.0) =T, +(T, -T,) 


0 


T(r,0) 


Joy x (0) = 4 Join, (0) = j 
far) ell f (z —8) +) -/4@-л) ә (ду fa Pate - 0 SE В) 


Љо (O) = foo y. x (9) + oo y x (0) + foy x (0) + оу х (0), ѕасһапідие },,(0) = f, (2z +0) 
J (0)+ f; (c - 0) f 70) * f. (0 —z) Jf (0)+ f; (z —0)— f/2,(-89) - f. (0 —z) 


збы (9) = 1 Poria (0)= ^ 
ut (9) = J10(0) - f(x - 0) 07 — f. (0 - x) Хат XO) Е f;4 (0) - f. (z — 0) H (-0)+ f., (0 — x) 


Хох [АӨ Kan, O): zcy |40 fios (0)Соз(2п0) с.х = [40 fios. s (O)Sin((2n 10) 
0 0 


0 
yx. = | dO for x, Cos((2n+1)0) с, = |40 fior x, (8) Sin(2n0) 
0 0 


Dya x+ = dëi fr: x (0)Sin((2n +1)0) 


40 er: x: (9) Co y+,xX+ T Í dO Zaart (8)Cos(2n0) C y4 x- — 
0 


о —Ə à 


DS c — 


C» y- X+ 


dO fo, x.Cos((2n+1)0) с,у_ = [40 уу (0) Sin(2n0) 
0 
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r 
Tirri m) 2 
T(r,0)= Bed us 


E > 2n -2n 
л «| l, ) др. La E 14: 
2 1, L5 


< Sin(2n0) + 


) 
Pays d | sl y 
` ue — 
I QT 
Ё; т | Ше» 
б 


EZ OK. 8). 
< (ey Ge" = III 


Sin(2n0) 
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Tableau synoptique des fonctions propres suivant la symétrie des conditions aux limites 


Propriétés de symétrie Symbole 
symétrie 


de 


Fonctions propres utilisées dans 
le cas d'un disque complet en 
angle creux ou plein 


Jo(0)= f(x -0) et f,(0)= f,(-0) Y+ X+ 


Cos(2n0) 


J(0)= 107-0) et f,(0)-—f,(CO) — Ys,X- 


Sin((2n +1)0) 


f,(0)=-f,G-0) et f(0)=f,C0) — |Y-X+ 


Cos((2n+1)0) 


f9(0)- —fo(r-0) et 7,00) = – 36-0) Y-,X- 


Sin(2n0) 
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Exemple : Solution de l'équation de Laplace à deux dimensions sur des corps de dimension infini, 
soit en coordonnées cylindrique (г, 9), la solution de l'équation de Laplace sur un secteur d'espace 
partant du rayon Ir avec des conditions aux limites de Dirichlet. 


AT (r,0) 
(r0). -0 
T(r.9),, = 0 
T(r.), , = 7,0) 
T(r,0) fini 


La solution se développe en série de fonctions propres sous la forme : 


T(r,0)=e+ fin(r)+ Sfr + gr^ Xe, Sin(4,8) * d,Cos(A,0)) 


n-l,oo 
nm 


CL.=d,=0 f,=0 À,= 
Ө, -0, 


c, = Í dx f,(0) Sin(A,0) 


2 сї” 
T(r,0)= > 22—Sin(A,(6 -0,)) 


| 1- (-1) (2n+1)x 4 — L^ Sin(2,(0-0,)) 
si f,(0)=1= с, = = À = T(r,0)= > r n 1 
»@®) Â 7 "^  0,-0, E gt A 2n+1 


Pour le problème 
AT(r,0) 
T(rj8), = 0 
T0). =0 


T(r,O), = f,(0) 
T(r,0) fini 
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La solution se développe en série de fonctions propres sous la forme : 


(2n+1)r f | 
„=> с, =| dx f,(0) Sin(2,0) 
2(0, SCH | : 
dope el Sin(A, (0 — 0.) 
2 0, — 0, k y^ n 1 
| 1 (2n +1)л 4 — L^ Sin(A,(0 - 0 ) 
si f,(O)=1=c, =— > À, ==> T(r,0)= r n f 
f, (0) pe eue = 2: 220 
et pour le probléme 
AT(r,0) 
Т,(",0) = 
T(r,0),., = 
T(,0),, = 7,00) 
T(r,0) fini 


La solution se développe en série de fonctions propres sous la forme : 


E (2n+1)r Е ü 

= "H SE | dx f,,(0)Cos(A,0) 
T(r,0) = 22. d Cos(1, (0 — 0,)) 

| Е _ (=D" _ (2n+1}r m 24 (0-0) 
si f,(0)=1— c, = S ы к: T(r,0)= => De FR 
et pour le problème 
AT(r,0) 
Tro) = 
T,(r.9), , = 
T(r,9)| , = 300) 
T(r,0) fini 


La solution se développe en série de fonctions propres sous la forme : 


An nT 


moe T NE ы," Cos(2,(0 — Ө, )) 
0,-0, 0,-0 т” 
si f,(0)=1=c„=0= T(r,0)-1 


e | 10 f,(0)Cos(4,0) F= Í dx f,(0) 
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Exemple : Soit en coordonnées cylindrique (r,8), la solution de l'équation de Laplace sur un secteur 
d'espace partant de l'origine avec des conditions aux limites de Dirichlet. 


AT(r,0) 
T(,0), , = 0 
T(r.9),, ы, =] 
T(r,0) fini 


La solution ne dépend pas du rayon, les solutions sont donc à rechercher dans la classe des 
fonctions linéaire de l'angle 


T(r,0)=a0+b=T(r,0)= SE 
Ө, —0, 


Lorsque l'angle est plat, à l'horizontale alors (demi-plan supérieur ): 


ааа? 
— NA эу 0,y 20 


6-0 62z T(rO)=Z=ST(x;y)=| 7 
С ТОН 
|r =Â Ê = 267; six<0,y 20 
л 
Lorsque l'angle est plat, verticalement, alors (demi-plan gauche ): 
1 
9+7 FA 5+2 sio elo Z 
бузы uw ep- ЖШН God à 
2 2 л 2: T 1 80-27 3 0 3л 
=——+ si 0 ‚2л 
2 л 2 m 2 
0 
x š Y 
+ si x> 0, y > 0 | ОН 
п À =>+ х 
3 2л + And Z | ааа ý 
=e VEZA q: si x > 0, y «0 
2 T 2 л 


Lorsque l'angle est droit : 


0 Aretan >) 
T x 
0, 20 шс e >Á 


